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Abstract 



One of the least accurate tested sectors of the Standard Model of elementary particles is 
the quark mixing mechanism. It encodes the transitions between different quark families in 
weak charged-current decays in terms of the Cabibbo-Kobayashi-Maskawa-matrix Vckm- 
A complex phase S in this matrix is the only possible source of the violation of the CP- 
symmetry in the Standard Model. This fact makes the study of Vckm mandatory: As 
soon as the measurements of CP-violation cannot be fitted with the single parameter 6, 
evidence of new laws of physics will be found. Today the only unambiguously measured 
CP-violating phenomenon is the indirect CP-violation present in the |AS| = 2 transitions 
inducing the mixing of the neutral Kaon states K'^ and K^. It is characterized by the 
well determined parameter sk- Even from this single experiment on CP-violation one can 
test the Standard Model, because together with the unitarity of Vckm it constrains some 
parameters of the Standard Model such as the mass of the top-quark. 
On the theoretical side the study of Vckm is challenged by the presence of the strong 
interaction, which may screen or enhance the weak transition amplitudes of quarks and 
which binds the quarks into hadrons. The short distance QCD effects are comprised in the 
parameters i]i,ri2 and 773 of the effective |AS| = 2 -hamiltonian: 
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b{ii)Qs2ilJ') +i^-c.. 



For the following outline we only have to explain that = VjfiV*^ encodes the relevant 
elements of Vckm and that rijb{fi) = 1 in the absence of QCD corrections. 
Now for the investigation of ck the coefficient 772 is most important, but a precise deter- 
mination of the subdominant term involving 773 is also required for a reliable analysis. The 
mass difference between the two physical neutral K-mcsons Kl and Ks, however, is dom- 
inated by r]i. When this work was begun, only the coefficient 772 was known beyond the 
leading log approximation [15]. This approximation is unsatisfactory, as it involves large 
errors due to renormalization scale ambiguities and as it does not allow for the use of the 
fundamental QCD scale parameter Aj^. 
The results of this thesis include the following: 

i) The calculation of 771 and 773 in the next-to-leading order (NLO) is presented in 
detail. 



ii) An improved determination of ImAt, which is the key quantity for CP-violation, 
using the measurement of Ek and the NLO values for 771, 772 and 773 is given. Results 



are tabulated as functions of the CKM-parameters Vcb and Kb/ Kb, the mass of the 
top-quark and of B^, which parametrizes the size of the hadronic matrix element. 

The bounds on the Standard Model parameters resulting from the unitarity of Vckm 
are studied. We have found that with the new values for rji and 773 a wider range of the 
parameters mentioned in ii) are allowed. In this sense the calculation has vindicated 
the Standard Model. 

The NLO value of 771 leads to an enhancement of the short distance contribution 

to the Kl — Kg-mass difference. Now the short distance QCD contributions explain 
the dominant part of the measured mass difference. The old leading order result has 
suggested a long-distance dominance. This has been a puzzle, because long-distance 
effects are expected to be suppressed by a factor of A.qcb/''^c with respect to the 
short distance contributions. 

A new feature of the NLO calculation of rji and rj^ is the presence of Green's functions with 
two insertions of operators (bifocal structures) both of which have a non-zero anomalous 
dimension. These Green's functions have required a subtle analysis of the renormalization 
of effective field theories. The discussion includes the following points: 

v) The RG of bifocal structures is analyzed in detail. Different methods for the solu- 
tion are presented. These include the framework of an inhomogeneous RG equation, 
which is best suited for formal proofs like those mentioned in vii). For the practical 
calculation a matrix formulation is given which in the case of rj^ only involves one 
8 X 8-matrix instead of four of them as in the old LO calculation in [45]. 

vi) The structure of the reduced effective |AS| = 1 lagrangian is discussed. In matrix 
elements with single operator insertions those operators which vanish by the field 
equation of motion or which are BRS-exact do not contribute to the matrix ele- 
ments. But one faces non-trivial contact terms, if the Green's function contains two 
operator insertions. It is well-known [52,54,58] that these terms can be absorbed into 
I AS| — 2 operators. In our case the contact terms have been identified to correspond 
to operators describing subleading effects in m^/mc. 

vii) Likewise |AS| = 1 evanescent operators appear in the effective lagrangian. Here we 

had to develop the correct prescription to deal with evanescent operators in bilocal 
structures. Just as in the case of the operators mentioned in vi) non-vanishing terms 
appear, which are also equivalent to the matrix elements of |AS| = 2 operators. But 
now they can be absorbed into a finite renormalization of the |AS| =2 operators. 

viii) The analysis of vii) has stimulated a related investigation: The definition of evanescent 
operators is not unique, because one can add {D — 4) times any physical operator 
to them. We show that any definition of them leads to an effective lagrangian, in 
which the physical sector is unaffected by the presence of properly renormalized 
evanescent operators. Hence one can drop the latter from the reduced lagrangian. 



iii 



Yet the arbitariness in the definition leads to a renormahzation scheme dependence 
in the physical sector. This result is important for practical NLO calculations, because 
the Wilson coefficients and the anomalous dimension matrix depend on the definition 
of the evanescents. Formulae to transform between the schemes are presented. 

ix) Last but not least there are about 40 two-loop tensor integrals up to rank six involved 
in the calculations. The calculation of rji has required to calculate their finite part. 
Here wc present a compact formula for the most general integral which can appear in 
two-loop calculations of light hadron systems. This is the two-loop vacuum bubble 
diagram with arbitrary tensor structure, small infrared regulating masses on the lines 
flown through by one loop momentum and arbitrary heavy masses on the other lines. 
The formula holds in D dimensions and is therefore a good investment into the future 
of higher order QCD corrections. 
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1. Introduction 



The aim of elementary particle physics is the discovery of new fundamental laws of na- 
ture and their classification within a theoretical framework which arranges the observed 
elementary particles into simple patterns and describes their interactions in a unified way. 

The stage for this theoretical description is the Rclativistic Quantum Field Theory incor- 
porating both Einstein's theory of Special Relativity and Heiscnbcrg's, Pauli's and Schro- 
dinger's Quantum Mechanics. Its story of success started in 1947, when the experiment 
of Lamb and Retherford [1] gave unambigous evidence for the splitting of the 2si/2 and 
the 2pi/2 energy levels of the hydrogen atom by 4 //eV. At the same time Dirac, Feynman, 
Tomonaga, Schwinger and others [2] developed the theory of Quantum Electrodynamics, 
which explains the Lamb shift by quantum fluctuations of the photon fleld. Since then 
many more actors have joined the electron and the photon fleld on the stage: Starting with 
the works of Gell-Mann and Zweig the partonic constituents of hadrons were successfully 
described by quark fields, whose strong interaction is mediated by gluons as encoded in 
the Lagrangian of Quantum Chromodynamics [3]. And with the unification of the electro- 
magnetic and weak interaction by Glashow, Salam and Weinberg [4] the scene was entered 
by three more vector bosons, which are massive owing to their interaction with the Higgs 
particle [5]. These elements form what is now called the Standard Model of elementary 
particle physics. It is the most successful theory of physics: It covers energy scales between 
those of the Lamb shift and the center of mass energy of today's particle colliders such 
as LEP at CERN, which are separated by 17 orders of magnitude! No experiment has yet 
contradicted its quantitative predictions, although many of them have been tested to an 
extraordinary precision (e.g. an accuracy of 10^^° for the anomalous magnetic moment of 
the electron on both the experimental and the theoretical side) . 

On the other hand some sectors of the standard model are poorly verified at present. One of 

them is the quark mixing sector, which describes the coupling of quark fields to the charged 
electroweak vector bosons in terms of the Cabibbo-Kobayashi-Maskawa (CKM) matrix [6]. 
Here the non-perturbative nature of the long-distance strong interaction prevents an easy 
extraction of the CKM parameters. 

Further we know that the Standard Model cannot be true for another 17 orders of mag- 
nitude in energy, as we then would reach the Planck scale, where gravitation becomes 
important requiring a new theoretical concept. Yet we know from other considerations 
such as vacuum stability that there must be new physics much below the Planck scale. 
How can we detect new physics apart from building accelerators with higher and higher 
energy? 

First one can look at rare processes which are forbidden in the tree level approximation. 
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Since they are suppressed in the standard model new physics can contribute sizeably to 
their transition rates. Moreover, even low-energy rare hadronic processes are sensitive to 

the masses of heavy particles appearing in intermediate states. This allows to derive bounds 
on the mass of the recently discovered top-quark [7] from these processes [8,9]. When the 
experimental determination of the top mass improves, these bounds provide a consistency 
check of the Standard Model quark sector over more than two orders of magnitude in 
energy. 

Second one can look for violations of symmetries which are respected or weakly broken by 
the Standard Model. Here the CP symmetry is of utmost importance, because the only 
possible source of its violation in the Standard Model is a single parameter, the phase 5 
in the CKM matrix.^ The CP transformation exchanges particles and antiparticles and 
reflects the spacial coordinates. If in the future we cannot fit all CP violating observables 
with this single phase, we will have discovered new physics. Further in the calculation of 
these observables the difficulties with the long-distance strong interaction are sometimes 
reduced, as the strong (as well as the electromagnetic) interaction respects CP symmetry 
indicated by the smallness of the neutron electric dipole moment. 

All these aspects make the study of hadronic physics mandatory. Here a key object is the 
K-meson system, to which we owe a significant part of our present day knowledge about 
the Standard Model quark sector: In 1964 Christenson, Cronin, Fitch and Turlay have 
found the first evidence for CP violation by observing decays of the long-lived neutral 
K-meson Kl into the CP even two-pion state [10]. There are two potential sources for this 
phenomenon: The Kl state is not a pure CP odd state but has admixtures of the CP even 
neutral Kaon state or the weak interaction triggering this decay violates the CP symmetry 
thereby allowing transitions from the CP odd neutral Kaon state to CP even final states. 
In the first case we speak of indirect CP violation in contrast to direct CP violation in the 
second case. The CKM mechanism of the Standard Model predicts both types, if the phase 
S is not equal to or tt. While precise measurements are available for the parameter ek [H] 
characterizing the size of the indirect CP violation in the neutral Kaon system, there is still 
a controversy about the direct CP violation [11]. Then in 1970 the suppression of flavour- 
changing neutral current (FCNC) decays of K-mesons lead Glashow, Iliopoulos and Maiani 
(GIM) to postulate the existence of the charm quark [12]. We finally mention the prediction 
of the charm-quark mass from the observed Kl — Ks -mass difference by Gaillard and Lee 
in 1974 [13]. On the other hand the theoretist's work to relate the measured quantities to 
the Standard model parameters is made difficult by the presence of the strong interaction: 
Its long range interaction confines the quarks into K-mesons and it modifies the weak CP 
violating amplitudes of interest. Since these transitions take place on short distance scales of 
the order of the inverse W-boson mass Mw, the QCD corrections can be reliably calculated 
perturbatively. The long-range effects, however, must be treated by other methods such 
as 1/Nc-expansion, sum rule techniques or lattice gauge theory. 

The Technical University of Munich has already largely contributed to the theoretical 
understanding of short distance QCD corrections to various rare hadronic processes [8,15, 



We do not consider a possible 6-term in the QCD lagrangian here. 
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17,41,43,44,46,49,61,62,66]. 

In this work we will focus on short distance QCD calculations to the mixing of the flavour 
eigenstates K'' and K° of the neutral K-meson. Since the mass eigenstates Kl and Ks are 
linear combinations of them, K° — K° -mixing encodes the information on the indirect CP 
violation. 

The hamiltonian describing this mixing reads: 
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Here |AS| = 2 denotes the change in the strangeness quantum number in the transition 
between K° and K°, Gp is the Fermi constant, \j — VjdV^^ comprises the CKM-factors, 
and Qs2 is the local four quark operator (see fig. 5.3 on p. 64) 

Qs2 = (Sj7^(l -75)o?i)(sfc7''(l -75)4) = {sd)v-A{sd)v-A (1-2) 

with j and k being colour indices. The Inami-Lim functions S{x) and S{x, y) [14] depend 
on the masses of the charm- and top-quark and describe the | AS| = 2-transition amplitude 
in the absence of strong interaction. They are obtained by calculating the lowest order box 
diagrams depicted in fig. 5.1 on p. 61. The short distance part is included in the coefficients 
rji, 7)2 and rj^ with a common factor 6(/i) split off. 

Here rji and 773 are most important for the size of Ek, while r/i gives the dominant contribu- 
tion to the short distance part of the Kl — Ks-mass difference. When this work was begun, 
a next-to-leading order (NLO) calculation for r)2 had been performed [15], but no QCD 
corrections beyond leading logarithms had been known for rji and 773. The leading log (LL) 
approximation is very unsatisfactory, because it is plagued with huge renormalization scale 
uncertainties and no use of the QCD scale parameter can be made. In this thesis the 
complete \AS\ = 2-hamiltonian in the NLO will be presented [67]. 

As the importance of Ek has been stressed already, let us now discuss the Kl — Ks~mass 
difference: Its measured value is 3.5 /xeV, which, by the way, is of the same order of magni- 
tude as the Lamb shift mentioned in the beginning of this chapter. The prediction of the 
charm-quark mass, mc ~ 1.5 GeV, from the Kl — Ks-mass difference in [13] was surpris- 
ingly good for an analysis which neglected strong interaction effects. The inclusion of the 
latter, however, made the theoretical prediction worse, as estimates of the hadronic matrix 
element {Qs2) by various nonperturbative methods resulted in a significant reduction of the 
naive vacuum insertion value used in [13]. Further the leading-log calculation of the short 
distance QCD corrections gave an additional suppression of the predicted Kl — Ks-mass 
difference due to r/i < 1, so that the calculation reproduced less than half of the mea- 
sured Kl — Ks-mass difference. On the other hand one also expects additive long-distance 
contributions corresponding to light meson poles which are not contained in (1.1). These 
contributions are poorly calculable, but should be suppressed with respect to the short 
distance part by power counting arguments. Therefore the low result of the short distance 



4 



1. Introduction 



calculation is somewhat surprising and in fact triggered some speculation to assign the 
deficit to new physics. Our NLO calculation for rji [17], however, together with present day 
values for the input parameters such as has established ?7i > 1 implying a sizeable 
increase of the short distance contribution to the Kl — Ks-mass difference, which is now 
roughly of the order of 70% of the observed value. This fits much better with the expec- 
tation, especially when contrasted with an earlier incorrect NLO calculation [18] finding a 
drastical decrease of the short distance part. 

In the following chapter we will summarize some field theoretic basics needed for the 
subsequent chapters. Chapter 3 will discuss the operator product expansion and renormal- 
ization group techniques. Here we will also develop the renormalization group formalism 
for Green's functions with two operator insertions. Chapter 4 contains new results about 
evanescent operators [19]. They comprise the correct treatment of bilocal structures in- 
volving evanescent operators, which had to be developed for the calculation of 773. Further 
a new scheme dependence associated with the definition of the evanescents is discussed. 
In chapter 5 the operator basis involved in the NLO |AS| = 2 lagrangian will be analyzed. 
In chapter 6 the NLO calculation of rji and rjs will be described in detail, and chapter 7 
is devoted to the phenomenological analyses. Here we will analyze the consequences of the 
new NLO prediction for ek on the Standard Model parameters. We will also look at the 
new prediction of the short distance contribution to the Kl — Ks -mass difference. 
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2. Quantum Field Theory and the 
Standard Model 

Our today's knowledge about the rich field of elementary particle physics is the result of 
many decades of research probing the laws of nature at small distances. The experimental 
effort has been paralleled by a similarly successful development of the theoretical framework 
to understand the guiding principles of the observed phenomena. 

Clearly we cannot summarize all results of this development here, so that we instead 
focus on those aspects of modern quantum field theory which are relevant for the actual 
calculations concerning the K'^ — K^ -mixing performed for this thesis. 

2.1. Lagrangian Field Theory 

A classical field theory like electrodynamics contains a set of functions {(f)j{x)} of the space- 
time variable x = {t,xy as its basic building blocks. The dynamics of these fields (j)j{x) is 
encoded in a scalar functional, the lagrangian (density) £ [(j)j{x), d^(f)j{x), x], from which the 
field equations of motions can be obtained. Here we have collectively summarized indices 
related to intrinsic properties like spin or charge quantum numbers and those distinguishing 
different physical fields into the single multi-index j. 

When passing from classical to quantum field theory the functions 4>j{x) are replaced 
by field operators 4>j{x). Hence one has to deal with a quantum mechanical system with 
an uncountable infinity of degrees of freedom labeled by the space-time coordinate x. 
Its hamiltonian density T-[{x) is a polynomial in the elementary fields <f)j{x). For non- 
interacting fields the field theory can be solved, i.e. one can construct a Hilbert space (the 
Fock space) whose elements are eigenfunctions of the Fourier-transformed hamiltonian 
density 7i{k). A short glimpse on the simple structure of the lagrangian and on the variety 
of particle phenomenology shows that it is hopeless to try to solve the interacting quantum 
field theory exactly. 

2.2. Symmetries 

Before discussing the dynamics of field theory it is useful to look for underlying patterns 
which may possibly restrict the particle spectrum and the interactions. This is indeed the 



Throughout this thesis we use natural units c = h = 1. 
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case, if the total action 

= j d^xC[(t>^{x),dMx),x\ (2.1) 

is invariant under some group of transformations of the arguments of C 
Consider first some group G of linear transformations of the space-time variable: 

x*" x'^ := with u^G (2.2) 

If we can find a (possibly multivalued) representation Vi^G) of G such that with 

- [d{u;)yeV{G) 
Mx) ^ 0;(x') = 0;(u;-x) := [d{u;)].Mx) (2.3) 

the action (2.1) satisfies S[4>'j] = S[4>j], we speak of an external symmetry. 
A symmetry group which only transforms the fields according to (2.3) but leaves the space- 
time coordinate unchanged is likewise generating an internal symmetry. 
The key importance of symmetries for a classical field theory is due to Noether's theorem: 
If a symmetry Lie group leaves S invariant, one can construct a conserved current from 
the fields and find a globally conserved charge. These charges can be used to label different 
solutions of the equations of motion. A quantized theory may or may not enjoy a certain 
symmetry of the corresponding classical theory. In the second case one terms the symmetry 
to be broken. Otherwise the operator corresponding to the conserved charge commutes with 
the hamiltonian so that we can find a common eigenbasis of both operators. Wc can use 
the quantum numbers of the conserved charge to label the eigenstates of the hamiltonian 
and transitions between states with different quantum numbers are forbidden. In fact this 
property also applies to discrete symmetries. 



2.2.1. External Symmetries 
Poincare— Invariance 

When probing the laws of nature at distance scales at which gravitation is unimportant 
they turn out to be invariant with respect to 

• translations in space and time, 

• rotations 
and 



• boosts into uniformly moving coordinate systems. 

The last two types of transformations generate the proper Lorentz group , all of them 
generate the Poincare-group. We can test the abovementioned symmetries in a twofold 



2.2 Symmetries 
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way: First we can repeat the experiments in different transformed coordinate frames^ and 
second we can check the associated conservation laws predicted by Noether's theorem. 
They include the conservation of energy, momentum and angular momentum. 
The fact that nature is invariant under (2.2) with G = greatly simpifics the theorist's 
work to guess the correct lagrangian describing the non-gravitational interactions: C must 
not explicitly depend on x and all fields must transform as in (2.3) according to some 
(possibly multivalued) representation of the proper Lorentz group. Since with (2.2) df'x — 
d^x', these fields must combine to L^-singlets in C. 

The proper Lorentz transformations leave the Minkowskian scalar product 

invariant. Yet there are more transformations which have this property: We may extend 

V by 

the parity transformation P: a; — > —x 
and 

the time reversal transformation T: xq — > —Xq. 

L {L\ PL\ TL\ PTL^} is the full Lorentz group. Unlike the transformations of we 
cannot carry out the discrete symmetries P and T by any active transformation of our 
experimental device, but still we can test them via the conservation or nonconservation 
of the associated charges. Since {1,P} and {1,T} are two-element groups, the associ- 
ated charges P and T can only assume two difi^erent values. We may choose P and T as 
multiplicative quantum numbers being +1 or —1. 

There is no reason why nature should respect these discrete symmetries. Nevertheless it 
was not before 1956 that Lee and Yang discovered the nonconservation of parity [20]. In 
the following we will need the transformation property of Dirac fields under P and T. They 
transform as 

P : ijjix) i/j'{x') = 77p7V(2;) (2.4) 

with an arbitrary phase rjp in any representation of the Dirac algebra. The current j'^(x) = 
ip{x)j^ip{x) is easily found to transform as a vector field, i.e. like x'^, under P. With 

75 := "^70717273 the axial current j"^^(x) = ip{x)j^jr^ip{x) likewise transforms with an extra 
factor of (—1) compared to j^{x). In the standard representation of the Dirac algebra one 
further finds: 

T : ^{x) ^'{x') = riTi-f^-f^ip*{x) (2.5) 
with rjT also being an arbitrary phase factor. 

Scale Transformations 

In the early 1950's Stiickelberg and Peterman as well as Gell-Mann and Low for the first 
time studied the behaviour of a quantum field theory under dilatations of the space-time 

^Theorists are allowed to neglect the technical problems, if this is tried for the boosts. 
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variable: 

X ^ x' = \-^x (2.6) 

or equivalently in momentum space 

V ^ p' = Ap. (2.7) 

One might expect that for large A masses in £ become negligible and any physical quantity 
scales according to its naive dimension determined by power counting. But this is not the 
case for a quantized theory due to the subtleties of renormahzation to be discussed in 
sect. 2.4.3. Here we find the powerlike scahng behaviour modified by so called anomalous 
dimensions. In the context of renormahzation the symmetry under scale transformations 
becomes a powerful calculational tool which will be looked at in detail in chapter 3. 



2.2.2. Internal Symmetries 
Gauge Symmetries 

In view of the discussion of Quantum Chromodynamics (QCD) wc first look at an SU(N)- 
symmctric theory containing Dirac fermions, which are chosen to transform according to 
the defining representation of SU(N): 

V'.' = C/i,(nV',-, C/ = e-^"^\ (2.8) 
Here the T^'s span the Lie algebra su(N). They satisfy 



■fabcrpc ^2.9) 



with the structure constants /"^'^. Continous internal symmetries are called gauge sym- 
metries. If the gauge parameters 9"' in (2.8) are independent of x, one speaks of a global 
symmetry. Local gauge symmetries, where the ^°'s are allowed to depend on x, play a key 
role in field theory: They incorporate a natural description of interactions between the 
fermions mediated by a vector field yl^(a;). The lagrangian 

c = ^.(i|).._^5,,)v,-^F;.F"'^'' (2.10) 

with 



is invariant under local SU(N) transformations, if the vector field ^^(a;) transforms as 



rpa j^la 



= u{e{x)) 



T^Al 



{e{x))d,u{e{x)) 



(2.11) 
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Charge Conjugation 

Charge conjugation is a discrete global internal transformation. The C operation exchanges 
particles and antiparticles. In the standard representation of the Dirac algebra a Dirac field 
transforms as 

C : ilj{x) ilj'{x') = rjcn^-f^'^'^ix) (2.12) 
with the usual arbitrary phase rjc- 

Soon after the discovery of parity non-conservation the charged weak lepton current was 
proposed to be of the type e(a;)7^(l — ■y5)i'e{x). This interaction also violates charge con- 
jugation symmetry, but it is invariant under the combined transformation CP. Finally in 
1964 it became clear that the weak charged hadronic current violates CP invariance [10]. 

2.2.3. CPT Theorem 

After quantizing a theory one wishes to find a representation of the symmetry transfor- 
mations of interest in the space of state vectors on which the field operators (f){x) act. A 
theorem of Wigner states that one can always find such a transformation, the correspond- 
ing operator is either unitary or anti-unitary and commutes with the hamiltonian. We 
have already alluded to this theorem in the discussion before sect. 2.2.1. In the case of C 
and P this operator is unitary, while it is anti-unitary for T. 

According to our present knowledge the electromagnetic and the strong interaction re- 
spect C, P and T separately, while the weak interaction violates all of them. Nevertheless 
Liiders and Pauli [21] have proven that any Poincare-invariant field theory must respect 
the combined anti-unitary transformation CPT. This has the consequences that parti- 
cles and antiparticles have the same mass and spin and that particles with zero quantum 
numbers are identical to their antiparticles. The implications of the CPT theorem are 
experimentally confirmed. 

2.3. Standard Model 

In the Standard Model the electromagnetic and the weak interaction are unified in a 
SU(2)l X U(1)y gauge invariant lagrangian £cw [4]. The strong interaction is comprised in 
the SU(3) invariant lagrangian of Quantum Chromodynamics, >Cqcd [3]- Both combine to 
the Standard Model lagrangian Csu- 

The weak sector is the richest field of the Standard Model. First the electroweak gauge 
group SU(2)l X U(1)y is chiral, because left-handed fermion fields transform differently 
from the right-handed ones implying parity violating couplings. Second the SU(2)l x U(1)y 
gauge symmetry is spontaneously broken to the electromagnetic gauge group U(l)pjj^. This 
means that the physical vacuum state does not have the SU(2)l x U(l)Y-symmetry of the 
lagrangian. In the Standard Model the symmetry breaking is provided by the Higgs-Kibble 
mechanism [22] introducing a scalar potential in C with degenerate minima which trans- 
form into each other under SU(2)l x U(1)y transformations instead of being individually 
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invariant. As a consequence three gauge bosons become massive due to their interaction 
with the Higgs field. The chiral structure of the Standard Model forbids explicit mass terms 
for the fermions. Fcrmion masses are therefore generated through Yukawa interactions with 
the Higgs field. After diagonalizing the resulting quark mass matrix the coupling between 
the mass eigenstates of the quarks and the weak charged vector bosons is found to be 
of the form: 



where Uj denotes the up-type quarks of the three families, u,c,t, and dk stands for the 
down-type quarks d,s,b. The Yukawa interactions in £sm do not only induce fermion 
masses, but also allow for family-changing couplings in (2.14) described by the unitary 
Cabibbo-Kobayashi-Maskawa (CKM) matrix {Vjk) [6]. A unitary 3 x 3-matrix can be 
parametrized by nine parameters, five of which are redundant in (2.14), as they correspond 
to physically irrelevant redefinition of the quark fields by U(l) phase transformations. The 
other four correspond to a complex phase S and three rotation angles, i.e. V is real for S — 0. 
If V is real the charged current interaction Cy^c in (2.13) is CP-conserving. A non-zero 
5 corresponds to CP violation. This can easily be seen by recalling from sect. 2.2.3 that 
the latter is equivalent to T violation which is an antiunitary transformation involving 
a complex conjugation. Hence since jCy,c is T invariant for real V, it must in general 
transform under T into an expression with complex conjugate V. Since further Jcc 7^ Jcc 
this is different from the untransformed Cu,c- Yet only a non-zero S is not sufficient for 
CP violation, we must also ensure that it cannot be absorbed into a redefinition of the 
quark fields. This means that the three up-type quarks as well as the three down-type 
quarks must have different masses. In the lepton sector any mixing matrix can likewise be 
absorbed into a redefinition of the neutrino fields, if the neutrinos are massless. 
For the weak neutral current such fiavour-changing vertices do not occur as a consequence 
of the unitarity of V. This phenomenon is called GIM mechanism [12] implying the sup- 
pression of fiavour changing neutral current (FCNC) processes. 

The QCD lagrangian enjoys an unbroken SU(3) gauge symmetry. It is obtained by adding 
a fiavour label k to the fermion field in (2.10): 



with the colour indices now being suppressed. In (2.15) the fermion fields ipk{^) correspond 
to quarks with mass and the vector fields A'J^{x) correspond to gluons. 

2.4. Perturbation Theory 

When quantizing a theory we want to preserve the symmetries of the classical theory. I.e. 
we wish to have the vectors of the Hilbert space on which our field operators act transform 




(2.13) 



(2.14) 




(2.15) 
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according to a representation of the symmetry group. E.g. for some internal symmetry Lie 
group with generators and structure constants /"^'^ one seeks to find a set of operators 
satisfying 



and 



kj 



(2.16) 



(2.17) 



Here the second condition imphes the consistency between the symmetry transformation 
in the Hilbcrt space (LHS) with the corresponding transformation of the fields (RHS). 
In (2.17) an additional unphysical phase factor is allowed, too. Further microcausality 
requires that the operators corresponding to the Noether currents commute, if the distance 
of their arguments is spacehke. All these conditions are fulfilled, if one postulates canonical 
commutation relations for boson field operators and canonical anticommutation relations 
for fermion field operators: 



5C 



5 (do^k) 



= iS^^\x-y) for = 



(2.18) 



The operators are then obtained by simply replacing the classical fields (pj by (pj in the 
expressions for the classical Noether charges Q"'. Finally for an unbroken symmetry the 
vaccuum state |0) is chosen to be invariant with respect to the symmetry transformations. 



2.4.1. Green's Functions, S— Matrix 

In this section we consider for simplicity a theory with a single real scalar field ^(a;) . The 
vacuum-to-vacuum amplitude of a product of n field operators is called n-point Green's 
function: 



Gn{xi,...Xn) = (O|T0(a;i)-...0(a;O|O). 



(2.19) 



Here T means time ordering of the fields such that the time coordinates increase from 
right to left. For a noninteracting theory all Green's functions can be easily expressed in 
terms of the simplest nonvanishing Green's function 



G2{XI,X2) = 



1 



(2.20) 



with the Stiickelberg-Feynman propagator 

G{k) = 



^2 _ _|_ ^g- 



(2.21) 
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(2.18) implies that the Fourier components of a free field 0(x) satisfy the commutation 
relations of annihilation and creation operators. Hence one can use them to create n- 
particle states from the vacuum state |0). These n-particle states form the canonical basis 
of the Fock space JF. 

For an interacting field theory we neither know the form of G{k) in (2.20) nor can we 
relate the n-point function in (2.19) to the two-point function (2.20). On the other hand 
we empirically know that largely separated particles essentially behave as free particles. 
Further in a scattering experiment the relevant distances \xj—Xk\ and a;°| are large. The 
Riemann-Lebesgue lemma states that for \xj — x'^\ oo only non-integrable singularities 
of G{k) can give a contribution to the integral in (2.20). Hence one assumes that also for 
an interacting theory G{k) develops a pole as in (2.21): 

~ R 

G2(k) = — ^ h less singular terms . (2.22) 

This makes the concept of asymptotic states and fields plausible: One assumes that the 
Hilbert space of the interacting field theory contains some subset of state vectors which 
for t — > — oo develop into states which correspond to n-particle momentum eigenstates of 
the free hamiltonian. I.e. for t — > — oo any of these instates tends to an element of the 
canonical basis of a Fock space These states describe the incoming particles in scattering 
experiments. The in-statcs are created and annihilated by the Fourier components of a free 
field operator 0in. On the subspace spanned by the in-states one has 

${x) ~ R^'%n{x) (2.23) 

i— >— oo 

in the weak sense. Here one has to allow for a real constant factor i?, because (2.18) fixes 
an overall normalization of <^in. Next one postulates a field (j)ont{.x) and a set of out-states 
which have the described property for i — > oo corresponding to the outgoing particles in 
the particle collider. It is reasonable to assume that every particle species which can be 
used to prepare the incoming beam can possibly be observed in the final state and vice 
versa, thus the out-states span the same Fock space T as the in-states. Lorentz invariance 
implies the invariance of the vaccuum, |0) interacting = |0)in = |0)out; and of one-particle 
states, |p) interacting = |p)in = |p)out- Hcncc with (2.23) One finally obtains an interacting 
two-point function which behaves asymptotically for | a:;? — 2^2! ~^ °° accordance with 
(2.22). 

The scattering experiment is described by the transition amplitude or S-matrix element 

S^p = (a, out|/3, in) = (a, in|,S|/3, in). (2.24) 

With the concept of asymptotic fields and states it is now possible to relate the S-matrix 
elements to the Green's functions. This relationship is given by the reduction formula of 
Lehmann, Symanzik and Zimmermann (LSZ) [23]. It is particularly simple in terms of the 
Fourier transformed {n + Z)-point Green's function (27r)^5(X]"=i ^i) ' Gn+i{ki, . . . , kn+i-i)'- 

(pi, . . out|?i, ...qiin) = 
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n I 

. n+l 



R—{2i:)%Y.Pi -Y.^i)- GUi ■ ■ - -Pn, ?i, . ..qi-i) , (2.25) 



where we have for simphcity assumed that pi ^ qj for all In (2.25) the truncated 

Green's function G* is defined by 



Gl{ki,...,kr-i) = lim ... lim, Grih,...,K^^] 



(2.26) 



If vector bosons or fermions are scattered the fields in (2.25) carry indices which are con- 
tracted with spinors or polarization vectors of the external particles. 

From (2.25) and (2.26) one realizes that Green's functions carry a lot of redundant infor- 
mation, only their on-shell values are related to physical observables. 
Finally one is left with the evaluation of the Green's functions Gn{xi, . . . Xn). They can be 
most conveniently expressed in terms of the path integral: 

In perturbation theory one expresses G„ by a power series in the coupling constants of the 
lagrangian. In (2.27) this is done by expanding the kernel 



„ n. LJ 



(2.28) 



with >Cint being the interaction lagrangian. The individual terms of this perturbation series 
are usually depicted in a collection of Feynman diagrams. With (2.28) and (2.27) one 
realizes that one can formally write for the S-matrix in (2.24): 

g — T e^-^ d-^x^intl<l>{x),dij.'t>{x)] ^2 29) 

Within perturbation theory the assumption (2.22) is indeed verified: In any finite order of 
perturbation theory the interacting two-point function is of the form given in (2.22). mpoie 
and the residue R of the single particle pole are also power series in the couplings. 
In gauge theories the path integral (2.27) includes the integration over field configurations 
which are related to each other by gauge transformations. Further even the free propagator 
(2.20) of a gauge boson when derived from (2.10) does not exist. A remedy is the addition 

of a gauge fixing term Cg{ = — (^(9^/1^^ / (2^) with gauge parameter ^ and the introduction 
of Faddcev-Popov ghost fields rj"- via Cpp = {d^Tj"') Dl^rj^ [24]. The graphs with FP ghosts 
correct for the unphysical polarizations of the gauge fields corresponding to gauge degrees 
of freedom. FP ghosts are scalar Grassmann fields meaning that they anticommute. 
Hence the QCD lagrangian suitable for perturbative calculations reads: 

>CpQCD = V^.(^^-m,)V'.-^F;,F'^'^^-^(a'^^::)' + (aVXry'' (2.30) 
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One terms quantities which do not depend on the arbitrary gauge parameter ^ gauge inde- 
pendent. S-matrix elements and physical observables are gauge-independent. In contrast a 
field polynomial or Green's function is called gauge invariant, if it is invariant with respect 
to the gauge transformation (2.8) and (2.11). 

After all in the case of QCD the standard reduction procedure sketched above does not 
seem very useful, because in the very beginning the assumption (2.22) is in contradiction 
with the observed confinement property. The poles in the Green's functions are not related 
to the particles described by the elementary fields in the lagrangian as suggested by (2.25), 
but rather by bound states, hadrons, formed by them. The asymptotic states are hadron 
states which we cannot simply relate to quark and gluon fields. We may then ask how to 
make use of perturbative QCD at all. To answer this questions one has to mention a further 
property of perturbative quantum field theory discovered by Wilson [25]: It is possible to 
expand a Green's function in terms of some ratio of mass parameters and each term of 
this expansion is the product of a Green's function involving a composite operator and 
a coefficient being independent of the external states. This Wilson coefficient comprises 
the short distance physics whose scale is set by the inverse heavy mass parameter. Further 
QGD enjoys the feature of asymptotic freedom ensuring that the QGD coupling is small 
when probed at small distances. Hence for the calculation of the Wilson coefficients one 
can use perturbation theory and can further take external quark states rather than hadron 
states. 

Sometimes people try to describe hadron properties in terms of off-shell quark Green's 
functions. They, however, are unphysical, as mentioned in the discussion after (2.25), e.g. 
they are gauge dependent. This point will be relevant in the discussion of the use of the 
field equation of motion in sect. 3.10.1. 

2.4.2. BRS Invariance 

The lagrangian (2.30) is no more gauge invariant due to the gauge-fixing term and the 
Faddeev-Popov term. It seems as if we have lost all the nice features of local gauge sym- 
metry. This is not so, because (2.30) is still invariant under Becchi-Rouet-Stora (BRS) 
transformations [26]. They are constructed by factorizing the gauge parameter in (2.8) and 
(2.11) into the product of the FP field and a Grassmann parameter 5X: 



The infinitesimal transformations of the various fields read: 



-ig{r]''6\)T''i: = igT'^rj^'ipSX, 
igi^T^r^^SX, 



(5,^" + gP'^rfAl) 5X 

-Igr^'v'v'sx 



■a 



(2.31) 
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The basic QCD lagrangian (2.15) is clearly BRS invariant, because for the quark and gluon 
fields the BRS transformation is just a gauge transformation. The ghost fields transform 
such that £gf + £fp is invariant. The BRS symmetry is a global supersymmetry, because it 
involves an ^-independent anticommuting parameter d\. Its magic is that it encodes the 
original local gauge symmetry which had to be destroyed by adding £gf + Cpp to Cqcd- 
A global symmetry of C induces relations between Green's functions called Ward-Takahashi 
identities [28]. The Ward-Takahashi identities implied by BRS invariance are Slavnov- 
Taylor identities [29]. We will need them in the form 

5brs(0|TXF|0) = (0|T (5brsX)>^|0) + (0|TX(5brs>^)|0) = 0, (2.32) 

where X and Y are product of fields whose arguments may coincide. When making use of 
(2.32) one has to anticommute SX to the right (or to the left) thereby picking up a relative 
"-"-sign between the terms involving 5brs-^ and 5brs^) if they contain an odd number of 
Dirac fields. 

2.4.3. Renormalization 

The perturbative evaluation of Green's functions involves integrations over loop momenta. 
Some of these integrals diverge when the loop momenta get large, they are UV divergent. 
To give a meaning to these divergent expressions one first has to regularize the divergence: 
An additional parameter D is introduced in the lagrangian such that the latter coincides 
for D = Dq with the original lagrangian. The loop integrals are finite in some range 
for D and analytic functions in D. For D ^ Dq one rediscovers the divergence in the 
form of some singularity. Throughout this paper we use dimensional regularization which 
conserves most of the interesting symmetries. Here D is the space-time dimension and 
Do = 4. The appearance of UV divergences originates in the fact that quantum field 
operators are singular distributions, see (2.18). In general their product is ill-defined when 
their arguments coincide. The renormalization process gives a meaning to the field products 
in the interaction lagrangian. 

The renormalization process is carried out recursively in the power of the coupling constant 
g or equivalently in the number of loops. It is encoded in Zimmermann's forest formula [27]. 
Starting from the one-loop order one finds that one can absorb all divergences into a multi- 
plicative redefinition of the parameters of the lagrangian such as masses, couplings and the 
normalization of the fields which is connected to the quantity R in (2.23). These multiplica- 
tive renormalization constants are Laurent series in D — 4 in dimensional regularization. In 
the n-th order all divergences of 1- to (n-l)-loop subdiagrams are already remedied by the 
earlier renormalization steps. The remaining overall divergence can then be absorbed into 
new n-th order terms in the renormalization constants. For this to work it is crucial that 
after the subtraction of sub divergences the divergences are local, i.e. polynomial in the ex- 
ternal momenta of the Green's function under consideration. Further the renormalization 
constants must be universal meaning they must remove the divergence of some subloop in 
whatever diagram it occurs. 
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Yet for a predicative theory we must also ensure that the renormalization process does not 
induce new couplings in any new order of perturbation theory, because then we would end 
up with a theory with infinitely many undetermined parameters. This can be ensured, if 
the lagrangian contains only couplings with non-negative dimension. In this case we term 
the theory to be renormalizable by power counting. This criterion applies to the Standard 
Model. 

In the following we will also be confronted with effective field theories which contain cou- 
plings with negative dimension. The parameters of such effective field theories, however, 
will be fixed by a matching procedure in which the effective theory Green's functions are 
related to Standard Model Green's functions, so that the theory is predicative. 
Since we will always work to lowest nonvanishing order in the weak coupling constant, we 
will only have to discuss renormalization of the QCD lagrangian. First we define 

A-D 



The QCD Lagrangian now reads in terms of the unrenormalized (bare) parameters: 

+ (a^'Tj^^)'") (dpS"'' - g^""^ P^^Ap) Ty^o)'*- (2.33) 
The bare fields and parameters are related to the renormahzed ones by 

^(o),a ^ zTa;, ^(°) = zyv, r,'^ = 

- Z,g^\ = ZAi, m(°) = Z^m. (2.34) 

Here an arbitrary dimension-1 parameter, the renormalization scale /i has been introduced 
to keep the renormahzed coupling constant g dimcnsionless. In dimensional regularization 
the n-loop diagrams involve poles in e up to order n. Hence the renormalization constants 
are of the form: 




(2.35) 



Now the renormahzed lagrangian reads: 

>CpQCD = -^ZA{dpAt-d,A;){d''A''''-d''An-^idPA';y 

+Z^7p^i^7pk - Z^Zmmktpkipk + Zrfd^rfdpT]"- 

-Z.zT^IJfr"' {dpAl - a.A;) A^M- + ZgZ^zTgpfi^i^T'^^l^ 

2 

-ZlZ\^H^'P^'f^A'l,AlA''PA'^'' - ZgZ„Z][''gn'P^' {d^rf) A^rf (2.36) 
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Sometimes it is useful to decompose >CpQCD in (2.36) into a sum £o + >Cct, where Cq has 
the form of (2.33) with the bare quantities replaced by renormalized ones. C^t is called 

counterterm lagrangian. 

Gauge invariance has implied that all interaction terms in i2pQCD involve the same coupling 
constant g. Dimensional regularization has the advantage that it manifestly preserves gauge 
invariance in the unrenormalized Green's functions. As shown by 't Hooft the subsequent 
renormalization procedure does not spoil gauge invariance, so that Zg is the same in all 
these interaction terms [30]. Finally also renormalized Green's functions obey the Slavnov- 
Taylor identities, but with a renormalized parameter 



5K = z-^^'^z-"H\. 



(2.37) 
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3.1. Renormalization Schemes 

When determining the renormahzation constants in (2.36) one calculates Green's functions 
obtained from Cq + and adjusts the renormalization constants such that the result 
is finite. The Green's function to some fixed order in g involves n-loop diagrams with 
interaction vertices only from Cq and subloop counterterm diagrams with fewer loops and 
with vertices also from Cct- The individual contributions contain divergent terms which 
depend non-polynomial on the external momenta pi, but they cancel in the sum of every 
n-loop diagram and its sub-loop counterterm diagrams [27]. Since these non-local poles 
cancel, one can absorb the remaining divergences into n-loop counterterms. 
Now one is free to subtract any momentum-independent constant Zg"^ together with the 
divergences order by order in perturbation theory. Any such choice of the Zq"^'s in (2.35) 
corresponds to a different renormalization scheme. In the minimal subtraction scheme (MS) 
[31] one picks z'q'^ = 0. All schemes in which the zjf'^^s do not depend on the masses are 
called mass independent. Due to Weinberg's theorem [32] they enjoy the property that also 
the divergent parts of the Z^^^'s are independent of the masses (see e.g. [33]). In other words 
in mass independent schemes the divergences are not only polynomial in the momenta but 
also in the masses. The main advantage of these schemes, however, is that they allow for 
the solution of the renormalization group equations to be discussed in sect. 3.2. 
Since in every renormalization scheme the divergences are removed, any two schemes can 
at most differ by a finite renormalization, for example an arbitrary scheme is related to the 
MS scheme via 



1 ^ 7(1) 9 ^ 7(2) 9 

16^^^' (16^ 



(3.1) 



From this one sees that the Z^"^'s are scheme-independent while the other Laurent coef- 
ficients in (2.35) depend on the lower order zjf^^s. In fact the Z.,^"-''s are simply related 



to the one-loop Z\^'^s. The finite renormalization 1 + Z^^^g"^/ (IGvr^) + Zq^^ g"^ / {IGn'^^'^ 
in (3.1) may be viewed as a perturbative redefinition of the fields, masses and coupling 
constants in C in (2.36). 

In dimensional regularization all loop integrals involve an additive term 7e — ln(47r) with 7e 
being Euler's constant. It is reasonable to subtract this term in every order of perturbation 
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theory. The resulting scheme is called MS scheme [34] . There is a simple way to relate the 
MS scheme to the MS scheme: One first replaces the renormahzation scale // in (2.34) by 



Recalling that /i rcsp. JI has been introduced into (2.36) to maintain integer dimensions 
for the renormalized parameters one realizes that each loop integral / d^q comes with a 
factor of Jl'^^. Then one expands (3.2) to the desired order in e and multiplies this power 
series with the Laurent series in e of the unrenormalized diagrams. In this way one modifies 
the pole part of the unrenormalized Green's functions such that 7e and ln(47r) vanishes. 
Finally one simply applies minimal subtraction to the remaining expression. This method, 
which will be used in this thesis, has the advantage that one can avoid Euler's constant and 
ln(47r) in intermediate steps of the calculation. Further one can easily translate relations 
proven in the MS scheme into the MS scheme. 

Yet there is another aspect of (3.2): When inserted into (2.36) one may view (3.2) as a 
redefinition of (7 by a (7-independent power scries in e, 1 + Y^k=i (^k^^, or alternatively as a 
redefinition of the field monomials multiplied by g. Here in (3.2) this series reads 



Clearly multiplying with such a power series does not commute with the renormalization 
process, so that it corresponds to a change of the renormalization scheme. Let us call this 
procedure evanescent redefinition of the coupling or of the field monomial. Such evanescent 
redefinitions appear naturally in the context of effective four-fermion couplings. We will 
investigate them in chapter 4 [19]. One task will be to translate an evanescent redefinition 
into a finite renormalization (3.1). We finally remark that the change from the MS to the 
MS scheme is equivalent to a change of the scale /i as evidenced from (3.2). 
Any S-matrix element and any physical observable is independent of the chosen renormal- 
ization scheme to the calculated order. This means, if one calculates some observable in 
scheme 1 to order (yfj" and afterwards expresses gj and the other parameters such as masses 
in terms of the coupling gu and the parameters of some other scheme II, the result differs 
at most by terms of order g^^^ from the calculation of the observable in scheme II. 

3.2. Renormalization Group, Aqcd 

The invariance of the S-matrix with respect to the change of the renormalization prescrip- 
tion implies a continous symmetry of the quantized theory as noticed first by Stiickelberg 
and Pctcrmann [35]. The group associated to this symmetry is the renormalization group 
(RG). Its power roots in the fact that it hnks different orders of the perturbative series: 
By applying a renormalization group transformation to some n-th order Green's function 
we can add some logarithmic term of all uncalculated higher orders to the result. When 




(3.2) 




(3.3) 
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the RG is applied properly, these logarithms are the dominant contribution of the higher 
orders and they are summed to all orders by the RG transformation. 
In mass independent schemes renormalized Green's functions depend on the renormaliza- 
tion scale fi. Changing fi corresponds to a very simple way of modifying the renormalization 
prescription as discussed after (3.2). The S-matrix is independent of /i, if the bare quan- 
tities in (2.33) are chosen //-independent. Hence the renormalization scale invariance is a 
special case of the renormalization scheme invariance. 



3.2.1. Renormalization Group Functions in QCD 

In this section we will collect the important RG functions of QCD in the MS scheme (see 
e.g. [36]). From the //-independence of the unrenormalized parameters in (2.34) one easily 
finds: 

^ f^-^ = -T^lm{g{lJ)) (3.4) 



with the beta-function 



mi^)) = -'9-Y^^9 (3.5) 



and the anomalous dimensions of quark mass and quark field 



^ dZm / / dZ^^'^ 



m in (3.4) is the current mass appearing in the lagrangian (2.36). It is also called running 
mass, because it depends on /i. The renormalization group functions in mass independent 
schemes depend on /( only through g{ij,) as indicated in (3.4). This feature allows to solve 
the renormalization group equation, which reads for a renormalized truncated connected 
n-quark Green's function: 



d 

^ djj, 



z:^,'''\g{^^))&:{p,,g{^l)MM).l4^ = 0. (3.7) 



Here /f^ means 



where the partial derivatives are performed with the bare quantities kept fixed. 
Let us briefly sketch the relation of (3.7) to the scale transformations alluded to in sect. 
2.2.1: Recalling that /i has been introduced to maintain an integer dimension d for &^ one 
easily flnds (see e.g [36, p. 199]): 

d \ ~ 

g^+X— + d\G'^{Xpi,g{li),m{g{ii)),ii) = 0, 
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which may be used to ehminate /x^ in (3.7) in favour of the scahng factor A of the external 
momenta. 

Now n enters Green's functions only in powers of e, thus after expanding in e only in 
powers of logarithms. Since (3{g) = 0{g^), 7^ = 0(5-^), 7^ = Oig"^) and /i^ = (3.7) 
connects the In /x's of different orders in g of G*^ to each other. 

We will need the renormalization group functions up to the next-to leading order (NLO) 
[36]. The beta-function reads: 

with 

/?W = = |jv2-^iV/-2C;./. (3.10) 

Here is the number of colours, Cp — {N^ — 1)/{2N) and / is the number of quark 
flavours. For the anomalous mass dimension one finds: 

with 

7i°^ = 6C^, 7^^^^) = ^^^(sC^ + fA^-y/). (3-11) 

Whenever we discuss leading order (LO) expressions, only the first term is kept in (3.9) 
and (3.11). The anomalous dimension of the quark field will only be needed in the leading 
order: 

7. - ^'^'^2+0{9') (3.12) 

with 

7^'^ = Cf^- (3.13) 

With these definitions it is easy to write down the solutions of the RG equations (3.4) for 
«(/.) = ^^(/.)/(47r): 

^ ^ 1 /3F^log[log(/iVA^)] ^ ^/ log^[...] \ .3,4. 

The mass parameter A/ has entered (3.14) as a constant of integration. Hence the massless 
coupling g of the classical lagrangian has become a function of the dimension-one parameter 
Ay in the quantized theory. (3.14) also exhibits the property of asymptotic freedom, «(//) — > 
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for — > oo. By inverting (3.14) one can define A/. For a defined in the MS scheme one 
has 



-■>!" 



'^r-M +0(a(/i)) 



(3.15) 



In the hmit // — > oo this is a precise definition of Aj^. Since the NLO coefficient /?i appears 
in an overall factor rather than in a radiative correction, one realizes that it is mandatory 
to go to the NLO to define Aj^. Conversely only in NLO expressions for some physical 
observable of interest one can make use of Aj^. This is an important reason to perform 
perturbative calculations beyond the leading order. Af has been defined in the MS scheme 
by Buras, Floratos, Ross and Sachrajda [37] and the MS scheme has been introduced by 
Bardeen, Buras, Duke and Muta [34]. 

Solving the RG equation (3.4) for the running quark mass yields for the mass at some scale 
II expressed in terms of m(/x = m): 

In (3.16) S£ and J^^ are defined by 

d^rP = ^ and J({) = +f ^ 

If only contains a single mass parameter apart from /x, for example some spacelike 
external momentum \/—q'^, Infi^ must necessarily appear in the form ln(— g^/^^) To make 
sense of a perturbative calculation wc must choose fi^ ~ —q^ in order to avoid a large 
logarithm multiplying the expansion parameter q;(/x). One may solve (3.7) to obtain G^^ at 
any other scale /x. This RG improved expression differs from a standard perturbative result 
such that it inchides the terms involving Q;(/i) ln(— g^//i^) to all orders in perturbation 
theory. The LO RG improved perturbation theory is therefore called leading log (LL) 
approximation. 



3.3. A First Look at K^-RO -Mixing 

Yet in most physical processes more than one mass scale is involved. To make use of 
RG improved perturbation theory more field theoretic tools are needed. To motivate the 
following sections let us have a first look at the main subject of this thesis, the |AS| = 2 
transition inducing K'^ — K'^ -mixing. 

Fig. 5.1 on p. 61 shows the lowest order AS = —2 transition amplitude in the Standard 
Model. One-gluon radiative corrections are depicted in fig. 5.2 on p. 64. Now this amplitude 
is not useful to describe physical | AS| —2 processes for the following reasons: 
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i) The diagrams of fig. 5.2 incorporate QCD corrections perturbatively. While this is 
appropriate for short distance effects, the long distance QCD interaction is non- 
perturbative and cannot be described by the exchange of gluons. 

ii) There are largely separated mass scales involved: The radiative corrections of fig. 
5.2 contain large logarithms such as a\n{mc/Mw) = 0(1), so that the radiative 
corrections have the same order of magnitude as the leading term. One would like to 
sum this logarithm to all orders in perturbation theory. RG techniques applied to the 
Standard Model amplitude, however, will not achieve this, because with them only 
logarithms of /j, can be summed to all orders. 

iii) For the same reason it is not clear at which scale the running coupling should be 
evaluated. Between /i — rric and /i — Mw the coupling a{ii) changes roughly by a 
factor of three. 

iv) The real external states are mesons rather than on-shell quarks. Off-shell quarks are 
an even worse description of mesons, because off-shell amplitudes are unphysical as 
mentioned in the end of sect. 2.4.1. 

Hence to tackle the problem one wishes to separate short distance effects from the long 
distance physics. This separation is provided by Wilson's operator product expansion [25] 
discussed in the following section. 



3.4. Operator Product Expansion 

In weak processes involving hadrons one is confronted with two couplings: While ordinary 
perturbation theory in the weak coupling g^j works well, we have already seen in sect. 3.3 
that we need a framework to include all-order effects in the strong coupling g and possibly 

non-perturbative features of QCD. 

For this we expand the kernel of the path integral cxp [i J Csu] in terms of the clcctrowcak 
(and Higgs) interaction lagrangian, but keep the QCD interaction in the exponential. I.e. 
we expand as in (2.28), but with the QCD action added to Si^ee instead of being kept in 
Cint- Hence with (2.13) a four-quark Green's function corresponding to a weak charged 
current tree level process is found as 

2 t 9^ — — 



-j^d^'y J d''^Jlc'iy)w^iy)Jcciz)w-iz) + h.c] |o)*^(3.i7) 



Here the quark fields arc free fields with respect to the electroweak interaction, but in- 
teracting fields with respect to QCD. (3.17) represents all diagrams with one W-boson 
exchange and an arbitrary number of gluons. The index 2 in Cf^^'^ means second order 
in the weak coupling constant. The weak current J^^ is the first example of a composite 
operator, as it contains interacting fields at the same space-time point. 



3.4 Operator Product Expansion 
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Now Wilson has shown that one can expand a product of field operators in a series of 
composite operators with increasing dimension. For the product of the two currents in 

(3.17) this reads^: 

Ji'WM = ^r(^ - ^)^. (4^) ' (3-18) 

with the Qj's being local four-quark operators as depicted in fig. 5.7 on p. 69. The Green's 
functions on the RHS of the operator product expansion (OPE) (3.18) are the matrix ele- 
ments of Qj. The Wilson coefficients Cj" are functions which are singular ai y = z. The 
coefficient multiplying the operator with the lowest dimension is the most singular one 
by power counting. They comprise the short distance behaviour of the current product in 

(3.18) when y approaches z. After inserting (3.18) into (3.17) one performs the path inte- 
gration over the W-boson fields. By this the open Lorentz indices in (3.18) are contracted 
with a W propagator iVj^^y — z, Mw, A*)- With 

C,(M^,/.) = -^/rf^C^r(C)^^r(C^H^>/^) (3-19) 

one reahzes that (3.18) turns into an expansion in the (inverse) W mass Mw- With (3.19) 
one finally obtains for (3.17): 

Q(MH.,/^)(0|TV^l(xl)V^2(a^2)V'3(a;3)V'4(a;4)^/A4•(?/)|0)*'■ (3-20) 



Here we have absorbed the weak coupling constant into Cj. When we pass to the operator 
basis for K° — K° -mixing in chapter 5 we will extract a factor of 

V2 " 8M^ ^"^-"'^ 

out of Cj to show that the operator product expansion in weak decays replaces the Standard 
Model weak interaction by a four-fermion interaction with the Fermi constant Gp as the 
coupling constant. We will also sometimes leave out the integral over the space-time coor- 
dinate y of the composite operator in (3.20). This merely leaves out an overall momentum 
conserving 5-function when (3.20) is Fourier transformed into momentum space. 
In momentum space (3.20) corresponds to an expansion in m^gh^/M^, where mught stands 
collectively for the light^ quark masses and the small external momenta and Mandelstam 
variables. The local operators in (3.18) may be obtained by expanding the momentum 
space W propagator^: 



— 1 —1 —p 



2 



^sum on repeated indices 
^lighter than 

^We use the Feynman rules of [36]. 
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While it is trivial from (3.22) that (3.20) works for tree level diagrams with Cj ' — Qw/^-, 
it is not obvious that (3.20) can be extended to diagrams with gluon loops, in which p in 
(3.22) is a loop momentum. From Wilson's work we know that (3.18) holds to any order 
in perturbation theory. Clearly for small loop momenta (3.22) is a good expansion, and 
the contributions from small loop momenta on each side of (3.20) will match with the tree 
level coefficient Cf'^ The contributions from large loop momenta, however, will modify Cj 

by adding terms a/ (47r) • C]^\ . . to Cj°\ Yet these UV momenta are not sensitive to the 
small internal masses and small external momenta. 

When the heavy top quark is involved in the amplitude of interest, it must also be integrated 
out like the W-boson. 

In the discussion above we have made plausible the following properties of the operator 
product expansion: 

i) The Cj 's do not depend on the infrared structure of the Standard Model amphtude. 

ii) The Wilson coefficients Cj depend only on the heavy mass scales Mw and mtop and 
on n, while the matrix elements contain the light masses. Large logarithms such as 
\n{mc/Mw) are split into ln{fi/Mw) + ln(mc//i). The former logarithm resides in the 
Wilson coefficient and the latter is contained in the matrix elements. 

iii) RG methods are applicable to Cj allowing for the summation of the logarithm in 
all orders of perturbation theory. With the boundary conditions of the RG evolution 
being /i = Mw and /i — rUc one sums the large logarithm \n.{mc/Mw) of the Standard 
Model amphtude. 

iv) The Cj's are independent of the external states. 

Hence we have overcome the problems i) to iv) mentioned in sect. 3.3. The short distance 
coefficients Cj are physically sensible quantities, which can be obtained in perturbation 
theory. Clearly in the end one has to cope with hadronic matrix elements, which must be 
evaluated by nonpcrturbativc methods such as lattice gauge theory, expansion or sum 
rule techniques. Yet one can sometimes find them from experimental data and insert the 
obtained values into predictions for other observables of interest. The dimensions of the 
Qj's in (3.18) now manifest themself in powers of the hadronic scale Aqcd- 
The price paid for the factorization of short and long distance effects in (3.18) is the 
inclusion of only finite powers of m^gh^/M^. These powers are modified by logarithms 
which come with powers of a and can be summed to all orders in a by use of the RG. 
Since the light quark masses are much smaller than Mw in most cases only the operator 
with the smallest dimension in (3.18) is sufficient. When ji passes the flavour thresholds 
777.6 or fric-i however, we repeat the factorization process. Below 777c the OPE is an expansion 
in Aqcd/"^c and higher dimension operators may become important. 



3.5 Effective Field Theories 
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3.4.1. Renormalization of Composite Operators 

Since composite operators involve products of interacting fields at the same space-time 
point, they require a renormalization in addition to that of the fields they are composed 
of. This operator renormalization has been worked out by Zimmermann [38], who called 
renormalized operators normal products. 

The unrenormalized field product composed of unrenormalized fields will be denoted by 
the superscript bare. We will frequently be concerned with four-fermion operators: 

gbare ^ ^'^^'"''g^^bare . ^bare .^^bare^ UO SUm OU A;, (3.23) 

where qk and qk are strings of Dirac matrices. We will only need four-fermion operators 
which are Lorentz singlets or pscudosinglets, i.e. all Lorentz indices in Qk = Qk ® qk are 
contracted. Colour and fiavour indices are suppressed in (3.23), we will make them explicit 
where necessary. Next in the matrix elements the fields creating the external states are 
conventionally not displayed. After expressing the bare fields in (3.23) in terms of renor- 
mahzed fields as in ^2.34), the matrix elements are still divergent. The necessary operator 
renormalization of can require counterterms proportional to other operators Q^^"^^. 

This phenomenon is called operator mixing. We say that Qj mixes into Qk, if QY^"^ re- 
quires a counterterm proportional to Qk ■ Hence the multiplicative renormalization involves 
matrices renormalizing vectors of operators Q. The renormalized operator reads 

QA^) = QTi^) = zjkQt'^i^)- (3.24) 

As in (2.36) a quantity without superscript is understood to be renormalized. In cases 
where this may be confusing (such as in chapter 4) we will mark renormalized quantities 
with ren. In the phenomenological sections also the caret on the operators will be omitted. 
A set of operators which renormalize each other in (3.24) is said to close under renormal- 
ization. 

In (3.24) the standard definition of Q^^ has been given. It would be more useful to include 
the field renormahzation into the defintion of Qf^., so that would be expressed in 
terms of renormalized fields instead of bare fields. Again, in ambiguous places completely 
renormalized (i.e. finite) Green's functions will get a superscript ren: 

(Qj) = {Qfy^^ = = (3.25) 

We will also have to discuss the mixing of four-fermion operators with operators containing 
only two external quark lines, but additional gluon or ghost legs. They involve the product 
of the corresponding field renormalization constants in (3.25) instead of four quark field 
constants Z^^"^. 

3.5. Effective Field Theories 

The operator product expansion in terms of some inverse heavy mass described in the 
previous section has a complementary interpretation in terms of an effective field theory: 
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Since the Wilson coefficients Cj on the RHS of (3.20) arc independent of the external states 
we may consider them as coupling constants multiplying four-fermion couplings in some 
effective Lagrangian: 

= CjZj^Zli)^{x)qkip2{x)il)^{x)qkil}i{x) (3.26) 



Clearly the S-matrix 



T exp 



^Id^x (/:qcd(x)+£^(x 



(3.27) 



derived from (3.26) yields to first order in the effective couphng constants Cj and to all 
orders in the QCD coupling constant g the RHS of (3.20). We may compare the effective 
Lagrangian (3.26) with the QCD Lagrangian (2.36): The last line in (3.26) corresponds to 
a renormalized interaction vertex expressed in terms of renormalized fields. We may view 
Zjj^ as renormalizing the effective coupling constant Cj instead of the composite operator. 
Yet the concept of an effective field theory must allow for going beyond the first order in 
the effective interaction in (3.27). This is possible as we know from the work of Witten [39]. 
The Green's functions of second order in the effective interaction involve 



{-\jd''xjd''yTC\x)C\y)). 



These bilocal structures are in general divergent even in the absence of QCD interactions 
as can be seen from fig. 5.4. Hence one must add counterterms proportional to local oper- 
ators to Cqcd + d in order to obtain a finite result. Although the effective theory is not 
renormalizable by power counting, this is possible in all orders of both the effective and 
the QCD interaction, because after subtracting subdivergences the remaining divergences 
are polynomial in the momenta. The non-renormalizable interaction, however, forces us to 
introduce new counterterms in every order of perturbation theory. Hence we cannot include 
all-order effects in the effective action as we can do with the RG summation of the QCD 
interaction. Yet ordinary perturbation theory works well for the weak interaction, so that 
working to finite order in the effective coupling is perfectly sufficient. 
In this thesis we will work up to the second order in the effective couplings. After adding 
the local operator counterterms to one ends up with: 

= £i + C,(/i)C,,(/x)Z,-,i//i)Q^^ + a(/x)Z,f (/.)Q|'-^ 
= Cj Zj^Zl ^{x)qki){x) i){x)qkip{x) 

+CkCk' Z,^,, Zl ii{x)q'r^{x)ii{x)q'r^{x) 

+Ci Zl^^Zl ^{x)q'r^P{x) i^{x)itr^{x) (3.28) 

with Qi — ijj{x)qj°'^ijj{x)ijj{x)qf"^ijj{x) , no sum on I, being the local operators needed as 
counterterms to the bilocal structures. In K° — K^^ -mixing the Q^s are |AS| = 1 operators 



3.6 Appelquist-Carrazone Theorem 
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and the Qi^s are |AS|=2 operators. In (3.26) and (3.28) we have for simphcity only dis- 
played four-fermion operator counterterms. If none of the quark fields in Q^, Qi has the 
same flavour quantum number as one of the antiquark fields, only four-quark countert- 
erms appear in C^^. Yet in the general case so called penguin diagrams lead to a mixing of 
four-quark operators into other operators as depicted in figs. 5.8 (p. 69) and 5.11 (p. 70). 
Prom (3.28) we reahze the advantage of the effective lagrangian picture: The extra terms 
in (3.28) exactly look hke renormalization factors for the effective couplings Cj, but now 
both in terms of the effective interaction and of QCD: 

/-ibare 'Z~l/^ i i 

'^l — '^rl -r Z.^j^, ;OfeOfe' -h 

This form resembles the QCD coupling renormalization in (2.36). Further the effective la- 
grangian form allows for a simple inclusion of the correct wave function renormalization, 
which is simply taken into account by the corresponding counterterm Feynman rules de- 
rived from >Cqcd + Finally the renormalization group equations for the C^'s and C;'s 
to be discussed in sect. 3.8 are very easily found from (3.26) and (3.28). 

3.6. Appelquist-Carrazone Theorem 

The decoupling theorem of Appelquist and Carrazone [40] deals with the effect of heavy 
virtual particles on low energy processes: If we can remove some heavy particle field from 
the lagrangian without spoiling the renormalizability^, the loop effects of the heavy particle 
can be absorbed in the renormalization constants of the low energy lagrangian (without 
the heavy particle field) and into additive terms which are suppressed by powers (times 
logarithms) of the heavy particle mass. 

The Appelquist-Carrazone theorem puts some hierarchical order into field theory: It has 
allowed to discover QED without knowing about its embedding into the Standard Model 
and it allows to study Standard Model predictions at LEP energies without knowing the 

physics at, say, 10^^ GcV. 

We will be concerned with the Appelquist-Carrazone theorem in a twofold way: First we 
will use that quark fields and clearly also the W field do decouple with respect to the strong 
interaction: QCD works with 5 fiavours as well as with 6, hence we can remove the top 
quark from the lagrangian when passing with the renormalization scale n below its mass. 
Wc repeat this with the W mass and then when crossing the b- and c-threshold. We could 
even incorporate these inverse power corrections by taking into account operators with 
higher dimension in (3.28) and we can use the RG to sum the logarithms accompanying 
these powers together with powers of the strong coupling constant. 

The weak interaction, however, is a prominent example for non-decoupling: The heavy 
top-quark is clearly needed as an isospin partner of the bottom quark and after removing 
the W-boson we are faced with non-renormalizable effective four-fermion couplings. The 

■^We will use the term renormalizability for two different properties: It either refers to the locality of 
counterterms, in this sense we also call an effective theory renormalizable, or it means renormalizability by 
power counting such that all couplings have non-negative dimension. 
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non-decoupling allows to study the weak interaction in hadronic processes, although their 
energy scale Aqcd is much lower than the scale M^, at which the weak transition takes 
place. 



3.7. Initial Condition for the Wilson Coefficients 

RG improved perturbation theory starts with the definition of the initial condition for the 
Wilson coefficients. This is done by matching the Standard Model amplitude to the effective 
theory as prescribed in (3.20). For definiteness we will consider some weak four-fermion 
amplitudes. 



3.7.1. Single Operator Insertion 

Let's look at some weak tree level process such as a | AS| = 1 transition, which is mediated 
by the exchange of one W boson. The corresponding Standard Model amplitude factorizes 
order by order in 1/M^ into a Wilson coefficient vector C — (Ci, . . .CnY and matrix 

— * 

elements {Q). For illustration we only display the dependence on the renormalization scale 
the W mass and a light quark mass m: 

-iGf'\Mw,^i,m) = Cj{Mw,l^){Qj){l^,m). (3.29) 

A |AS| = 1 four-quark operator is depicted in fig. 5.7 on p. 69. 

To obtain (3.29) one has to calculate the Standard Model amplitude and the matrix element 
in perturbation theory to the desired order in the QCD coupling g. By comparing both 
results one finds Cj(Mvy,/x). For a meaningful perturbative result for Cj the factorization 
scale /i = /io must be chosen to be of the order M^i/ to keep the logarithm a\n{fl/M^r) 
in Cj small. Thereby one obtains the initial value for Cj{Mw, f^o) with respect to the RG 
evolution. With the RG evolution of Cj from /j, = /iq down to = Hm ~ m described 
in the following section we want to sum Indian/ IJ'o) ~ ln{m/Mw) present in the LHS of 
(3.29) to all orders. As indicated one is not forced to choose fj,m = m, /iq = Mw exactly. 
The two large logarithms ln(//,m//xo) and \n{m/Mw) differ by a small logarithm of the 
order of magnitude of the remaining non-summed radiative correction. Changing /iq and 
lira modifies the result by terms of the order of the uncalculated higher order corrections. 
These scale dependences are inherent to RG improved perturbation theory and may be 
viewed as an estimate of the accuracy of the truncated perturbation series. For the rest of 
this section we choose fiQ = Myy and Hm = ^ for clarity. 

From the matching procedure one can count the powers of the summed logarithms: In LO 
the matching is done in the order Ivl' {m / Mw) and the LO RG will sum 

777 

a^l^n 71 = 0,1,2,.... 

Mw 

In NLO we will sum the log's with one additional power of a. Hence the matching is done 
in the order h^{m/Mw), which requires the calculation of the one-loop diagrams of fig. 
4.1 on p. 43. 



3.7 Initial Condition for the Wilson Coefficients 
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3.7.2. Double Operator Insertion 

Next we will look at a process which appears in the second order of the effective couplings. 
Clearly we have in mind the |AS|=2 transition with light internal quarks. Its Standard 
Model amphtude is described in LO by fig. 5.1 (p. 61) and the effective theory matrix 
element is shown in fig. 5.4 (p. 66). Wc will often refer to fig. 5.4 as the prototype for a 
matrix element with two insertion of four-fcrmion operators, although we will later also 
discuss matrix elements containing one other operator such the one depicted in fig. 5.9 
(p. 69) or 5.11 (p. 70). 
Here the matching condition reads 



T QkQk' 



ren\ ren 



+ Ci (Qi) . (3.30) 



Here the square brackets around the two operators shall denote the inclusion of the local 

operator countertcrms in C^^ in (3.28). 

Next let us comment on the powers of the logarithm. There are two different cases, both 
of which appear in K^ — K*^ -mixing: 

Case I: 

This generic case is the following situation: 

The diagram in fig. 5.4 is divergent and requires local operator counterterms depicted in 
fig. 5.3 (p. 64). In the example of K° — K° -mixing these local operators Qi are |AS|=2 
operators. With the divergence comes ln(/x/m) in the finite part even in the absence of 
QCD corrections. Its twin in the Standard Model amplitude is \n(m/Mw). In the LO we 
can simply do the matching from these logarithms and the leading log summation comprises 

a"ln"+'-^, n = 0,l,2,.... 
Mw 

Consequently in NLO we can still do the matching from the one-loop graphs figs. 5.1 
and 5.4, but now from the a^ln'^{m/Mw) parts. Here we must also take into account the 
coefficients of the Qi^s. 

Case 1 will appear in the calculation of 773 in (1.1) in chapter 6. Cf. also the NLO calculation 
of rare K decays in [41]. 

Case II: 

The insertion of two local operators into the diagram of fig. 5.4 is convergent. This situation 
appears when there are several contributions of the form in fig. 5.4 combining such that the 
divergences of the individual contributions cancel. With the cancellation of the divergences 
also the In //'s will disappear, so that there is no large log in the corresponding one-loop 
amplitude of fig. 5.1. Hence the leading log summation will comprise 

777 

aMn'^-— , 77 = 0,1,2,... 

Mw 
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and the matching has to be done in the order a°ln {m/Mw)- Consequently we are forced 
to perform the NLO matching in the order \vP{m/Mw)- This requires the calculation of 
the finite parts of the diagrams of fig. 5.2 (p. 64) and of those of fig. 5.5 (p. 66). 
Case II appears in the calculation of rji. The cancellation of the large logs in the corre- 
sponding Standard Model amplitude is caused by the GIM mechanism: The divergence of 
the diagram in fig. 5.4 is multiplied with the product of two columns of the CKM matrix, 
which is zero due to the unitarity of the latter. The finite part is of the order m^/M^ and 
one speaks of a hard GIM suppression. One could even term this suppression superhard, 
because there is even no ln(m^/M^) multiplying m'^/My^r and the RG enhancement is 
thereby suppressed by one power of a compared to case I. 



3.8. Renormalization Group for Operators and Coefficients 
3.8.1. Single Operator Insertion 

In this section we will set up the RG formalism for matrix elements with single insertions 
of composite operators, which can be derived from {!} in (3.26). 



Anomalous Dimension Matrix 

The anomalous dimension matrix 7 is defined by 



lijidifJ')) 



7-1 d 7 



(3.31) 



From (3.24) one simply obtains the RG equation for the renormalized operator: 



0. 



(3.32) 



We expand Zjh as usual in g and e as shown in (2.35). The anomalous dimension matrix 



is related to Zj^ as follows: 



.(0) 



(0) 9'^ (1) 5'^ 



(167r2)^ 



+ o(/) 



(3.33) 



= -2Zi^^ - 2eZl 
1^'^ = -4Zf)+2{zW,zW}-2/3oZ«+0(5). 



(3.34) 



In (3.33) and in (3.34) we have only kept the terms relevant for NLO calculations. We have 
allowed for a finite renormalization Z^^^ as well. The structure of (3.34) reveals that one 
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obtains the correct 7*^^^ by inserting the finite one-loop counterterms with a factor of 1/2 
instead of 1 into the one-loop counterterm diagrams. (3.34) can be easily derived from the 
insertion of (3.33) into (3.31), see e.g. [43]. Alternatively one can recall that is related 
to the coefficient of In /i in {Q) . In dimensional regularization the In /i's enter the finite 
parts of Green's functions via 

{^Ji^'-l)- = 21n/., 

e 

where the —1 stems from counterterm diagrams. Hence e.g. in lowest order the coefficient 
of In /X is simply twice the coefficient of the divergence. 
In the same way one finds for (5q in (3.9): 

/3o = -2Z^gl. (3.35) 

As mentioned before, the different coefficients of Z are not independent. In NLO one has 
the relation 

4Zf^ + 2/3o4^^ -2Zf^zf^ = 0. (3.36) 

(3.36) is a consequence of the renormalizability, i.e. of the locality of the counterterms 
[42]. Since in mass independent schemes the divergence of each diagram plus its subloop 
counterterm diagrams is polynomial in the momenta and masses, any dependence on In jjL 
also cancels from the divergent parts. The two-loop diagrams involve ^'^^Z^^ je^ ^ while the 
one-loop counterterm diagrams involve ^^Z^^ je times either Zi'^ /e or z'^gj^/e = — /9o/(2e). 
The vanishing of In/i enforces (3.36). We will use this argument in a different context in 
a formal proof on evanescent operators, which vanish in D = A dimensions, in chapter 4. 
Relation (3.36) is tightly related to the finiteness of anomalous dimensions, which has been 
used for an elegant proof of (3.36) in [43]. 



Solution of the RG Equations 

Let us first derive the R( 
one immediately obtains: 



Let us first derive the RG equation for the Wilson coefficients: Prom n4-C^ = in (3.26) 



d 



a 



0, 



(3.37) 



thus the RG equation for the Wilson coefficient vector C involves the transpose of 7. The 
solution of (3.37) is given by 



Cjil^i) = [U{i^i,i^o)]jkCk{i^o) 
where C/(//i,//o) is the evolution matrix: 

" r9(Ml) ,7'^(C/') 



Tg exp 



/ dg 



P{g') 



(3.38) 



(3.39) 
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Here T g means that the matrices "y{g'),'y{g") ... in the expanded exponential are ordered 

increase from right to left for g{ixo) < g{l^i) \ 
decrease from right to left for g'(//o) > fl'(A*i) J 
The LO solution of (3.39) is 



such that the couplings 



C/(°^(//i,//o) = exp 



a! In 



a(/fo) 



(3.40) 



with 



7 



(0) 



2/3o 



(3.41) 



In the diagonal basis of d the LO evolution matrix U^'^^ is simply a diagonal matrix with 
diagonal elements [a{fio) / a{fii)]'^^ , where the dj's are the eigenvalues of d. 
The NLO solution of (3.39) reads: 



t/(/^i,/^o) = (l + ^j)t/(°n/^i,/^o) 



a(/fo) 
An 



J 



(3.42) 



Here clearly the NLO running a has to be used in U^^\ J in (3.42) is a solution of the 
matrix equation [44]: 



J + 



7 



(0)T 



7 



(i)r 



+ 



2Po 



2 ' 



(3.43) 



In practical calculations it is not useful to solve (3.43) by transforming into the diagonal 
basis. On better solves (3.43) directly: Since it is a system of hnear equations for the 
elements of J, the entries of J are simply rational functions of N and /. 
For the LO evolution matrix it is useful to transform to the diagonal form, as this saves 
computation time, because U^^^ has to be calculated newly for every new pair (/io,/ii). In 
our phenomenological analysis of chapter 7, however, we have simply exponentiated the 
matrix in (3.40), because the computer algebra system Mathematica provides a fairly 
fast numerical matrix exponentiation algorithm. 



3.8.2. Double Operator Insertion 

The local operator counterterms Z^^ ^{ji) Qi in C^^ do not influence the RG evolution of the 
coefficients C^, but they modify the running of C/. We will investigate this in the following. 
Often one does not insert the full set of operators in JC} into both places in fig. 5.4, but 
limits oneself to a subset which closes under renormalization. Hence k and k' in ; may 
run over different ranges. This will be so in the case of 773. 
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Anomalous Dimension Tensor and RG Equations 



The evolution of the Wilson coefficients in is fixed with (3.38). Now C^^ contains new 
local operator counterterms with coefficients Q. It has become standard to define the 
operators Qi with two inverse powers of g, e.g. Q — rn} j {g'^Jt'^^) • 7/x (1 — 75) (8) 7'' (1 — 75), 
to avoid operator mixing in the order g'^. Then Z^^i — 0{g^) and the mixing matrices start 
in the order g^. 

From fi-^C^ = and = one finds with (3.28): 



This can be compactly rewritten as 



d 



/^^c'K/^) = Ci>{ii)^i>i + Ck{ij)Ck'{ij)^kk',i 



with the anomalous dimension tensor 



^kn,l 



J— 



' 9^ y (1) 



— — [ikk'^nn' + ^k'lnn'] ^k'n',l'^l' 



VI- 



The analogue of (3.34) is found as 



ikn,l 
/ kn.l 



-1,(1) 



kni 



^1,(2) 



kn,l 
(1) 



+ 2e 
+ 2/^0 



-1,(1) 



knI 



-1,(1) 



kn,l' 



1,(1) 



- kn,l 

- 2 



^ J kn,l' 



zi 

-2 [^0"' 
-2{[^o-^'^^ 

-2 { [^r''^'^]^^, + kk' [^0"''^' 



-1,(1) 



VI 



-1,(1) 



i' J L 



k'n',l 

+ Oie). 

k'n',l ^ ' 



(3.44) 



(3.45) 



Solution of the RG Equations 

The solution of the inhomogeneous RG equation (3.37) for the running Wilson coefficient 
C(//) reads: 



a(5(/x)) = uf{g{ii),g^)Cv{g^) 

J dg'ul!i}{g{i,),g') 



+ 



OlV + -.^ n 'JlV 



Snn' + 



90 



J2 



Skk' — 



Jkk' 



,/2 



-,'2 



167r2 
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,/2 



x i> J 

"vv' 157^2 



^(0) 1 

inm,k' ^ 



+ 



A (0) 



2 ^nm,k' 



/ nm,k' 



167r2 



(3.46) 



Here f/(°) and J are the RG quantities related to single insertions of Q as defined in (3.40) 
and (3.43). Further the Wilson coefficients and evolution matrices have been labeled with 
g{lj,) rather than to comply with the integration variable g'. The coupling at the scale /^o 
at which the initial condition is defined is denoted by go = g{fJ^o)- 

The first term in (3.46) is solely related to matrix elements with single insertions of Q. 
There are no factors involving J here, because the initial coefficients Ci{go) start at the 
order g^. 

(3.46) nicely reveals the structure of the RG in bifocal matrix elements: First the two Wilson 
coefficients C;/ and C,^' run from /iq down to the intermediate scale /i' with g{iJ,') = g' . 
Then they are linked by the anomalous dimension tensor to a single coefficient (7^/, which 
runs further down with the NLO evolution matrix [/(//,//'). The integral sums over all 

intermediate scales ^' . 

For formal analyses as those in chapter 4 the form (3.46) is well suited. In a practical 
calculation, however, the solution of the integral in (3.46) requires to transform the Qis 
to the diagonal basis of 7. The same must be done for the Q^'s of at least one of the two 
operator bases corresponding to the double operator insertion. 

Consider the case that the two operator bases feeding fig. 5.4 (p. 66) have K and K' 
elements, while L linearly independent operators Qi are needed as counterterms to fig. 5.4. 
In the calculation of 773 we will have K = 2, K' = Q and L = 1. The simplest way to 
solve (3.44) is to label the product CkCk' with a single index running from Ito K ■ K'. We 
obtain the form as in (3.37) and can proceed as in the case of single operator insertions. 
The evolution matrix has K ■ K' + L rows and columns. Yet this method hides the different 
sources of RG admixtures to Ci. 

Next one can transform to the diagonal basis for one set of operators Q^'- The RG evolution 
then decomposes into K separate sectors each involving {K' + L) x [K' + L) matrices. 
This is similar to the method used in the LO analysis of K° — K*^ -mixing of Gilman and 
Wise [45]. If this method is used beyond the leading order, however, the mixing matrices of 
different orders must commute, if the standard formalism of sect. 3.8.1 shall apply: Only if 
^(0)^ ^(1) = 0, the NLO matrix J does not induce a mixing of the different sectors of the 
^O evolution. In the calculation of 773, where K — 2, this is indeed the case for the relevant 
LO and NLO mixing matrices given in the literature [43] . 

The discussion above stresses the advantage of commuting LO and NLO mixing matrices, 
especially for the case of multiple operator insertions. It is noticeable that one can indeed 
always achieve this situation by a finite renormalization, because 7(1) is scheme dependent. 
In the context of four-fermion operators scheme transformations involving a set of continuos 
real parameters naturally enter the scene due to the presence of evanescent operators. We 
will look at this in chapter 4. 
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Of course one can further decompose the RG evolution by transforming the remaining 
operator sets to the diagonal basis, too. Yet this does neither significantly improve the 
algorithm nor does it provide more insight into physics. 

The K sectors each having a (fC' + L) x {K' + L) evolution matrix contain some redundancy, 
because all K of them encode the running of the K' + L other operators. This redundancy 
can be removed by replacing the Q's by Cf'^''s which are obtained by dividing the Cis by 
the Cfe's of the diagonahzed basis. This results in a matrix evolution for a coefficient vector 
with K' + K ■ L components. We have used this method in the calculation of rj^ thereby 
reducing the problem to a single 8x8 evolution matrix instead of using four such matrices 
as in [45]. 

We finally remark that the evolution matrix has a block-triangular form, because the 
coefficients Q of the local operators Qi cannot mix into the products • C^/. Since such 
block-triangular matrices appear in many contexts, their discussion deserves a separate 
section. 



3.8.3. Block— triangular Mixing Matrix 

We will frequently encounter block-triangular anomalous dimension matrices: 

7 = 




(3.47) 

with submatrices 7a, 7b and 7c. This form leads to a block-triangular LO evolution matrix: 

C/(°)(//i,/io) = 



/ f/io)(/xi,/xo) 



exp 



4 In 4^ 



exp 



41n4^ 



(3.48) 



with da,dc defined analogously to (3.41). Hence the diagonal blocks of the LO evolution 
matrix exponentiate separately. The RG evolution of the Wilson coefficients of the subspace 
Sa corresponding to 7^ is not affected by the coefficients of the other subspace Sc. The 
Wilson coefficient vectors in Sc, however, satisfy a RG equation to which an inhomogeneous 
term is added involving the solution of the RG equation of subspace Sa. 
This invariance of Sa persists in NLO: (3.47) leads to a block-triangular J matrix in (3.43) 
protecting Sa from admixtures. 

We will be confronted with the case that 7c is a 1 x 1 matrix. Then one can also solve U^^l 
in terms of 7^°^ , the submatrix n x n matrix 7^°^ and the column vector 7^' 



(0). 



USI = db-exp C/f (//i,//o)-l„C/i°n)"i,/^o) 



{da - l„dc) In 



a(/ii) 



(3.49) 



which is a row vector due to d^ — 7^°^^/ (2/3o)- In is the n x n unit matrix. 
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A similar relation for 

" (t I) 

can be derived from (3.43) for the case that Jc is just a number: 

[(1 + 4) In -4]"' (3.50) 

The relations (3.49) and (3.50) are useful to check the RG evolution in the calculation of 
7^3 in chapter 6, because the subspace Sa is related to the Wilson coefficients for |AS| = 1 
transitions which is known to NLO from [44] . 



3.9. Scheme Dependence of Coefficients and Anomalous Dimen- 
sion 

General aspects of renormalization scheme dependences have already been discussed in 
sect. 3.1. In the framework of effective field theories scheme dependences appear in a new 
context: Even when keeping the scheme in the full Standard Model lagrangian fixed there is 
still the freedom to invoke a finite renormalization on the operator renormalization matrix 
(3.24) in analogy to (3.1). Prom now on we will only discuss such scheme dependences 
related to the effective theory. 

Collecting some results of [44] let us write the matrix element in the form: 



9 



167r2 



(Q) 



(0) 



(3.51) 
(3.52) 



If the matrix r obtained in two schemes a and b differs by Ar — Vf, — ra, one finds the 
corresponding finite renormalization as 



Zn — Zb 



9 



+ 0{e) 



This leads to the scheme independence of 7^°^ , while 7*^^^ transforms as 



76 



(1) 



7«+ [Ar,7(°)l +2/3oAr. 



(3.53) 



(3.54) 



This yields for J defined in (3.43): 



(3.55) 



so that J + r-^ is scheme independent. 
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Prom (3.51) also the scheme transformation of the Wilson coefficient is clear: 




(3.56) 



With (3. 38), (3. 42), (3. 56) and (3.55) one finds the scheme dependence related to the upper 
end Ho of the RG evolution canceled in the RG improved Wilson coefficient (3.38). The 
scheme dependence of the lower end h has to cancel with a corresponding dependence in 
the hadronic matrix element. 

Now one is often confronted with the situation that one calculates the matrix element 
(3.51) with different renormalization prescriptions. This is related to the fact that in many 
cases mathematical structures such as the Dirac algebra defined in 4 dimensions cannot 
unambiguosly be extended to D dimensions. This phenomenon is related to the evanescent 
redefinitions mentioned in sect. 3.1: A priori one can modify any D-dimensional quantity 
by terms of order e without changing its limit for D — > 4. One then has to check which 
of the D-dimensional definitions leads to a consistent renormalizable theory complying 
with the observed symmetries of the described physics. Hence one has to prove that two 
different renormalization prescriptions leading to different r's in (3.51) really correspond to 
renormalization schemes. A necessary condition is that the change in the renormalization 
prescription is equivalent to a finite renormalization as in (3.53). This has been done for 
the important question of how to treat 75 in D dimensions in e.g. [43], [44] and [46] and 
will be done in chapter 4 for the question of how to define evanescent operators. 

3.10. Unphysical Operators 

Effective lagrangians also contain unphysical operators. We have already alluded to the 
presence of evanescent operators, which will be analyzed in detail in chapter 4. 

3.10.1. Equation of Motion and BRS Exact Operators 

Another set of unphysical operators is related to the Euler-Lagrange equation of motion 
(EOM) derived from the lagrangian. In classical physics one can simply drop terms which 
vanish by the field equation of motion. In a quantized field theory the issue is more difficult, 
because the path integral involves the integration over all field configurations, not just those 
which obey the EOM. The second subtlety stems from the fact that in Green's functions 
of operators containing derivatives the latter must always be understood to act on the 
time ordering T as well. This is necessary for the Green's function to be a correctly 
defined operator-valued distribution. (Otherwise it could not be Fourier-transformed to 
momentum space giving the familiar —ik'^ for each d'^.) The third difficulty concerns the 
renomalization process, as two operators whose difference vanishes by the EOM may have 
bare Green's functions with a different forest structure. 

The correct implementation of the EOM in quantum field theory also depends on the 
regularization, in BPHZ renormalization the use of the EOM requires oversubtr actions, 
while this is not so for dimensional regularization. 



40 



3. Renormalization Group and Operator Product Expansion 



First analyses of the subject used the EOM derived from the lagrangian (2.15). Yet for 
perturbative calculations the lagrangian (2.36) containing a gauge fixing term and Faddeev- 
Popov (FP) ghosts should be used. In the following we will collect the results obtained with 
this modern approach [51-53,58]. We recommend [54] for a nice explanation of the subject. 
Hence the EOM we use will always include the terms from £gf and Cpp 
The necessary theorems we will be exemphfied with operators of the forthcoming calcula- 
tion of K° — K° -mixing. We will be confronted with the following | AS( = 1 operator: 



Qeom = s-f^LT'^d 



(3.57) 



with L — 1 — 75. This operator with the EOM of the gluon field attached to the |AS| = 1 
current contains the couplings depicted in fig. 5.9 (p. 69), and in the right pictures of 
fig. 5.10 (p. 70) and fig. 5.11 (p. 70). Four-fermion Green's functions with one insertion of 
Qeom vanish [52] . This statement holds for both bare and renormalized Green's functions. If 
external gluons are involved, the matrix element of Qeom equals a sum of matrix elements 
of s^ij,LT"'d in each of which one of the external gluon fields is missing. These contact 
terms vanish after the LSZ reduction (2.25), because they miss a gluon pole. Hence Qeom 
is termed on-shell equivalent to the zero operator. The same, of course, is true for any 
operator vanishing by the EOM for the fermion fields. 

Another class of unphysical operators is related to £fp and the BRS invariance discussed 
in sect. 2.4.2. We will involve the following operator: 

Qbrs = ^^s-f^^Ln- d^d^A^ + s-fi^LT^d- {d''f),)r)Jabc. (3.58) 

Qbrs is the BRS variation of some other operator: 

Qbrs = SbrsQ' = Sbrs (^^s'j^LT'^d ■ 8^7},^ (3.59) 

i.e. it is BRS-exact. From (2.32) one obtains: 

(0|fes(xi)(i(x2)V^(x3)V^(x4)|0) = -(OIQ'-'^^^^'^BRS {s{xi)d{x2mxs)iP{x^)) |0),(3.60) 

which also does not survive LSZ reduction thereby being zero on-shell. 
From the discussion it is clear that Qeom and Qbrs do not contribute to the matching of 
Green's functions with single operator insertions, because the matching can (though need 
not) be done on-shell. 

3.10.2. Mixing 

In order for the coefficients of Qeom and Qbrs to be irrelevant one has to verify the block- 
triangular form of the mixing matrix to ensure that Qeom and Qbrs do not mix into 
physical operators. 
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To discuss the mixing it is useful to introduce three classes of operators: 

P: The set of gauge invariant operators which do not vanish by the EOM. Hence P contains 

physical and evanescent operators. 

S_: The set of operators vanishing by the EOM such as Qeom- 

B: The set of BRS-exact operators. The operator Qbrs belongs to this class. 

Now the renormalization matrix has the following form [52,58]: 



( * 


V 



* 
* 





* \ 

* 
* 



corresponding to 



B 

V ^ / 



(3.61) 



i.e. the renormalized version of (3.60) involves operators from B (B S as counterterms. 
The block-triangular form of (3.61) allows to ignore the operators from £ and B, because 
they neither contribute to the matching nor do their coefficients mix into the coefficients 
of the operators from V. The lagrangian with the renormalized operators of B(BS removed 
is sometimes called reduced lagrangian. In a calculation beyond one-loop, however, the 
operators oi B Q)S appear in subloop counterterm diagrams. On can even avoid operators 
from B altogether by using a background field gauge. In our case of K° — K° -mixing the 
use of a background field gauge does not simplify the calculation. 

Of course there are other non-gauge invariant operators which belong to none of the three 
classes. But since the initial values of their coefficients are zero and none of the operators 
from V ® B ® 8 mixes into them, they clearly do not affect any physical coefficient. 



3.10.3. Double Insertions 

For double insertion the issue is more complex, because now the operators va. B ® £ give 
a non-zero contribution. For a four-fermion matrix element containing one insertion of 
Qeom one finds [52] : 



{jj d''xd''yTQ^OM{x)Q{y)) = 

i E (-l)'^M // d^'xd^yTKls^.Lnix)] -j 



[X 



Q{y)). (3.62) 



Here the sum extends over all gluon field derivatives present in Qeom- -Pa means the number 
of derivatives in A. The main point in (3.62) is that the variation of Qiy) with respect to 
A^(x) is proportional to 5''^\x — y). Hence the matrix element in (3.62) is identical to 
the matrix element of a single local operator Qk{y)- Since this holds for both bare and 
renormalized operators, we can reduce the effective lagrangian C^^ by substituting Qk for 
Qeom- If we kept both of them, the anomalous dimension tensor defined in sect. 3.8.2 would 
have a degenerate eigenvalue. The product of the two Wilson coefficients corresponding to 
the operators of the RHS of (3.62) would evolve parallely to the coefficient Ci of Qi under 
the RG flow, and we could do the reduction of the effective lagrangian at any scale. The 
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method to perform the reduction, i.e. to construct the local operator is decribed e.g. 
in [54]. In a practical calculation, however, it is sufficient to know that we can drop Oeom 
from the operator basis, because we have to do the matching and mixing calculation anyway 
with all operators of the desired dimension in including the QiS. 
Next for the operators from B one finds with (2.32): 

+on-shell vanishing terms. (3.63) 

Here the RHS immediately vanishes, if Q G P © because gauge-invariant operators are 
BRS-invariant. For Q e B we involve 5|j^s on the RHS, which is either zero or yields an 
operator vanishing by the EOM of the ghost field. RecaUing further that the renormahzed 
version of (3.63) involves only operators from B ® £ one finally realizes that one can also 
trade the renormahzed operators from B for suitable QiS. 
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4.1. Motivation 

When effective field theories containing four-fermion interactions are regularized dimen- 
sionally, evanescent operators, which vanish for D = 4, enter the scene. A prominent 
example is connected to the Fierz transformation: Since it is allowed only for D = A, the 
difi^erence of some /^-dimensional operator and its Fierz transformed one is an evanescent 
operator. Since the lagrangian is a D-dimensional object, it also contains these evanescent 
operators. The whole issue would be trivial, if the evanescent sector was decoupled from 
the physical sector in the matching and mixing process. 

Yet when calculating radiative corrections in a SU(N) gauge theory to some matrix element 
of a physical operator one has to face evanescent operators both in the matching procedure 
(3.29) and among the operator counterterms in (3.24). Let us exemplify this with the 
current-current operator Q corresponding to the Dirac structure 

Q = 7m(1-75)®7''(1-75). (4.1) 

We will use this example frequently in this section together with an anticommuting 75 
(NDR scheme). The consistency of this scheme in NLO QCD short distance calculations 
has been demonstrated in [43], [44] and [46] for cases without closed fermion loops. In 
the diagrams of the efi^ective field theory the latter can be avoided by passing to a Fierz 
transformed operator basis. The use of an anticommuting 75 in a renomalizable field theory 
has been justified in [47]. The general arguments in this chapter, however, are not restricted 
to the NDR scheme. 

Figure 4.1: One-loop current-current type radiative corrections to (Qk)- 

Now consider the one-loop matrix elements of Q depicted in fig. 4.1. The Dirac structure 
of the result is easily expressed in terms of a linear combination of Q and an evanescent 
Dirac structure: 

Ei[Q] = 7M7.7^(l-75)«)7V7'^(l-75)-(4 + a£)7^(l-75)(8)7'^(l-75)-(4.2) 

From the basis decomposition of the four-dimensional Dirac algebra we know that (4.2) 
vanishes in four dimensions. Nevertheless (4.2) appears with a factor of 1/s in bare matrix 
elements and therefore in counterterms to physical operators. Depending on the fiavour 
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structure of Q one must take into account also the penguin-type diagrams of fig. 5.10 (p. 
70). They do not contribute to the coefficients of evanescent operators in the one-loop 
order. Hence Q mixes into some evanescent operator After inserting into 

the diagrams of fig. 4.1 another evanescent operator -E2[Q] with five Dirac matrices in 
each string comes into play. Hence one realizes that one must deal with an infinite set 
of evanescent operators, because the Dirac algebra is infinite dimensional in dimensional 
regularization. With the real parameter a in (4.2) we have displayed the arbitrariness in 
the definition of -EifQ]. One can a priori add e times any physical operator (not only 
proportional to Q) to any evanescent operator. Indeed, in the literature one finds different 
definitions of the evanescents. 

Let us display the evanescent operators in the effective lagrangian (3.26) in the following 
way: 

rl r< y-l/nbarc I ^ /v— 1 771 \r\ ibare 

+CE,,Z-^],iQ''r + Ce,,Z-^],eA [Qt''' . (4.3) 

Here and in the following we will distinguish the renormalization constants related to some 
evanescent operator Ej[Q„l\ by denoting the corresponding index with Ejm- 
Buras and Weisz were the first to recognize that evanescents cannot simply be neglected but 
infiuence physical quantities [43]. They focused on the evanescent operators' contribution 
to the matching (3.29) and found that the one-loop matrix elements of Ei[Q] involve finite 
components proportional to physical operators. In [43] it has been proposed to cancel these 
components, which are local, by a finite renormalization of -E'i[(5]- By this the evanescent 
operators do not contribute to the matching in (3.29) anymore and their undetermined 
coefficients Cg^^ are irrelevant at the matching scale Hq. From (3.34) the authors of [43] 
found the influence of this flnite renormalization on the physical part of the NLO anomalous 
dimension matrix 7^-'^) . 

To insure that evanescent operators remain irrelevant at any scale one must insure that the 
evanescents do not mix into physical operators. As discussed in sect. 3.8.3 such a block- 
triangular anomalous dimension matrix protects the physical coefficients Ck in (4.3) from 
admixtures of the Ce^^- This has been noticed first by Dugan and Grinstein in [48]. They, 
however, have used a definition of the evanescents different from that of [43], i.e. a different 
a in (4.2). We do not need to repeat the definitions used in [43] and [48] in this introductory 
section to understand the issue. The authors of [48] proved that the anomalous dimension 
matrix indeed has the desired block-triangular form using their special definition of the 
evanescents, if the finite renormahzation of [43] is performed. Yet it is not clear at all that 
these features holds for any definition of the evanescents: An evanescent redefinition of e.g. 
Ei[Q] in (4.2) by adding e ■ 5a ■ Q to it modifies the components of the divergent parts 
of the bare two-loop matrix elements proportional to physical operators. Hence we are 
immediately confronted with the question: 

Does the definition of the evanescents employed in [43] also yield a block-triangular anoma- 
lous dimension matrix? 

We will answer this question affirmatively in sect. 4.3. There we will also generalize the 
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method of [43] . Especially one is not forced to use the definition of the evanescent operators 
proposed in [48], whose implementation is quite cumbersome. 
Then the next point to discuss is: 

Is the dependence on a in (4-2) spurious or does it influence the physical Wilson coefficients 
Ck and the physical portion of the anomalous dimension matrix 7? 

We will find in sect. 4.4 that the latter is the case. Here we will also find that evanes- 
cent redefinitions corresponding to a change of a in (4.2) imply a renormalization scheme 
dependence in the physical sector, i.e. we will be able to find a quantitity Ar such that 
the Wilson coefficients and the mixing matrix transforms as in (3.56) and (3.54). This 
touches a non-trivial aspect: It is not clear from the beginning that the change of the 
definition of the evanescents does not spoil renormalizability of the theory, i.e. the locality 
of the divergences. Yet in NLO it is trivial that renormalizability holds for any definition 
of the EjkS, because a redefiniton as in (4.2) only reshuffles local terms in the divergences 
of the bare matrix elements. When going beyond NLO terms of order e"^, etc., in (4.2) 
become important, and these terms might be constrained by the condition of maintain- 
ing renormalizability. We will not discuss this point. In fact we will extend the discussion 
only beyond NLO when addressing the all-order argument of [48] for the block-triangular 
mixing matrix. 

A main subject of this thesis is the correct renormalization of bifocal structures with two 
coloured operators. Clearly for this we have to answer the question: 

What is the correct method to handle evanescent operators in matrix elements with two 
operator insertions? 

This point will be the subject of sections 4.5 and 4.6, where four-fermion processes and 
inclusive decays will be discussed. 

4.2. Preliminaries and Notation 

Let us consider a set of physical four-quark operators {Qfc = V'^feV' " '^Qki', k — 1,2,3, . . .} 
as shown in (3.23). For simplicity we restrict the basis to dimension-six operators. As 
usual we are interested in the Green's functions of a SU(N) gauge theory with insertions 
of Qk renormalized in a mass-independent scheme. For simplicity of the notation we pick 
minimal subtraction (MS). The arguments are easily generalized to other schemes like MS. 
The dirac structures Qk — Qk Qk corresponding to Qk are considered to form a basis of 
the space of Lorentz singlets and pseudosinglets ior D — A. Neither the Lorentz indices 
of Qk and Qk are displayed nor any flavour or colour indices, which are irrelevant for the 
discussion of the subject. [F (g) 1] Qk [1 (X> F'] means Tqk ® g^F'. 

For the purpose of this chapter we define the perturbative expansion of the matrix elements 
of Qk differently from (3.51): 




(4.4) 
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4. Evanescent Operators 



We have extracted the wave function renormahzation here, because it is not related to the 
issue of the mixing of evanescent operators. 

Now the insertion of Qk into the one- loop diagrams of fig. 4.1 yields a linear combination 
of the Qis and a new operator with the Dirac structure Q'j^ = Yip^a ® 1] Qk W ® 7'^7^]: 

&rt'' = 4?(Q.)^°^ + 4!i',(Q'.)^°^ no sum on A;, (4.5) 

where (. . .)*^°^ denote tree level matrix elements. Both coefficients have a term proportional 
to 1/e and a finite part. Q'^ is now decomposed into a linear combination of the Q^'s and 
an evanescent operator: 

= {hi + auie)QY''-rE^{Qk^-'''-rO{e'). (4.6) 

Here the /jt/'s are uniquely determined by the Dirac basis decomposition in D = 4 dimen- 
sions. The flfei's, however, are arbitrary, and a different choice for the a^/'s corresponds to a 
different definition of -EifQfc] = -^/iIQa;, {o-rs}]- When going beyond the one-loop order new 
evanescent operators £^2[Qfc], -E'alQfc], • • • will appear. Their precise definition is irrelevant 
for the moment and will be given after (4.16). 

Now in the framework of dimensional regularization the renormahzation of some physical 
operator requires counterterms proportional to physical and evanescent operators: Prom 
(4.3) we find 

Zl&D = ZMT) + Zk,E,^{^%{Q^r-) . (4.7) 

(4.5) and (4.6) imply that depends on the a^^'s, while zj^^^^^ is independent of them. 
The definition of the coefficients in the expansion of Z in terms of the gauge coupling 
constant g has been given in (2.35). 

Next we describe the method of Buras and Weisz [43] for the treatment of evanescent 
operators: They have determined the Ofc^'s by choosing some set of Dirac structures M = 
1^(1) (g) . . . ^(10) (g) which forms a basis for D = A, and contracting all elements 

in M with Q'^, and Qi in (4.6): 

tr (7^"^^g^7^"^Vj = ifki + aMs) tr {^^'^^qn^'^^qi) + 0{e'), 

no sum on k and on m—l,. . . , 10. (4-8) 

The solution of the equations (4.8) uniquely defines the fki + aus. In other words, Ei[Qk] 
obeys the equations: 

^i[g.].,..7ir^7]r^ = 0{e') form=l, ...,10, (4.9) 

where i,j,r, s are Dirac indices. 

Our arguments will not depend on the scheme used for the treatment of 75. The use of 
a totally anticommuting 75 does not cause any ambiguity in the trace operation in (4.8), 
because all Lorentz indices are contracted, so that the traced Dirac string is a linear 
combination of 75 and the unit matrix. 
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E.g. the choice of 
M 



{1 1, 1 (8) 75, 75 (g) 1, 75 (8) 75, 7^ ® j^, 7^ (8) 7^75, 



l5lu. ® 7", 757m ® 7''75, 



(7'^'^} (4.10) 
gives for Q in (4.1) 

Q' = 7p7a7M (1 - 75) «) 7^7^^ (1 - 75) = (4 - 8e)Q + Ei[Q] + ©(s^) (4.11) 

as in [43]. We remark that this choice a = —8 respects the Fierz symmetry, which relates 
the first to the second diagram in fig. 4.1. 

A basis different from M in (4.10) yields the same fk/s, but different Ofc^'s. For example by 
replacing the sixth and eighth element of M in (4.10) by ^alpls ® ^""^^^^ and 75707^75 ® 
7"7^7'^ one finds 

Q' = 4Q + 16^(1 + 75) 0(1- 75) + £^([Qfc] + C'(£2), 

instead of (4.11), i.e. a difi^erent evanescent operator. The Dirac algebra is infinite dimen- 
sional for non- integer D and is spanned by M and an infinite set of evanescent Dirac 
structures. Hence one can reverse the above procedure and first arbitrarily choose the a^i's 
and then add properly adjusted linear combinations of the evanescent structures to the 
elements of M such as to obtain the chosen a^^'s. 

Yet the so defined evanescent operators do not decouple from the physics in four dimensions: 
In [43] it has been observed that their one-loop matrix elements generally have nonvanishing 
components proportional to the physical operators Q^. 



^0 
1 

+- 

£ 



(1) 



1 



mQk]) 



(0) 



(1) 



E-iu.Eo 



(4.12) 



Here a second evanescent operator E2, which will be discussed in a moment, has appeared. 
Clearly no sum on k is understood in (4.12) and in following analogous places. In (4.12) 
[Zo^^]e-i^^,i is local, because it originates from the local 1/e-pole of the tensor integrals and 
a term proportional to e stemming from the evanescent Dirac algebra. For the same reason 
there is no divergence in the term proportional to (Qi)^^^- Now in [43] it has been proposed 
to renormalize Ei by a finite amount to cancel this component: 



EilQk] 



ibare 



Ei[Qk]"- + 



167r2 



(1) 



E\k,E-[ 



EilQk 



Eik,l 



Qi 



(4.13) 
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4. Evanescent Operators 



With (4.13) the renormahzed matrix elements of the evanescent operators are 0{e), so that 
they do not contribute to the one-loop matching of some Green's function C^^^ in the full 
renormalizable theory with matrix elements in the effective theory: 

iG- = Ci{Qiy^^^ + CE,AEi[Qk]y'''' + 0{g'), (4.14) 

i.e. the coefficients C^;^^. are irrelevant, because they multiply matrix elements which vanish 
for D = 4. In [43] it has been further noticed that Zq'^ in (4.13) influences the two-loop 
anomalous dimension matrix of the physical operators, so that the presence of evanescent 
operators indeed has an impact on physical observables. 

Next we discuss E2[Qk\, which has entered the scene in (4.12): When inserting Ei[Qk\ 
defined in (4.6) into the one-loop diagrams of fig. 4.1, one involves 

= [/ + ae]l^ Qi + {hi + aus) [Qi] 

+ [7^7.®l]Ei[Q,][l®7V] (4.15) 
= {[f + ae]li + hie]Qi + E2[Qk]+0(e''), (4.16) 

which defines E2[Qk] — E2[Qk,{0'rs},{brs}]- Only the last term in (4.15) can contribute 
to the new coefficients hki- If the projection is performed with e.g. M defined in (4.10), 
one finds hki = 0^. In our discussion we will keep hki arbitrary. Clearly, one has a priori to 
deal with the mixing of an infinite set of evanescent operators for each physical 

operator Q^-, where denotes the new evanescent operator appearing first in the 

one-loop matrix elements of Ej[Qk]. 

With the finite renormahzation of E-i[Qk\ in (4.13) the evanescent operators do not affect 
the physics at the matching scale, at which (4.14) holds. In order that this will be true at 
any scale ^, however, one must also ensure that the evanescent operators do not mix into 
the physical ones, as noticed first by Dugan and Grinstein in [48]. In their analysis they 
have introduced another way to define the evanescent operators, which is also frequently 
used: It is easy to see that one can restrict the operator basis {Qk\ to the set of operators 
whose Dirac structures qk, qk are completely antisymmetric in their Lorentz indices. Dirac 
strings being antisymmetric in more than four indices vanish in four dimensions and are 
therefore evanescent. Operators with five antisymmetrized indices correspond to Ei in our 
notation, and E2 would be expressed in terms of a linear combination of Dirac structures 
with seven and with five antisymmetrized indices. Clearly this method also corresponds 
to some special choice for the a^^'s and h^is in (4.6) and (4.16). Now in [48] the authors 
have proven that with the use of those definitions and a finite renormahzation analogous 
to (4.13) the anomalous dimension matrix indeed has the desired block-triangular form, 
so that the evanescent operators do not mix into the physical ones. While the anomalous 
dimension matrix is trivially block-triangular at one-loop level, the proof for the two-loop 

^This is the case for any basis M in which for each 7^'") ® 7^"*^ e M the quantity ^p^crl^^^^l'^l'' ® 7^™^ 
is a Unear combination of the elements in M. 
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level was given in [48] by the use of the abovementioned special definition of the evanescent 
operators. The latter, however, has some very special features, which arc absent for the 
general case with arbitrary a^/'s and bki^s, e.g. the definition used in [48] automatically 
yields an anomalous dimension matrix which is tridiagonal in the evanescent sector. 
This immediately raises the question whether the more general method of [43] also yields 
a block-triangular anomalous dimension matrix for any chosen projection basis, i.e. any 
chosen set {aki}- In the following section we will prove that this is indeed the case and, 
more generally, that one may also choose the bki^s in (4.16) completely arbitrary. 



4.3. Block Triangular Anomalous Dimension Matrix 



Consider some set of physical operators {Qk} which closes under renormalization together 
with the corresponding evanescent operators {Ej[Qk] : j > 1}. Their 0(£)-parts a^s, brs, ■ ■ ■ 
are chosen arbitrarily. We want to show that the block of the anomalous dimension matrix 
describing the mixing of Ej[Qk\ into Qi equals zero. 



0, 



(4.17) 



provided one uses the finite renormalization described in (4.13). 

Our sketch will follow the outline of [48], where (4.17) has been proven by complete induc- 
tion. At the one- loop level (4.17) is trivial, and the induction starts in two-loop order: The 
next-to-leading order contribution to the anomalous dimension matrix has been stated in 
(3.34). The nonzero contributions to (4.17) in two-loop order are 



7 



(1) 



-4 



(2) 



Ejk,l 



-2f3o 



Ejk,l 



+2 Zl 



.(1) 



Ers J 



+ 



Ejk,m 



zi 



(1) 



ml 



(4.18) 



Here (4.18) contains terms which arc absent when the special definition of the evanescent 



operators in [48] is used: In [48] one has Z^^^ 



for j > 2 contrary to the general case, 



where any evanescent operator can have counterterms proportional to physical operators. 



Next we look at 



Z[ 



(2) 



, which stems from the 1/e-term of the 0(5f^)-matrix elements 

Ejk,l 

of Ej[Qk]- As discussed in [48], these 1/e-terms originate from 1/e^-poles in the tensor 
integrals multiplying a factor proportional to e stemming from the evanescent Dirac alge- 
bra. Now in each two-loop diagram the former are related to the corresponding one-loop 
counterterm diagrams by a factor of 1/2, because the non-local 1/s-poles cancel in their 
sum [42]. For this it is crucial that the one-loop counterterms are properly adjusted, i.e. 
that they cancel the 1/e-poles in the one-loop tensor integrals. In the one-loop matrix el- 
ements of evanescent operators the latter are multiplied with e originating from the Dirac 
algebra. Hence the proper one-loop renormalization of the evanescent operators must be 
such as to give matrix elements of order £, as shown for Ei[Qk\ in (4.13). 
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4. Evanescent Operators 



Prom the one-loop counterterm graphs one simply reads off: 

which yields the desired result when inserted into (4.18). Here the first two terms stem 
from insertions of physical and evanescent counterterms to Ej[Qk\, while the term involv- 
ing the coefficient of the one-loop /^-function f3{g) = — gf^/(167r^)/?o originates from the 
diagrams with coupling constant counterterms. The terms involving the wave function 
renormalization constants cancel with those stemming from the factor in (4.7). 
The inductive step in [48] proving (4.17) to any loop order does not use any special defi- 
nition of the evanescent and therefore applies unchanged here. 



4.4. Evanescent Scheme Dependences 



In this section we will analyze the dependence of the physical part of 7^^^ given in (3.34) 
and of the one-loop Wilson coefficients on an and bu. In practical next-to-leading order cal- 
culations one often has to combine Wilson coefficients and anomalous dimension matrices 
obtained with different definitions of the evanescent operators and it is therefore important 
to have formulae allowing to switch between them (see e.g. appendix B of [49]). 
We start with the investigation of the dependence of 7^^^ on an. The bare one- loop matrix 
element 



{Qk 



bare\ (1) 



1 r 



^Z\ 
+0{e) 



(1) 



kj 



+ 



4 



(1) 



kj 



zi 



(1) 



{EilQk]) 



(0) 



(4.19) 



is independent of an, which is evident from (4.5). Ei[Qk] depends linearly on an through 
its definition (4.6) with the coefficient 



_d_ 

dai 



Ei[Qk\ — —sSkiQi, 



(4.20) 



so that (4.19) gives: 



d 
dan 



d, 



(1) 



kj 



(1) 



^ki^lj, 



(4.21) 



while z[^^ is independent of an- 

(2) 

In the same way on can obtain the afei-dependence of Z{ . (4.7) reads to two-loop order 
(cf. (4.4)): 



(0) 



(4.22) 
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Prom (4.16) we know 



_d_ 



MQr 



and from (4.19) one reads off: 



(1) 



(4.23) 



(4.24) 



These relations and (4.21) aUow to calculate the derivative of (4.22) with respect to an. 
Keeping in mind that the evanescent matrix elements are 0{e) the 0(l/£:)-part of the 
derivative yields: 



d 
dan 



(2) 



kj 



+ 



ki 



k,E2r. 



i,Eu 



k,Eu 



{Smifij + fmiSij) , no sum on i. 



(4.25) 



Wc pause here to stress that (4.25) is local showing that an evanescent redefinition does 
not spoil the renormalizability of the physical sector in NLO. 



Next we extract 
ceding section: 



.(2) 



^(2) 



^(2) 



k,Eu 



k,E2-n 



from the one-loop counterterm diagrams as described in the pre- 



ki 



.(1) 



Z\ 



(1) 



z\ 



(1) 



1 

2 



(1) 



Ski, 



i,Eu 

no sum on i 



EiitEii 



Jki 







z?' 




k,Eik 





Elk,E2k 



^km-i 



no sum on k. 



(4.26) 



After inserting (4.26) into (4.25) we want to substitute the last term in (4.25). For this we 
derive both sides of (4.12) with respect to an giving: 



Zl 



(1) 



Elk,E2k 



i^kiflj + fki^lj) — 



_d_ 

da:, 



Eik,j 



z\ 



(1) 



I] 



Ski — 



z\ 



(1) 



Elk, Elk 



SkiSij, no sum on k. (4.27) 



Finally one has to insert the expression for (4.25) obtained by the described substitutions 
into 



_d_ 

dai, 



7 



(1) 



-4- 



d 



Z\ 



(2) 



kj 



+ 2 



' dan 
which follows from (3.34). The result reads: 
d 

kj 



z\ 



(1) 



dan 



-2 



Z'^\. z\ 



(1) 



_d_ 

k,Eik dai 



4 



'(1) 



Eik,j ' 



no sum on k, 



+2/3o Zl 



(1) 



i,Eu 



i,Eu 
Ski^lj, 



i,Eu 



no sum on i. 



(4.28) 
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4. Evanescent Operators 



Since the quantities on the right hand side of (4.28) do not depend on an, one can easily 
integrate (4.28) to find the desired relation between two 7's corresponding to different 
choices for aki in (4.6). To write the result in matrix form we recall the expression for the 
physical one-loop anomalous dimension matrix 



(0) 



7 



and introduce the diagonal matrix D with 



-2 



(1) 



Hence 



^W{a') = 7^(0)+ D-(a'-a),7W +2PoD-{a'-a 



Dim 



(1) 



-1 i,Exi 



no sum on i. 



(0) 



(4.29) 



(4.30) 



where the summation in the row and column indices only runs over the physical subma- 
trices. 

(4.30) exhibits the familiar structure of the scheme dependence of 7*^^'* as given in (3.54). 
The matrix J in (3.55) consequently reads: 



J(a') = J{a)-[D-{a' -a)f 



(4.31) 



Hence we have found a scheme dependence which is caused by an evanescent redefinition 
of a hare operator basis (i.e. of the bare evanescent operators). 

The dependence of the one-loop matrix elements on a can be found easily from (4.24) and 
(4.21): 



(g)--(a') = 



1 + 



167r2 



D-ia' 



{QY^-{a) + 0{g% 



(4.32) 



Since in (4.14) G does not depend on a and the evanescent matrix clement is 0(e), the 
corresponding relation for the Wilson coefficients at the matching scale reads: 



9' 



167r2 



D • {a — a) 



(4.33) 



Thus also the Wilson coefficient transforms correctly as in (3.56) and the scheme depen- 
dence cancels between coefficients and anomalous dimension matrix. 

In the same way one can investigate the dependence of 7^^) on hu given in (4.16): While 
z'l^ and z'q^ depend on hu, this dependence cancels in (3.34). Hence neither 7*^^^ nor the 
one-loop Wilson coefficient are affected by the choice of hu. Thus in NLO the hus do not 
induce a scheme dependence. By keeping hu arbitrary in a practical calculation one has a 
non-trivial check of the correct handhng of the evanescents: The has disappear from the 
result only after inserting the evanescent counterterms with a factor of 1/2. 
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The nice feature of (4.32), (4.33) and (4.30) is that the summed matrix indices only run 
over the physical subset. Prom (3.54) and (3.56) we could have expected a formula where 
the summation runs over the whole operator basis including the cvancscents. This is why we 
could not simply deduce (4.30) from (4.32) using (3.53) and (3.54). Possible contributions 
from summations over evanescent operator indices cannot be inferred from (4.32), because 
there they would multiply vanishing matrix elements. 

We have already stressed in the end of sect. 3.8.1 that it is advantageous to have a 7*^^^ 
which commutes with 7^*^). Now one can use (4.30) to simplify 7^^): By going to the diagonal 
basis for 7*^°^ = diag (li^^^ one can easily find solutions for a' — a in (4.30) which even 
W(n'\ = provided that all Zk^Ei^^s are nonzero and all eigenvalues of 7^°-* satisfy 
^ 2j3o. We will exemplify this in a moment. 
A choice for an which leads to a 7*^^) commuting with 7^°^ has been done implicitly in [43]: 
There the mixing of the two operators (5+ = Q (l + l) /2 and Q_ = Q (j. — Tj /2 has 

been considered, where 1 and 1 denote colour singlet and antisinglet and Q was introduced 
in (4.1). Now (5+ is self-conjugate under the Pierz transformation, while Q- is anti-self- 
conjugate, so that 7^°) is diagonal to maintain the Pierz symmetry in the leading order 
renormalization group evolution. As remarked after (4.11), the definition of in (4.11) 

is necessary to ensure the Pierz symmetry in the one-loop matrix elements. Consequently 
with (4.11) also 7*^^^ has to obey the Pierz symmetry preventing the mixing of Q+ and 
i.e. yielding a diagonal 7''^-*. A different definition of Ei[Q] would result in non-Pierz- 
symmetric matrix elements, but in renormalization scheme independent expressions they 
would combine with a non-diagonal 7*^^^ such as to restore Pierz symmetry. 
Let us consider the example above to demonstrate that one can pick a' such that 7(^^(0') = 
0: In [43] the definitions 

Ei[g±] = (±1-1) [7,7.7,(1 -75)® 7"7V(i- 75) 

-{4-8e) 7:^(1 -75)® 7^(1-75)], 

i.e. a++ = 8(l/iV — 1), a = 8(1/A'"-|- 1), a+_ = a_+ = 0, were adopted yielding a diagonal 

7^-^^ (a) —dia,g(j^\a) , {aj^ . Prom the insertion of Q+ and Q- into the diagrams of 

fig. 4.1 one finds Z^^Ei+ = Z_^Ei- = 1/4. Hence if we pick a'^^ — a±± — 2/^2\a)/(3o and 
a'^^ = 0, (4.30) will imply -f^^^{a') = 0. 



4.5. Double Insertions 



4.5.1. Motivation 

Prom the preceding sections we know that the coefficients Ce^i^ of properly renormalized 
evanescents in (4.3) are irrelevant for physical quantities and remain undetermined. In the 
following we will extend these results to 4-fermion Green's functions with two insertions of 
local operators. 
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4. Evanescent Operators 



For the discussion of Green's functions with insertion of two local operators as depicted in 
fig. 5.4 we first display the evanescents in in (3.28): 



bare 



bare 



+CEjkCE.,y l^B.ifeE^,,,,;^^''' + ZE-kE.,y,Eri^r [Qi 

+CkZ^l'Q\'^'"' + CkZ^Eri^r Ql 
+'^Ejk^EjkM''"' + ^E^k^E^kEri^r Ql 



(4.34) 



bare 



bare 



(4.35) 



Since the Q^'s have dimension six, the Qi^s are dimension-eight operators. For simplicity, 
we assume the Qi^s to be linearly independent from the QkS, e.g. the QkS represent 
AS" = 1 operators and the Q/'s denote AS — 2 operators. The Er [Qz] in (4.35) are defined 

analogously to (4.6) with new coefficients fki, aki, bki, etc. Hence new arbitrary constants 
Sfcz, bki potentially causing scheme dependences enter the scene. 
Clearly the following questions arise here: 

1. Are the coefficient functions C^^.^ irrelevant also for the double insertions; i.e. do 

(^J TEQ^ and (^J T EE^^ (4.36) 

contribute to the matching procedure and the operator mixing? 

2. Does one need a finite renormalization in the evanescent sector of double insertions; 
if yes, how does this affect the anomalous dimension tensor? 

3. How do the Ci and anomalous dimension matrices depend on the Uki, bki, o,ki, bki ? 

4. Are the RG improved observables scheme independent? 



4.5.2. Scheme Consistency 

In this section we will carry out the program of section 4.3 for the case of double insertions 
to answer questions 1 and 2 of page 54. 

We have mentioned already in sect. 3.7.2 that two cases have to be distinguished: The 
matrix clement of the double insertion of the two local renormalized operators can be 
divergent or finite: 
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Therefore we need or do not need a separate renormalization for the double insertion 



^kk' I 








case 1 
case 2 



(4.38) 



Let us start the discussion with the matching procedure 3.30: 



T^QkQk' 



+ CkCE,,., 



T QkEii [Qy 
CilQi, 



Qi 



(4.39) 



where G"^^° corresponds e.g. to the box-function shown in fig. 5.1 and the radiative correc- 
tions depicted in fig. 5.2. 

Since the coefficients Ce^^. must be irrelevant for this matching procedure, one must have 



T Ej [Qk] Qk' 



' O (e°) in case 1 (LO) 
O {s}) in case 1 (NLO and higher) 
O ie^) in case 2 



(4.40) 



and the analogous relation for two insertions of evanescent operators. Recall from sect. 
3.7.2 that in case 1 the LO matching is performed by the comparison of the coefficients of 
logarithms of the full theory amplitude and the effective theory matrix clement (4.40) (the 
latter being trivially related to the coefficient of the divergence), while the NLO matching 
is obtained from the finite part and also involves the matrix elements of the local operators. 
In case 2 the matching is performed with the finite parts in all orders. In both cases the 
condition (4.40) is trivially fulfilled in LO, because the evanescent Dirac algebra gives an 
additional e compared to the case of the insertion of two physical operators. Therefore a 
finite renormalization for the double insertion turns out to be unnecessary at the LO level. 
This statement remains valid at the NLO level only in case 2, in case 1 condition (4.40) no 
longer holds if one only subtracts the divergent terms in the matrix elements containing a 
double insertion. With the argumentation preceding (4.13) one finds that in this case the 
finite term needed to satisfy the condition (4.40) is local and therefore can be provided by 
a finite counterterm. 

The operator mixing can be analyzed analogously to section 4.3. Using the finite renor- 
mafization enforcing (4.40) and the locafity of counterterms, one shows for the anomalous 
dimension tensor (3.45): 



= 



and 



(0) (1) 

Egi,n 



0, 



(4.41) 



i.e. a double insertion containing at least one evanescent operator does not mix into physical 
operators. Together with the statement that evanescent operators do not contribute to the 
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matching this proves our method to be consistent at the NLO leveL As in the case of single 
insertions one can pick the aki, bki,. . . completely arbitrary and then has to perform a finite 
renormalization for the double insertions containing an evanescent operator in (4.36). This 
statement remains valid also in higher orders of the SU(N) interaction, which can be proven 
analogously to the proof given by Dugan and Grinstein [48] for the case of single insertions. 
Now we use the findings above to show the nonvanishing terms in (3.45) explicitly for the 
physical submatrix: 



/ kn.l 



■1,(2) 



kn,l 
1,(1)1 



z\ 

^0 



■1,(1) 



1,(1) 



-1,(1) 



- 2 



Zl 



Eii'l 
1,(1) 



nE, 



-1,(1) 



(4.42) 



The last equation encodes the following rule for the correct treatment of evanescent op- 
erators in NLO calculations: The correct contribution of evanescent operators to the NLO 
physical anomalous dimension tensor is obtained by inserting the evanescent one-loop coun- 
terterms with a factor of | instead of 1 into the counterterm graphs. Hence the finding 
of [43] for a single operator insertion generalizes to Green's functions with double inser- 
tions. Here the second term in (4.42) corresponds to the graphs with the insertion of a 
local evanescent counterterm into the graphs depicted in fig. 4.1, while the last to terms 
correspond to the diagrams of fig. 5.4 with one physical and one evanescent operator. 



4.5.3. Double Insertions: Evanescent Scheme Dependences 

In this section we will answer questions 3 and 4 on page 54. Let us first look at the 
dependence of the anomalous dimension tensor ^kk',i on the coefficients a^g. First one 
notices, that the LO 7^^,°! j is independent of the choice of the a^s- Similarly to the procedure 
of section 4.4 one derives the following NLO relation: 

tS,. {a') = tS,. (a) + [D ■ {a' - a)],, + [D ■ {a' - a)],,, (4.43) 

with the diagonal matrix D from (4.29). Note again that the indices only run over the 
physical subspace. 

The variation of the anomalous dimension tensor ^kk',i with the coefficients a^s again van- 
ishes in LO, in NLO we find the transformation 
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(0) 



(4.44) 



As in the case of single insertions, up to the NLO level there exists no dependence of 7 on the 
coefficients brs and also no one on the brs ■ This provides a nontrivial check of the treatment 
of evanescent operators in a practical calculation, when the brs, brs are kept arbitrary: the 
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individual renormalization factors Z each exhibit a dependence on the coefficients hrsi Ks 
but all this dependence cancels, when the Z's get combined to 7. This has been done in 
the calculation of 773 in chapter 6. 

Next we verify the scheme independence of RG improved physical observables. From the 
solution of the inhomogeneous RG equation (3.46) we find the local operators' Wilson 
coefficients Ci independent under the transformations (4.33), (4.31) and (4.43). 
In a similar way one treats the scheme dependence stemming from the coefficients aki- Here 
some work has been necessary to prove the cancellation of the scheme dependence connected 
to g''^Jkk' and ^^^^ y in (3.46): Although it is not possible to perform the integration in (3.46) 
without transforming some of the operators to the diagonal basis, one can do the integral 
for the scheme dependent part of (3.46), because the part of the integrand depending on 
cikis is a total derivative with respect to g. There is one important difference compared to 
the dependence on the a^/'s: A scheme dependence of the Wilson coefficient related to the 
lower end of the RG evolution ^ in (3.46) still remains. This residual ciki dependence must 
be canceled by a corresponding one in the hadronic matrix element. If the matrix elements 
are obtained in the parton model, the dependence of the Q's on Uki cancels due to (4.32). 
Finally, as in the case of single insertions [44] , one can define a scheme-independent Wilson 
coefficient for the local operator 



Q' + ^ • rnm,i ■ C^^ in) in) + O (g^) , 

(4.45) 



which multiplies a scheme independent matrix element defined accordingly. It contains the 
analogue of f in (3.52) for the double insertion 



/ - \ (0) 

TQnQm) = ^■'^nm,l-{Ql) . (4-46) 



4.6. Inclusive Decays 

Inclusive decays are calculated either by calculating the exclusive process and performing 
a subsequent phase space integration and a summation over final polarizations etc. (ref- 
erenced as method 1) or by use of the optical theorem, which corresponds to taking the 
imaginary part of the self-energy diagram depicted in fig. 4.2 (method 2). This figure shows 
that inclusive decays are in fact related to double insertions, but in contrast to the case 
of section 4.5 they do not involve local four-quark operators as counterterms for double 
insertions. In fact, even local two-quark operator counterterms would only be needed to 
renormalize the real part, but the imaginary part of their matrix elements clearly vanishes. 
The only scheme dependence to be discussed is therefore the one associated with the a^/'s, 
bki^s, etc., as there are no a^i^s b^i^s, etc. involved. 

To discuss the dependence on the a^j's it is nevertheless advantageous to consider method 1, 
i.e. the exclusive process plus the subsequent phase space integration. From section 4.4 we 
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already know most of the properties of the exclusive process: At the upper renormalization 
scale the properly renormalized evanescent operators do not contribute and the scheme 
dependence cancels. Further we know from (4.31) the scheme dependence of the (RG im- 
proved) Wilson coefficients at the lower renormalization scale. What we are left with is 
the calculation of the properly renormalized operators in perturbation theory, i.e. with on- 
shell external momenta. Clearly the form of the external states does not affect the scheme 
dependent terms of the matrix elements, they are again given by (4.32) and therefore can- 
cel trivially between the Wilson coefficients and the matrix elements, because the scheme 
dependent terms are independent of the external momenta. Since we now have a finite 
amplitude which is scheme independent, we may continue the calculation in four dimen- 
sions and may therefore forget about the evanescent operators. The remaining phase space 
integration and summation over final polarizations does not introduce any new scheme 
dependence, therefore we end up with a rate independent of the a^/'s, bki^s, etc. ^ 
Alternatively one may use the approach via the optical theorem (method 2). Then one 
has to calculate the imaginary parts of the diagram in fig. 4.2 plus gluonic corrections. Of 

Figure 4.2: The lowest order self-energy diagram needed for the calculation of inclusive 
decays via the optical theorem (method 2). 

course the properly renormalized operators have to be plugged in: 



One immediately ends up with a finite rate. What is left to show is the consistency of 
the optical theorem with the presence of evanescent operators and with their arbitrary 
definition proposed in (4.6), (4.16): The result must not depend on whether the evanescent 
operators are kept in the basis inserted into the diagram of fig. 4.2 or whether they have 
been removed knowing the above result of method 1. This means that evanescent opera- 
tors must not contribute to the rate, i.e. diagrams containing an insertion of one or two 
evanescent operators must be of order e 

Im (Ei [OkY'"' or) = O (£) and Im (e, [0^]^^^^ Ej [O,]"-^") ^ O {e) . (4.48) 

As in the previous sections one can discuss tensor integrals and Dirac algebra separately 
leading to (4.48). 

4.7. Summary and Outlook 

Let us fist the results of this chapter: 

^We discard problems due to infrared singularities and the Bloch-Nordsiek theorem. At least in NLO 
one can use a gluon mass, because no three-gluon vertex contributes to the relevant diagrams 




(4.47) 
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i) It is allowed to redefine any bare evanescent operator by {D — 4) times any physical 
operator without affecting the block-triangular form of the anomalous dimension 
matrix, which ensures that properly renormalized evanescent operators (as described 
in [43]) do not mix into physical ones. 

ii) We have analyzed the renormalization scheme dependence associated with the redefi- 
nition transformation in the next-to-leading order in renormalization group improved 
perturbation theory. The formulas transforming between different schemes have been 
derived. It is meaningless to give some anomalous dimension matrix or some Wil- 
son coefficients beyond leading logarithms without specifying the definition of the 
evanescent operators used during the calculation. In physical observables, however, 
this renormalization scheme dependence cancels between Wilson coefficients and the 
anomalous dimension matrix. One may take advantage of this feature by defining the 
evanescent operators such as to achieve a simple form for the anomalous dimension 
matrix. 

iii) We have extended the work of [43] and [48] and the results of i) and ii) to the 
case of Green's functions with two operator insertions. This analysis is necessary for 
the correct treatment of evanescents in particle-antiparticle mixing and rare decays. 
Scheme-independent Wilson coefficients have been defined. 

iv) Finally we have analyzed the role of evanescents in inclusive decay rates. 

The issue of the renormalization of composite operators in dimensional regularization is 
not restricted to four-quark operators. We suggest that one may investigate to what extent 
evanescent redefinitions are useful to gain insight into the renormalization scheme depen- 
dences of other operators. If, for example, some current j is handled with two renormaliza- 
tion prescriptions denoted by ji and j2, one could try to reabsorb the scheme dependence 
by an evanescent redefinition ji ^ j[ = {1 + ae + . . .)ji + e ■ (other currents). I.e. j[ yields 
with prescription 1 the same finite parts of physical amplitudes as j2 renormalized with 
prescription 2. 
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5. Effective Lagrangian for Flavour 
Changing Processes 

In this chapter wc will set up the operator basis needed for the calculation of the K'' — K*^- 
mixing described in chapter 6. We will also explain the general structure of the effective 
lagrangian and of the calculations which are necessary to obtain its RG improved coef- 
ficients. The basic framework is very similar for the general case of particle-antiparticle 
mixing. We will sometimes comment on differences between the considered |AS|=2 tran- 
sitions and |AB| = 2 or |AC| = 2 amplitudes. 

First the Standard Model diagrams contributing to the K*^ — K° -mixing will be analyzed. 
Since | AS| = 2 transitions involving light quarks contain | AS| = 1 processes as substructures, 
sect. 5.3 will discuss them in detail. Here we will notice that the definition of the correct 
I AS I = 1 operator basis involved in |AS| =2 amplitudes is much more complicated than the 
one needed to describe the four quark amplitudes appearing in |AS| = 1 weak decays. The 
analysis of the operator basis requires the careful application of the theorems of sect. 3.10 
on the use of the field equations of motion and on BRS-exact operators. 

5.1. |AS| = 2 transition in the Standard Model 

The I AS| =2 transition is a flavour-changing neutral current (FCNC) process and therefore 
forbidden at tree level. The lowest order contribution to it is depicted in flg. 5.1. 

Figure 5.1: The lowest order box diagrams mediating a |AS| = 2 transition. The zigzag lines 
stand for W-bosons or ficticious Higgs particles. The diagrams rotated by 90° must also 
be considered. 



5.1.1. Notations and Conventions 

Before writing down the result for the diagram of fig. 5.1, we set up the notations used in 
the following. 

The different contributions from the internal quarks involve different CKM factors A-,- = 
VjdV*g. The GIM mechanism + Ac + A„ = allows for the elimination of A„. Hence we 
can split up the Standard Model Green's function as 



G = A^G^ + A^G* + 2AeAtG'^*, 



(5.1) 
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which is understood to be truncated, connected and Fourier-transformed into momentum 
space. We will sometimes use the abbrevation Xi — ml/M^ for the squared ratio of some 
quark mass and the W mass. For the W propagator the 't Hooft-Feynman gauge will be 
used, while we keep the QCD gauge parameter ^ arbitrary. As only open fermion lines 
appear in the calculation, we can safely use a naive anticommuting 75 (NDR scheme) with 
L = 1 — 75 and R — 1 +75. is the number of colours, and 1 and 1 denote colour singlet 
and antisinglet-*^, i.e. (L (8) i?) • 1 stands for Sj(l — 75)^^ • + 75)(ii with j and k being 
colour indices. The SU{N) Casimir factor involved will be Cp — {N'^ — 1)/{2N). 
The result of fig. 5.1 is proportional to the tree-level matrix element of 

Qs2 ^ {sd)v-A{sd)v-At = (s7^(l-75)d)(s7''(l-75)rf)-l (5.2) 

shown in fig. 5.3. In four dimensions the Fierz transformation maps Qs2 onto an operator 
with the same Dirac structure, but 1 replaced by 1. Hence if we substituted 1 by 1 or 
by (1 + i)/2 in the definition (5.2), we would redefine Qs2 by an evanescent operator 
proportional to (s(i)v-yi v-a{^ — !)• If the evanescent operator Ei[Qs2] is defined as 
in (4.11), the Fierz symmetry is maintained on the loop level and the replacement 1 ^ U 
does not change the two-loop anomalous dimension. 

The expansion of matrix elements in terms of a has already been defined in (3.51). Anal- 
ogously we will expand the G's in (5.1), 

& = + + 0{a^), (5.3) 

and the Wilson coefficients: 

C, = cr + ^C« + 0(a^). (5.4) 

The G^' 's involve infrared (mass) singularities, which will be regularized by small quark 

masses. 

Finally the weak coupling constant gw will be traded for the Fermi constant Gp according 
to (3.21). 

5.1.2. Zeroth Order Amplitude 

In the leading order of ^ught/'^^hcavyj where mhcavy stands for nit or Miy and miight denotes 
any other massive parameter, one can neglect the external momenta in (5.1). 
One finds for the mixed top-charm contribution in (5.1): 

= ^M^S{x,,xt){Qs2)^"l (5.5) 

Here the Inami-Lim function [14] S{xj,Xk) is defined as 

S{xj,Xk) = S{xj,xk)-Sixj,0)-S{xk,0) + S{0,0), (5.6) 

^This is clearly a misnomer from a group theoretical point of view. 
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where the result of the box diagram with internal quarks j and k is denoted by S{xj,Xk) 
and the u-quark mass is set to zero. For j — c,t one finds: 



167r2 



M^S{xj){Qs2) 



(0) 



(5.7) 



with S{xj) = S{xj, Xj). Here one realizes that the efi^ect of the GIM mechanism is not only 
to forbid FCNC's at tree level, but also to cancel the constant terms in the <S's and to 
nullify K° — K° -mixing in the case of degenerate quark masses. 

The reminder of the thesis will deal with RG improved short distance QCD corrections in 
(5.1), which are parametrized by 

771 on the RHS of (5.7) for 

772 on the RHS of (5.7) for G* and 
?73 on the RHS of (5.5) for G^. 

I.e. we want to derive the low energy hamiltonian given in (1.1). 

Let us first look at the three contributions (5.5) and (5.7) to (5.1) with emphasis on large 
logarithms: 



S{xt) = xt 
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(5.8) 



clearly involves no large logarithm because of Inxj 
a{Mw)l{^T^) this is of the order 10"^. 

S{xc) = Xc + 0{xl) 



1.5. After multiplying with 



(5.9) 



Here we have only kept terms which are larger than those of order {msfric) / neglected 
by setting the external momenta to zero. As we will see in the following section one would 
naturally expect a large logarithm In Xc multiplying Xc here. Its absence is due to the GIM 
mechanism. 



S{Xc, Xt) 



- 8xt + 4 , 3 Xt 



4(1 -xt)^ 



In + - 



4xt-l 



+ 0{x^ In Xc 



(5.10) 



Here we encounter a large logarithm | In^d ~ 8. According to our discussion in sect. 3.7.2 
we can sum In Xc(q; In Xc)", n = 0, 1, 2, . . . with the help of factorization and RG techniques. 
S{xt) is much larger than S{x^ and S{xc,Xt), but it is CKM suppressed in (5.1). Yet this 
is not so in |AB| = 2 transitions, where S{xt) and its radiative corrections comprised in 772 
are clearly dominating. 
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Figure 5.2: The classes of diagrams constituing the 0(a)-contribution to G in (5.1); the 
remaining diagrams are obtained by left-right and up-down reflections. The curly lines 
denote gluons. Also QCD counterterm diagrams have to be included. Diagram Fg equals 
for zero external momenta. 



Figure 5.3: The diagram for the matrix element of Qs2 in the effective three-quark theory 
to order q;°. The cross denotes the insertion of the effective A^* = 2 operator. 
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5.1.3. 0{a) Radiative Corrections 



Since we are worried about the absence of the logarithm in S{xc), let us look at the one-loop 
radiative corrections to G'^ depicted in fig. 5.2. 

These diagrams have been calculated for arbitrary internal quark masses in [15]: 



In (5.11) new operators have emerged: 
f 
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(5.11) 



(5.12) 



We have written U and T in a manifestly Fierz self-conjugate way. The functions in (5.11) 
are: 
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(5.13) 



Let us discuss (5.11) in more detail: First (5.11) is obviously unphysical, because it is gauge 
dependent. This is an artefact of the use of small quark masses to regularize the infrared 
singularities while at the same time using zero on-shell quarks for the external states. For 
the same reason we encounter the new operators T and U. When factorizing (5.11) to 
extract the Wilson coefficients C^, these unphysical parts will match their counterparts in 
the effective theory with the zeroth order coefficient Cf \ Hence they do not contribute 
to the Cfc's. Yet we can anticipate this by noticing that the terms involving T, U and ^ 
are independent of Mw-, which sets the scale for the short distance physics. An alternative 
method to treat the infrared singularities is to keep the quark masses exactly equal to zero 
and to rely on dimensional regularization to cope with both IR and UV singularities. This 
gives a different result, as the limit m — > is non-uniform owing to the fact that m^ /e 
has a zero mass limit, while the expanded version 1/s + Inm diverges in this limit. If we 
had used such an on-shell IR regularization, the additional operators T and U would be 
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absent. Since further the wave function renormahzation constant equaled zero then, 
(5.11) would turn out to be gauge independent, too. The use of small quark masses for 
regularization has the advantage that the vanishing of the otherwise invisible IR terms 
from the C^'s provides a check of the calculation. Yet one may argue that this method 
requires a more complicated calculation. In fact this is not so as we will see in chapter 6. 
Having in mind the RG evolution down to scales //c of order rric we have arranged the 
logarithms in h in (5.13) such that one can easily distinguish those which are small for 
\x ~ vric from the large logarithm log(m^/M^). For this we have abstinated from using 
Xc, but have turned to m^. We will do this again without warning. Due to the magic 
of factorization the coefficient of In (m^/M^) is a linear combination of the anomalous 
dimensions of composite operators of some effective field theory. 

Let us finally remark on the //-dependent terms in h: The colour octet term proportional to 
(A^ — 1)/(2A^) is /i-independent. In the colour singlet part proportional to Cp one identifies 
7^-* = 6Ci? as the coefficient of In /i^, which stems from the running charm quark mass in 
(5.7). Likewise 2CF^\\i\j? originates from the wave function renormahzation. 

5.2. General Structure of the RG Improvement 

After integrating out the heavy top quark and the W-boson one is left with an effective 
field theory containing a lagrangian of the form (3.28). Now we extract the Fermi constant 
out of the Wilson coefficients to render them dimensionless: 

^|AS|=2 _ ^^T/ n ni^) ^^P\ri H^)n /'c; 

Here Vckm,f and Vckm,f stand for products of CKM elements and F is a flavour index. 
The light internal quarks in fig. 5.1 and fig. 5.2 require to match the corresponding parts 
of G to matrix elements with two insertion of |AS| = 1 operators Qfe. Fig. 5.4 and fig. 5.5 
show such diagrams for the case of two insertions of current-current operators. The Qi% in 

Figure 5.4: Diagram Dq in the effective five-quark theory. The cross denotes the insertion 
of a A^* = 1 current-current operator. 

Figure 5.5: The classes of diagrams giving the Oia) contribution to {Oii) and (T^jj) for 
J = 1,2. The other members of a given class are obtained by left-right and up-down 
reflections. Also QCD counterterms have to be included. Diagram Dg = for zero external 
momenta. 

(5.14) correspond to |AS| =2 operators. In general their matrix elements also participate in 
the matching at the initial scale \iw We will be interested in the contribution of the lowest 
dimension operators: The Qfe's in (5.14) have dimension six and the QiS are dimension- 
eight operators. The remaining sections of this chapter will concern the correct reduced 
operator basis needed in (5.14). 
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Figure 5.6: The classes of diagrams in the effective three-quark theory contributing to {Qs2) 
to order a. The other members of a given class are obtained by left-right and up-down 
reflections. QCD counterterms have to be included. 

Next one has to carry out the program of sect. 3.10, 3.7 and 3.8 consisting of the following 
steps: 

i) Find a reduced operator basis for (5.14), from which renormalized operators of the 
types described in sect. 3.10 and properly renormalized evanescent operators have 
been removed. 

ii) Match the Standard Model Green's function G in (5.1) to the matrix elements de- 
rived from (5.14) as shown in sect. 3.7. This is to be done at a scale satisfying 
Mw < ^ rrit, but being arbitrary within this range. 

iii) Follow sect. 3.8 to perform the RG evolution down to /Xc ~ tUc thereby summing 
ln(|Lic/ ^w) to all orders. When passing the b-quark threshold remove this degree of 
freedom from the calculation. This simply modifies the running coupling and triv- 
ially changes the penguin operators to be introduced in sect. 5.3.2. Strictly speaking 
one matches (5.14) to another lagrangian with the b-quark integrated out, but the 
structure of both lagrangians is essentially that of (5.14). 

Hence in (5.14) describes the physics between the scales jic and jiw- At the scale 

11 c one has to match the matrix elements derived from (5.14) to those of a new effective 
lagrangian with the c-quark integrated out and to repeat steps i) to iii). This 

amounts to the matching of diagrams of the type of fig. 5.4, fig. 5.5 to those with a single 
local four-quark operator as depicted in fig. 5.3 and fig. 5.6. In fact the reduced C^^2 
contains only a single physical dimension-six operator, which we have already met in (5.2): 

-^ifffJ ^ = {psl^l + ^sl^l + C^S2^c^t) QS2 (5.15) 

In the matching we equate dimension-eight Green's functions derived from C^^^^ with di- 
mension-six Green's functions from ^C^g^J"^. Hence by power counting the C'52's in (5.15) 

pick up two powers of mc(/ic)- The subsequent RG improvement of the (552 from 
down to some scale // will sum In /ic/ The final scale /i must be large enough to trust 
into perturbation theory, this is roughly of the order of 1 GeV. The matrix element of 
Qs2{fJ-) must be evaluated between neutral K-meson states. The calculation of such QCD 
binding effects can be done in the framework of lattice gauge theory, 1/A^-expansion or 
with sum-rule techniques. 

The observation that we need two powers of rric in the matching greatly simplifies the 
finding of the correct set of dimension-eight operators Qi needed in C^^^l~'^. The only 
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candidate is 

2 

777 

Qi = -^s^i^Ld-s^Ld, (5.16) 

which we aheady know from an example in sect. 3.8.2. To understand this recall that 
any other dimension-eight operator lacks the factor of and involves derivatives and/or 
gluon fields instead^. The matrix elements between s-quarks and d-quarks evaluated for 
the matching at /ic cannot produce a factor of m^: The c-quark can only enter matrix 
elements with single insertions of such operators through loops in gluon propagators, but 
the result is always proportional to k'^g^" — k^^k" rather than m^, where /c^ is the momen- 
tum flowing through the propagator. Hence the QCD gauge structure keeping the gluon 
massless provides us with unexpected help here. For the same reason other dimension- 
eight operators cannot mix into during the RG evolution from /iw to fj,c- In principle 
operators with 'p~^^m.l/ g^ multiplying other Dirac structures can also appear, but we will 
see that this is not the case due to the flavour structure of the problem. Let us drop some 
words on the factor /Z"^*^ in (5.16): Bare quantities must not depend on /j,. Therefore the 
bare operator reads 

777^ 

fl'bare 

With (2.34) and (3.2) one is therefore lead to (5.16). It is also clear from fig. 5.4 that 
the dimension of is larger than the dimension of the QkS by 2 — 2e, because the loop 
integration gives a power of D — 2. We have refrained from defining powers of yU^^ into the 
Cfc's, because this would lead to an unnecessary complication in the RG equations (3.44). 
The origin of these simplifications of jO,lf[J~ in (5.14) is the fact that the internal c-quark 
is heavier than the external s-quarks and d-quarks. The situation is therefore completely 
different for the case of the c-quark contributions to |AB|=2 transitions or the effect of 
internal s-quarks in |AC|=2 amplitudes. Here the external quarks are heavier and the 
dimension-eight operators containing derivatives dominate over the effect of /I^^^m^/gf^ ■ 
7^L (g) ■j^L in |AB| = 2 amplitudes resp. Ji~'^^m1/g'^ ■ 7^L ® ■j^^L in |AG| = 2 amplitudes. 
Fortunately the former are phenomenologically uninteresting, as they are suppressed with 
mil ml, and in the case of |AG| = 2 transitions there are other badly controllable long- 
distance effects in the game which most hkely make short distance calculations of D° — D° - 
mixing useless. 

5.3. |AS| = 1 -Operator Basis 

Let us now find a basis for the |AS| = 1 operators Qk in (5.14), which closes under renor- 
malization. 

^Operators containing only one power of nic can be removed by the use of the EOM for the fermion 
field according to the procedure described in sect. 3.10 giving msirnc/g^ ■ 'j/j.L (g) j^^L. 
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5.3.1. Current-Current Operators 

Consider first the matching procedure for the LO diagram fig. 5.1 with only light quarks in 
the intermediate state. Contracting the W boson hnes yields the diagram of fig. 5.4 with 
two insertions of the current-current operator of fig. 5.7: 

= s-ff,Lk-h^Ld-l, (5.17) 

where k and I are u or c. 



Figure 5.7: |AS| = 1 current-current operator 



For k 7^ I only mixes with 

= s'y^Lk-hi.Ld-l. (5.18) 

The corresponding diagrams are shown in fig. 4.1 on p. 43. The 2x2 mixing matrix in 
LO and NLO has been derived in [15]. If the corresponding evanescents are defined as in 
(4.11), the NLO anomalous dimension matrix commutes with the LO one. These matrices 
are diagonal in the basis containing the operators g^' and Q'^ defined by 

Q'i = l(g^'±gf). 

Their anomalous dimensions in the NDR scheme with (4.11) are 

7f = ±6^, = ^(-21±'-l^'lN±'-f). (5.19) 



5.3.2. Penguin Operators 

We will focus first on the LO operator mixing. For k — I the current-current operators mix 
with penguin operators. The simplest divergent one-particle irreducible penguin diagrams, 
into which Q\'^2 ^e inserted, are depicted in fig. 5.8. These diagrams require the gluon- 



Figure 5.8: The mixing of a |AS| = 1 current-current operator into gluon-foot penguin 
operators. Diagrams with crossed gluons compared to G2 and G3 must be included, too. 

foot penguin operators Qgi, Qg2 and Qg^ of fig. 5.9 as counterterms. We will encounter the 



Figure 5.9: The gluon-foot penguin operators 



basic penguin diagram Gi later as a subdiagram, so that we comment on it here: Although 
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it is a dimension-two diagram, its result contains no term proportional to the square of 
the internal mass. One has instead 



Gi oc s^^.LT^d ■ [p'^p'' - p'g"'') , (5.20) 

where u and a are the open Lorentz and colour indices of the gluon line and p is the 
momentum flowing into the gluon vertex. (5.20) holds to all orders in D — 4. Especially 
the inserted operators Qi'/2 do not mix into an operator of the form rTn?s^"'LT°'d. Likewise 
the dimension-one diagram G2 does not give a positive power of the internal mass. 
In order to find a closed operator basis one now has to investigate into which operators 
Gi, G2 and G3 mix. One readily finds that e.g. Gi mixes via box diagrams into new four- 
fermion operators. A box diagram with two fermion lines containing one Qgi vertex and 
three QCD gluon-fermion vertices is divergent, because Qgi involves two powers of the 
loop momentum due to (5.20). These new operators are fermion-foot penguin operators: 

q=d,u,s,c,b 
q=d.,u,s,c.,b 

Qe = s-f^Ld- Y Q^'^Rq-^ (5.21) 

q=d,u,s,c,b 

We will be more interested in the mixing of Q1/2 into them, which is shown in fig. 5.10. 
At this point we realize that our simple minded approach to the operator mixing leads 



Figure 5.10: The mixing of a |AS| = 1 current-current operator into fermion-foot penguin 
operators depicted on the right. 

to problems, as we could renormalize the divergence in fig. 5.10 either by a counterterm 
proportional to Qgi or by one involving fermion foot operators. Further we obtain mixing 
into a third type of operator as depicted in fig. 5.11: 

Qfp = s^^LT-d-{d^rJ,)rjJabc. (5.22) 



Figure 5.11: The mixing of a |AS| = 1 current-current operator into ghost-foot penguin 
operators depicted on the right. 

Fortunately we can put order into the penguin zoo with the theorems of sect. 3.10. For this 
we first construct the dimension-six operators (with total ghost number zero) from class 
B: The only possibility is given by Qbrs from (3.58). It is related to Qfp via 



<5brs — Qfp + Qgi, 
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where Qgf is the term on the LHS of (3.57) stemming from the gauge fixing part of the 
QCD lagrangian: 

= ^^s^^LT^d-d^'d^Al. (5.23) 

Next we arrange the operators found so far into combinations belonging to class S. Here 
only the EOM for the gluon field can give dimension-six operators related to the penguin 
zoo. The relevant operator has been introduced in (3.57). In terms of the operators defined 
in (5.21), (5.22) and (5.23) it reads: 

Qeom = ys^^.Ln-D,,F^'"'-^{Q, + Qe) + ^{Qs + Q5)-Qsf + QFp{5.24) 

Now all other operators involved in the mixing must combine to gauge invariant operators, 
because they belong to class V. We have already listed Qf, . . .Qq, the only other gauge 
invariant dimension-six operator is 

yS^^LT-d-D^aF'"''' = Qgl + QgS + QgS: (5.25) 

which we can drop from the basis, because it is a linear combination of Qeom, Qbrs and 
the other physical operators. 

Hence the unreduced |AS| = 1 operator basis consists of Q^^ ■ ■ - Qe, Qeom, Qbrs and of 
course of evanescent operators. This structure implies that the mixing of any operator 
into Qgi, Qg2 and Qg^ must be the same. Indeed, the divergent parts of Gi, G2 and G3 all 
involve the same coefficient, as required by (5.25). The theorems of sect. 3.10 also explain 
the observation of (5.20), because there is no dimension-six operator in V®S®B involving 

Due to sect. 3.10.2 and chapter 4 the reduced lagrangian for |AS| = 1 transitions is simply 
obtained by dropping the renormalized Qeom, Qbrs and the renormalized evanescents from 
the effective lagrangian. 

The on-shell vanishing of (Qeom) can be easily understood in terms of diagrams [54]. The 
Feynman rule for the one-gluon piece of Qeom depicted in Qgi of fig. 5.9 involves an inverse 
gluon propagator: 



9 



(5.26) 



Hence in any diagram in which the gluon line of Gi ends up in an internal vertex the gluon 
line is effectively shrunk to a point, because (5.26) is contracted with a gluon propagator. 
The resulting diagram is identical with a diagram containing one of the other operators in 
(5.24) and these contributions cancel. If the gluon of Gi is an external fine, the diagram 
vanishes after LSZ reduction. Care is necessary, however, if renormalization constants are 
calculated off-shell, because then Qeom and Qbrs give non-zero contributions, which must 
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be properly taken into account. For example in fig. 5.10 one finds off-shell terms involving 
s^i^Ld • (p^p^/p^) /s which require a counterterm proportional to Qbrs- 
Now for the case of the |AS| = 1 substructure of |AS| = 2 diagrams the situation is more 
complicated, because here the gluon of Gi can end up in another Q^i-vertex. Then we have 
one propagator hitting two inverse propagators and the result is non-vanishing. But from 
sect. 3.10.3 we know that this contact term is identical to the matrix element of a local 
dimension-eight |AS| =2 operator Qi. Hence we can still restrict the renormalized QkS in 
(5.14) to the set {Qf,...Q6}- 



5.4. Operator Basis for |AS| = 2 transitions 



I I 2 

Let us now assign the correct CKM factors to the effective lagrangians in (5.14) 

and >Cgg2~^ ill (5.15). 



5.4.1. Above Hc 

The reduced |AS| = 1 lagrangian can be taken from [44]: 



£|AS|= 



G 



F 



V2 



E E v:,ViaC,Q^'-Xtj:C,Q, 

i=l k,l=u,c j=3 



(5.27) 



Here the coefficient At originates from the fact that the matching of the penguin operators 
and the mixing of any four-quark operator into penguin operators involves the diagram of 
fig. 5.10 and radiative corrections to it all being proportional to A„ + Ac 



-A. 



With the findings of sect. 5.2 the effective lagrangian C^^^i of (5.14) is found as 



,|AS|=2 
'eflf,! 



£|AS|=1 _ 



G' 



9' 



C7Q7 



G F 



V2 



E E v:,viac,Q^'-\tj:cjQj 

1=1 k,l=u,c 



j=3 



G\ 



9' 



We may ask where the contact terms from the reduction of the |AS| =2 lagrangian have 
gone. Since diagrams with two insertions in which the two | AS| = 1 vertices are connected 
by a gluon line or by a ghost-loop cannot produce a factor of m^, these terms are absorbed 
into other dimension-eight operators not containing m^, which have been found to be 
irrelevant in sect. 5.2. 

Since there are no more derivative couplings involved, we can pass to the more common 



hamiltonian formalism by simply flipping the sign: H = — £ 



AS|=2 

cfr,i 



We will do further and 



deflne the |AS| =2 hamiltonian such that it only contains terms of order Gjp. For this we 



equate the first order term of the Gell-Mann-Low series exp[— i Jd'^xH['^^^-\x)] with the 
corresponding series involving £gg^^'~^: 



i J d^xH[^^^-\x))\^^s\=2 = (e 



\\AS\=2 + 0{Gp) 
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This notation is a bit sloppy, but useful. Let us write 



|AS|=2 



ct 



One then obtains 

H%x) 

H\x) 



H'Hx) = 



>^cXt^ E J2C^CJn,,ix) + K\t^CrQr{x) 
^ i=+-j=i ^ 



(5.28) 



with 



-gr(^)Qr(i/) - Qr(^)Qr(2/) - QT{^)Qr{y) 



Tlij{x) = < 



/A T -Qn^)Q'^^{y) - QT{x)QY{y) - Qrix)Qny) 



-Qnx)QY{y) + 2Qr{x)QY{y) 



iforj=l,2, 



\jd'yT [{Qnx) - QY'^{x))Q,{y) 



(5.29) 



+g,(x) (gr(t/) - gr(y))] 

The structure of (5.28) requires some explanation: 



for j=3,. . . 6. 



The top top contribution H* 

First there is only a |AS| =2 operator involved in H*. Wc have only listed the terms relevant 
for the leading term in the expansion in l/rriheavy In the case of if* this contribution comes 
from the internal top-quarks in fig. 5.1, which are matched at iiw to the local operator 
g52- Prom (5.7) and (5.8) we know that Cg\ is proportional to Xj, while the double penguin 
contributions can at most give a factor of Xb- Also the dimension-eight operators stemming 
from the contact terms in (3.62) and (3.63) can only give contributions proportional to 
m^gijj/M^, because there is no more top-quark in the effective theory. The calculation 
of 772 therefore only requires the consideration of matrix elements with a single |AS|=2 
operator. Yet the NLO matching has to be done from the finite parts of the two-loop 
diagrams in fig. 5.2. This calculation has been carried out by Buras, Jamin and Weisz [15]. 
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The charm— charm contribution H*^ 

Second there is no local |AS| = 2 operator involved in This feature is related to the 
absence of a large logarithm in (5.9). The GIM combination in Oij in (5.29) leads to 
the cancellation of the divergences and thereby of the ln(///Mv^)'s in the diagrams (5.4). 
The same holds true for the diagrams of (5.5) after the inclusion of counterterms for the 
|AS| = 1 operators and in fact to all orders in perturbation theory. This has been first 
observed by Witten [39]. Hence if^ belongs to case 11 of sect. 3.7.2. Yet we may conclude 
from sect. 3.7.2 only that there is no mixing into a local operator, but there could be a 
local operator involved in the matching at /iw, which could evolve down to /ic unaffected 
by the presence of the bifocal structures Oij. But in the case at hand the matching at fi\y 
can be made by the replacement 1/(M^ — k"^) l/M'^r in the diagrams of fig. 5.1 and 
fig. 5.2. This is allowed, if the renormalized Standard Model diagram does not become 
UV-divergent by shrinking the W line to obtain the effective diagrams of fig. 5.4 and 
fig. 5.5, because then the contribution of the UV momenta is suppressed. In H'^ the GIM 
mechanism cancels the UV divergences in the effective theory, so that there is no room for 
a non-zero coefficient of Q-j. This fact becomes very clear if one splits the loop integral 
connecting the two | AS| = 1 vertices into a part with momenta smaller than rric and an UV 
part comprising loop momenta above rric- For the first part we can use the factorization 
of the |AS| = 1 substructure yielding (5.29) up to terms of order m^/M^, and the UV 
momenta cannot give non-negative powers of to spoil this factorization, because the 
integral from rric to oo vanishes in the limit My/ oo after forming the GIM combination 
(c, c) — (c, u) — {u, c) + {u, u). Hence we only need to know the | AS| = 1 coefficients C+ and 
C_ , which have been determined in [43] . No extra | AS | = 2 matching in necessary. From 
sect. 3.7.2, however, we know that we need the finite parts of the diagrams of fig. 5.5 for 
the NLO matching at the scale /Xc, so that we can check the absence of a local operator in 
H'^ in (5.29) by comparing the matrix element of H'^{iiw) in (5.28) with (5.7) and (5.11): 

iG'^{f,w) + 0{xl) = {H^fiw)) = Xl^Ci{i,w)CjM{OijM). (5.30) 

On the other hand the operator mixing is simple and does not require new calculations. 
We have put the current-current operators in 0y into the diagonal basis to simplify the 
RG evolution. The |AS| = 1 Wilson coefficients C+ and C_ in (5.30) simply evolve down 
to fjL — He according to [44] . 

The charm— top contribution H*^* 

H'^^ involves the familiar large logarithm in (5.10), so that it belongs to case I of sect. 3.7.2. 
Hence both the LO and the NLO matching can be done from the one-loop graphs fig. 5.1 
and fig. 5.4. Again we have chosen the diagonal basis in TZij for the current-current vertex. 
For the other vertex we have chosen the operator basis , . . .Qeto use the results of [44]. 
Now the LO matching at // = fj,w is performed from the large logarithm In Xc in G"^' . In 
NLO we also equate the nonlogarithmic part of C^' to the diagrams of the effective field 
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theory to obtain the initial value C7(/xvk)- Since the initial value of the | AS| = 1 coefficients 
Cfe are of order a for A; = 1, 3, 4, 5, 6 and Ck — 1 + 0{a) for k — 2,+, — , the matching at 
the scale /i — /iw reads 



0(x^ InXr 



(0) 



(5.31) 



Here we have accounted for the fact that Cr starts in order a due to the inverse power in 
the definition in giving (g7(/x))(°) = mK^i) / {4na{fi)) ■ {Qs2{f^))^^^ ■ Hence in LO there 
is nothing to adjust in (5.31) and the coefficients of the large logarithms must be the same 
on both sides. This is indeed the case. Further for the NLO matching only the diagrams in 
fig. 5.1 and the current-current diagrams of fig. 5.4 are needed. For later convenience we 
up: 

H'^i^i) = //^*'+(/i) + i/^*'-(/i) 



split i/'^* 



G 



^ AcAj 



ffct,+ H^^- evolve independently under the RG fiow. Also C7{n) is split into C^{ij) + 
C^{lJ,), both of which receive different admixtures during the RG evolution. It is arbitrary 
which part of C7(ijiw) is put into C^(iJ,w), as long as their sum fulfils the matching condition 
(5.31). We will choose C^'^^\fjLw) = 1/2 • Ci'^\tjLw) 

The RG evolution is much more involved than in the case of H"^. With the well-known 
initial values of the |AS| — 1 coefficients [43,44] and with CyI/iw) obtained from (5.31) one 
knows 



C+{lJ'w)C2{(J'w) 

C+{lJ^w)Cz{nw) 
C+{(^w)Ca{ij.w) 
C+{lJ^w)Cr>{^-w) 



C_{lxw)C2{(J'w) 
C-{iiw)C^{^w) 

C-{lXw)CA{iJ'w) 
^7 {l^w) 



(5.32) 



to order a. These are the initial conditions for the NLO RG evolution. We will give the 
relevant expressions for Ck{fJ>w) and C^{fj,w) in chapter 6, where the actual calculation is 
described. 

In NLO all coefficients in have a non-zero initial value. The structure of the RG 
evolution has been described in sect. 3.8.2. In (3.46) the summation indices m,m',v and 
v' assume the values + and — and t,t',n and n' run from 1 to 6. The two terms of the 
summation over n correspond to and H'*'" and the two coefficient vectors in (5.32). 
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The first term in (3.46) involving the initial value C7{fj.w) can obviously be arbitrarily 
distributed over and H'^^~. 

The structure of the RG evolution is probably more transparent from (3.44): The anoma- 
lous dimension tensor '~fkn,i describing the mixing of |AS| = 1 coefficients into the |AS|=2 
coefficient 6*7 contains 12 elements ■jn+,7,'ln-,7, ^ = 1, • • - 6, which we summarize into the 
row vectors 

7+,7 = (7l+,7 , 72+,7 , 73+,7 , 74+,7 , 75+,7 , 76+,7) 

7-,7 = (7l-,7 , 72-,7 , 73-,7 , 74-,7 , 75-,7 , 76-,7) (5.33) 

Now (3.44) reads in the case at hand: 

I^^C^ = 7C7 + 7fc+,7 Ck C+ + 7fc_,7 Ck C- (5.34) 

which is split into 

l^-^C+ = 7(57+ + 7fc+,7CfcC+ 

/x-^(5f = 7Cf + 7fe_,7 CfcC_. 

Here 7 is the anomalous dimension of Q^. This yields the RG equation for C^"^^ defined in 

/.AcW(;.) = ( ^7^' McW(/^), (5.35) 

where 7 is the 6x6 | AS| = 1 anomalous dimension matrix and 7+ and 7_ are the anomalous 
dimensions of (5+ and (5.35) leads to a block-triangular evolution matrix of the form 
(3.48) and one again recognizes that Cii^ij) = Ct{lj) + C^{lJ) does not depend on how the 
initial value C'!{ixw) is distributed over C^{^w) and Cj{^w)- 

The solution of (5.35) involves two 7x7 evolution matrices. They, however, encode redun- 
dant information, because each of them contains the full 6x6 |AS| ~1 evolution matrix. 
Now one can do better: Define the vector 



D = 



(5.36) 



where C = {Ci,...Cq) is the |AS| = 1 coefficient vector. From (5.34) one finds the RG 
equation for the seventh and eighth component of D: 

I^^Dr = (7-7+)L>7 + 7^,7-^, 



f^-^Ds = (7 - 7-)^8 + 7-,7 • ^■ 



(5.37) 
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Hence D satisfies the following RG equation: 



/ 7 



T ~ 

r+,7 7 



7+ 






1-1- J 



D. 



(5.38) 



— * 

After obtaining the 8x8 RG evolution matrix in the standard way one gets D{ij) which 
has to be multiphed with C±{ijl) to obtain C^^\ijl) and 

H'^'Hi^) = -^AA E (7^l,n(/x),...,7^6,„(/x),g7(/^))•(^("H)") 



n=+,- 



Hence we have switched from a 7 x 7' representation to an 8 x 1 x 1' representation of the 
renormalization group. We have programmed both the solution of (5.35) and (5.38) and 
checked that we obtain the same result. Of course, the program for (5.38) is much faster. 
Further the upper left 6x6 block-submatrix of the 8x8 evolution matrix is the |AS| = 1 
evolution matrix of [44], which provides an additional check of the RG evolution. 
Now the calculation of the elements 'jk±,7 oi the LO and NLO anomalous dimension tensor 
has required the calculation of the divergent parts of the diagrams of fig. 5.4 and fig. 5.5, 
but now with insertions of TZij for j = 1, 2, rather than Oij and of the diagrams of fig. 5.12 
and fig. 5.13 corresponding to matrix elements of TZij for j' = 3, . . .6. Further there are 



Figure 5.12: Diagram yielding the LO anomalous dimension tensor 7 for k = 3. . .6. 
The white cross denotes the insertion of a AS* = 1 current-current operator. The dark 
cross stands for one of the | AS| = 1 penguin operators Q3, . . . Qq. 

diagrams in which the s-quark foot of the penguin operator is connected with the current- 
current vertex and likewise diagrams with the d-quark foot involved. An example is shown 
in fig. 5.14. For insertions of Q5 and Qq these diagrams vanish by chirality. For insertions 
of Qs and Q4, the diagrams are nonvanishing, but their sum cannot give a factor of in 
the divergent part. This follows from a subtle analysis exploiting (5.20) and the theorems 
on the mixing into gauge invariant operators. 



Figure 5.13: The classes of diagrams contributing to {TZij) for j = 3, . . .6, to 0{a). The 
other members of a given class are obtained by left-right and up-down refiections. Also 
one-loop counterterms have to be included. 



Figure 5.14: An example of a diagram involving the s-quark foot of the penguin operators. 
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5.4.2. Below /^c 

At the scale /^c ~ wc remove the c-quark from the theory. Hence we are left with the 
single local operator Qs2- We have hsted the low energy hamiltonian already in (1.1) and 
(5.15): 



H'^il') = ^K\cf^{lAQs2{lA = -^M^2AA?73^(^c,^t)fe(/^)Q52(/^).(5.40) 

The only effect of the charm threshold for Cg\ is the transition to the three-quark running 
a. 

The coefficient Cgiil^c) (5.40) is determined by the matching of the matrix element of 
//^(/Xc) in (5.28) to those of //^(//c) in (5.40): 

rni^c)) = ^A,2 E a(A^c)Q(/^c)(0..(/^c)) = 3^A^42n/^c)(g52(A^c))(5.41) 

This yields C''gl{^c) in terms of the RG improved |AS| = 1 Wilson coefficients Ck{^c)- Here 
we verify the findings of sect. 3.7.2 that in the NLO the finite parts of fig. 5.6 and of the 
two-loop diagrams in fig. 5.5 are necessary. 

The determination of Cgf{iJ,c) requires the same procedure with i?"^* of (5.28) and (5.40). 

Consider ffist the LO matching: The coefficient Crdj-c) is nonzero due to admixtures from 
C2^\f^w) picked up during the RG evolution from to //c- Hence the LO matching is 
done from the inverse power of q;(/Xc): 

Cg'(^J = ^J^JL.CM. (5.42) 

In the NLO matching all coefficients participate. We define 

C^M) = ^2r,,,52(/.)(g52(/.))(°). (5.43) 

Here the definition of ?"jj,52 differs from the one of rnm,i given in (4.46) by a factor of 2. In 
NLO one therefore finds: 



1 47i ^ 



0^/„ nC'7(/^c) + Yl T.^'ijMl^c)Ci{lIc)Cj{lIc) 



(5.44) 



On the side corresponding to the three-quark theory one expects the diagrams of fig. 5.3 
and fig. 5.6 corresponding to (^52)^°^ and ((552)*^^^ to be necessary. They, however, cancel 
with the terms accompanying (Q^) on the side of the four-quark theory. 



5.5 Why Working Beyond Leading Logarithms? 
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Now one is left with the RG evolution from Hc down to some scale n at which the hadronic 
matrix element (K°|Q52(a*)|K°) is obtained. The RG evolution of Qs2 is well-known [15,43], 
its anomalous dimension is 7+ given in (5.19). In terms of (3.41) and (3.43) the NLO 
solution for CglilAi 3 — ^1 ^1 reads: 



^(3) a(/xe)^-a(/x) \ 



«(//) 

Finally we can express the NLO riiS in (5.40) in terms of the coefficients 

^ ni<')( \ r ( \^<if ( ^ 7(3) "(/^c) 



The //-dependence in (5.45) is absorbed into h{jj), which equals 

m = [.(/.)]-<' (l + jf^^) (5.47) 

in the NLO. The index (3) denotes that one has to set / = 3 in the corresponding quantities. 
Prom (5.46) it is evident that the Ty^'s depend on the definition of the masses. We will 
adopt the convention of [15, 17] that the running masses in (5.46) are evaluated at the 
scale at which they arc integrated out, i.e. trie = mc{fic),^t = 'm-tifJ'w) ■ Whenever the 
TjiS are defined such that they multiply S{xc{mc)), S{xt{mt)) and S{xc{mc),Xt{mt)) in the 
effective hamiltonian in (5.40) we mark them with a star: rj*. 
The hadronic matrix element {Qs2) is parametrized as 

{K^\Qs2{l^)\K') = lBK{pi)fWK- (5.48) 

Here and fx are mass and decay constant of the neutral K meson. In (5.48) Bxi/J') 
measures the deviation from the vacuum insertion approximation Bk{ij) = 1. It must 
combine with h{^) in (5.47) to the RG invariant 

Bk = BMKpl), (5.49) 

because physical observables are scale independent. 




5.5. Why Working Beyond Leading Logarithms? 

After all we may ask whether it is necessary to do the effort of a NLO calculation of 
K*^ — K° -mixing instead of contenting oneself with the leading log approximation. Let us 
therefore close this chapter by listing the compelling reasons for passing to the next-to- 
leading order: 
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i) Only in NLO calculations it is possible to use the NLO quantity Aj^ defined in 
(3.15). 

ii) The dependence of f]3S{xc,Xt) on the top-quark mass first enters in the NLO. 

iii) The LO results show a huge dependence on the scales fiw and fic at which particles are 
integrated out. Such scale dependences are inherent to any RG improved calculation 
working with truncated perturbation series. It is caused by the fact that we can always 
subtract an arbitrary small logarithm (such as \n{iiw/Mw)) from the summed large 
logarithm. By this one pushes small terms from the LO summation into the NLO. In 
the NLO the scale dependences therefore diminishes. 

iv) One must go to the NLO to judge whether perturbation theory works, i.e. whether 
the radiative corrections are small. Further the corrections can be sizeable. 
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6. Calculation of t/i and 7/3 



The following chapter is devoted to the description of some details of the NLO calculation 
of 773 and Tji and to the presentation of the analytical results. 



The problem of RG improved short distance QCD correction to K — K'^ -mixing has been 
first addressed by Vamstem, Zakharov, Novikov and Shifman and by Vysotskii [55] . They 
did not use the operator product expansion approach, but extracted the coefficients of the 
leading logarithms in a complicated way from the Standard Model amplitudes. Pioneering 
work in the development of the factorization method has been done by Oilman and Wise 
[45] . Their leading log results for rji and 772 agreed with those of Vysotskii, but they found 
his result for 773 to be incorrect. In these early works the top-quark has been treated as 
light. Later Flynn [56] has redone the calculation for the case of a heavy top-quark. A 
good approximation neglecting flavour thresholds has been derived by Datta, Frohlich and 
Paschos [57]. 

Let us comment on the LO calculation of 773 of Oilman and Wise here: Their approach to 
the RO evolution is equivalent to the one described first in the paragraph on if"^* in sect. 
5.4.1, i.e. they consider the mixing of TZij into Q7 and work with coefficient vectors similar 
to (5.32). Yet their bilocal structures are defined in an inconvenient way, so that they 
need four evolution matrices instead of two. Further their operator basis is ovcrcomplctc 
resulting in 8 x 8 matrices rather than 7x7 ones. Indeed, the evolution matrices in [45] 
have a double eigenvalue. 

With (3.45) wc can simply read off the LO anomalous dimension tensor from the divergent 
parts of (TZij)^^^ depicted in fig. 5.4 and fig. 5.12. One finds 



6.1. The LO Analysis 



/ -4(7V + 1) \ 



/ 4(iV-l) \ 



-8 





8(iV- 1) 





-8{N + 1) 




(6.1) 



-16 



8(iV + l) 

V 16 ; 



-8(A^- 1) 

V / 



We have expressed the result of [45] in terms of (6.1) and found agreement. 
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6. Calculation of 771 and 773 



6.2. NLO Strategy 

From the discussion of chapter 5 we know that we need the divergent parts of two-loop 
diagrams for the calculation of the NLO anomalous dimension tensor needed for 773 and the 
finite parts of two-loop diagrams for the matching calculations in 771. According to sect. 3.7 
and 3.8 we can use any external state to perform the calculation and we choose massless 
external quarks with zero momentum. The QCD gauge parameter ^ is kept arbitrary. Its 
absence from the anomalous dimension tensor and the Wilson coefficients provides a check 
of the calculation. The same has to be true for the IR regulators, which are small masses 
rus and nid for the internal quarks. The problem therefore requires the calculation of the 
diagrams of fig. 5.5 and fig. 5.13. We remark here that the insertion of the right-handed 
penguin operators Q5 and Qq into the diagrams of fig. 5.13 involves different integrals than 
that of the left-handed ones Q3 and Q4. 

A further check of the result will be the absence of non-local terms involving In /j, after the 
inclusion of subloop counterterm graphs. 

Yet an incorrect treatment of the evanescent operators would pass all the checks mentioned 
above. In the calculation of 7^3 we will therefore keep Si and bi parametrizing the order-e 
term in the definition (4.6) and (4.16) of -Ei[(57] and E2[Q7] arbitrary. While the individ- 
ual Z-factors in (4.18) depend on 61 's in (4.16) this dependence cancels in 7*^^). We will 
insert the evanescent counterterms with a factor of A and with A = 1/2 the dependence 
on bi correctly vanishes. Conversely the non-local terms in the two-loop renormalization 
constants zjp^^^ responsible for the mixing of physical operators into evanescents vanish 
for A = 1. 

Clearly the task at hand is simplified by an algorithm for the calculation of the two-loop 
diagrams. For the Dirac algebra the computer package Tracer [58] has been used. Further 
two-loop tensor vacuum bubble integrals up to rank six are involved. An algorithmic way 
to calculate them is described in the following section. 



6.3. Master Formula for the Two— loop Integrals 

We define the general two-loop vacuum bubble tensor integral as 

n]^M"'''-''^{{mj},{Mij},{M2j}) = 

^4£g27Be rr rj p^^-'-p"- ■k^'^■■■k^''■ 

JJ ^ {ml - A;2) • • • {ml - k^) ■ {M^, - P^) ■ ■ ■ {Mfj, - p^) 



1 



{Mi-{p+kf)---{Mi^-{p+kn 



(6.2) 



Here the Mi^'s and M2j's are arbitrary masses, while the m^'s are small masses to regulate 
IR singularities. This is the most complicated two-loop integral which can appear in the 
calculation of QCD corrections in fight hadron systems. In our case of K° — K° -mixing all 
cases with < L < 4 show up in the two-loop diagrams. Let us set T small masses exactly 
to zero while keeping L — T of them nonzero. In practice one chooses T just as large that 



6.3 Master Formula for the Two-loop Integrals 
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Figure 6.1: The vacuum bubble diagram with two light masses mi and m2- The tensor 
structure is arbitrary. 



all the remaining L — T masses suffice to regulate the IR singularities. In our calculation 
we needed at most two nonzero light masses, which were quark masses and 171^. But 
also the case of on-shell dimensional IR regulators will be contained in the result, it simply 
corresponds to L = T. The heavy masses My and in (6.2) correspond to the remaining 
masses in the problem, i.e. M^y, rrit, rric and in the Standard Model diagrams and rric 
and m„ in the effective theory diagrams. We have only needed the latter, but for different 
combinations for the internal masses, so that we also keep them arbitrary. 
For the scalar case r = s = the result for (6.2) is known for all masses being arbitrary [59], 
it involves a F4 function and is not very useful for the situation at hand. 
Finally the prefactor in (6.2) has been introduced for later convenience. 
Let us now sketch the solution of (6.2): The first step consists of a decomposition into 
partial fractions of the propagators involving heavy masses: 

n'NM"''-'''{{m,},{M^,},{M2,}) = 



M N 

EE 

1=1 1=1 



M 



n 

k=l 



Ml - Ml, 



N 

n 

fe=i 



rpUl...Us,fJ,l. 



{{mj},Mu,M2i) 



I 



For equal heavy masses one finds derivatives of TjX\ 



V\...Vs,i),\...ii.r 



with respect to Mu or M; 



2Z 

instead. jg depicted in fig. 6.1. We refrain from using partial fractions also in 

the line fiown through by the loop-momentum /c, because this would yield inverse powers 
of the infrared regulators ray 

The first step to solve 2^^||^--''s,Ati---A'r ^j^^ expression of the sub-loop integral over p in (6.2) 
in terms of a Feynman parameter integral. The corresponding expression for a general one- 
loop tensor integral has been derived first by Davydychev [60]. Clearly only one Feynman 
parameter x is involved. Now one is faced with a sum of tensor integrals corresponding to 
all possibihties to build a rank-s tensor out of A copies of the metric tensors and s — 2A 
copies of the other momentum h. Each of the remaining integrals over k is proportional to 
the totally symmetric tensor 



_^...Vp^,jX\...ilr 



where {p2A+i ■ ■ ■ Ps) is a subset of {1, . . . s). The result reads 



rpV\...Vs,ll\...ilr ,, 

-'lh — a* 



[i] 

E E M^Y' 

A=0 {pi,...p2A}U{p2A+l,-Ps}={l,-s} 



(-1) 



3 o_i_r 
2 



r(-A + £)e2"^E 



r(^ + 2- A 







-2A 



\x{\-x)\ 



X-6 
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6. Calculation of 771 and 773 



il+l±l_A-T,A(^(^)' "^T+l, ... mi). 



Here the second sum runs over all ( ^ ) possibilities to select p^A numbers from the set 

V P2\ J 

{1, . . .s}. Further M'^[x) = M^Jx — M|;/(l — x) and the integration over z = \/ —k'^ i 
contained in 



IS 



Here the original IR singularities show up at the lower end 2; = of the integration in 
the limit 0. We have put T masses exactly to zero with T < 1 + [(r + l)/2], so 

that the remaining masses regulate the IR singularity at z = 0. Wc can disentangle the 
IR regulators from the integral by noticing that the integration path in (6.3) connects two 
branch points of the integrand. Hence we can express the integral in (6.3) by a contour 
integral 

/ dzz^-^f{z) = y / / dz{-zy-^f{z), 
Jo 2z sm(7r£) Jc 

where the contour C starts below the real axis at 2; = +00, encircles the point z = 
clockwise and returns to z = +00 above the real axis. Here C is arbitrary as long as the 
poles z = —m^j^^, . . .z = —m\ and the third branch point z = —M^ are to the left of 
C. Next one deforms C such that it encircles the other branch cut in (6.3) connecting 
z — — with ^ = 00 along the negative real axis. When doing this one has to pull C 
over the poles thereby picking up 27ri times the sum of the residues of the integrand at 
z = —m'^_^i, . . . z = —m\. One then obtains: 

jUM,mr^.,...m,) = {-If^^ dz + + z) ■ ■ ■ {-ml + + z) 

^ sin(7r£),i;, n^^^^ {ml-mf) ' ^ " ^ 

k^l 

In the first integral we can safely set the small masses to zero and express the integral in 
terms of Euler's beta-function. This integral represents the result for the case of dimen- 
sional IR regularization. All IR regulators rrij are contained in an additive term. This must 
be so, because we know that IR terms in the full and effective theory factorize with the 
LO Wilson coefficient. If one kept small external momenta as IR regulators, they would 
appear in the same way in an additive term to the first integral (with rrii = 0) in (6.4). 
Now we can do the final integration over the Feynman parameter x yielding: 

mr"'''-''{{mj},M^j,M,j) = 
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[f] 

A=0 {pi,...p2A}U{p2A+lvPs} = {l,...s} 



^5 A 1 



~2~ 



-1) 



r — s 
2 



r(-2 + L - ^ + 2e)e^'^^ 



(2 



A 



"'^+L-2-A,- 



+ 7^LitlR({^.},M,,M2,) 



(6.5) 



Here T^n^R is the part comprising the IR regulators. In (6.5) {k)i = k-{k + l) ■ ■ ■ (k + l — l) 
is the Pochhammer symbol. The function h in (6.5) reads: 

hl^(M,,Mj) := i,'^['dxx'+%l-xr+^[x(Mf-M^) + M^Y~''^hl,(Mj,M,) 

J 



'Ml 



-2s 



(M2)f ( --f ) B(l + A; + £, 1 + n + £) 



M.2 



2Fi(-p + 2£, 1 + A; + £, 2 + A; + n + 2£, 1 - — |). 



(6.6) 



The second expression in terms of GauB' hypergeometric function is useful for the case 
of one small mass Mj, e.g. = m^, in which case one can use Kummer's relations 

to transform to 2Fi{., ., .; /Mj). In general the integral representation is more useful, 
because it can be more easily expanded in terms of e. It further exhibits more symmetries 
and functional relations such as 

^ ;hl^(M„Mj) = (p-2s)hl;^\,(M„Mj). (6.7) 



For equal masses one simply has 
hljM,M) = 



H^' (^M^y B{l + k + e,l + n + e). 



For different masses the expansion of the integral in (6.6) in terms of e turns out to be the 
only cumbersome task to be done. One finds the familiar dilogarithm in the finite part and 
higher polylogarithms in the higher order terms in e. Yet one only has to expand a few 
of the h^s and one can obtain the others by the use of the functional equations (6.7). The 
expanded /i-functions have been stored in a Mathematica program to build a database. 
Now for a dimensional regularization of the IR singularities one finds T^fi'iR = 0. For the 
use of small masses as done by us one obtains: 



■1) 



2£7E,,4£ L-TL-2 2" 



r+i-r-A r(— A + £)e p \ - \ - 

sin(7r£) r(^ + 2 - A - £) ^ 



+A 



l+£±r_T-A+i-£ 



n 



L-T 
n=l 



mt 
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6. Calculation of 771 and 773 



' I \ / -,\k+l 



k 

k=0 



E i (-1 



(A + 1 - £) s+k+l-2X 



2 ..2 ■ (6-8) 



In practice the summation indices only run over a small range, in the calculation of rji and 
?73 only situations with L — T < 2 and k = I = appeared. 

We remark that the master formula (6.5) is a good investment into the future of short 
distance QCD corrections: First the formula holds in D dimensions. It can therefore be used 
for the matching calculations in a NNLO calculation and also for the two-loop counterterm 
diagrams in a n-loop calculation, which requires two-loop integrals up to the order n — 3 
in e, if the divergent part is needed. Probably all calculations beyond NLO involve so many 
diagrams that a complete computerization is unavoidable. Then formulas of the type (6.5) 
also for higher loop integrals are indispensable. Further (6.5) can simply be extended to 
more powers in the small masses. Such integrals are needed in NLO calculations for the 
coefficients of higher dimension operators. 



6.4. Result for the NLO Anomalous Dimension Tensor 



In this section we hst the result for the NLO anomalous dimension tensor 

From (5.28), (5.29) and (3.28) one readily finds the relation between the two-loop renor- 



malization tensor 



Z and the bare 0(a) matrix elements of TZij-. 



1.3, 



(1) 



(1) 



33 



1 (0),bcire 



■1,(1) 



^3,1 



.(0) 



■1,(1) 



(6.9) 



where wc have used the notation (3.52) and (2.35) for the expanded matrix elements and 
Z-factors. Now one has to take into consideration that 



(Qr) 



bare 



,m 



{- 



c, bare ^ \ bare 



5'bare 



Qs2r 



y H' 



bare 



(6.10) 



Hence 



(Qr) 



(l),bcire 



(2Z« + /3o) -{Qs2r + -^-{Qs2)^'^''-. (6.11) 



Now the divergent parts of the diagrams of fig. 5.13 and fig. 5.5 and the corresponding 
subloop diagrams yield the terms in the first line of (6.9) and the last term in (6.11). The 
remaining divergences therefore correspond to 



+ 



■1,(2) 



1 



(6.12) 
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Hence we can simply read off ^'^^^ from the 1/s-divergences of the diagrams of fig. 5.5 
and fig. 5.13 after the inclusion of subloop counterterms. We have also found that the 

dependence on ^ vanished from the 1/e^-terms in (6.12) after subtracting the term involving 
Further it has been checked that Z2 fulfills the relation to the \ je divergences of 
the one-loop renormalization constants. 

Now 7([.2,7 is related to Zi^'^^^ via (3.45), i.e. the evanescents must be inserted with 
a factor of 1/2. 

To show the scheme dependence found in chapter 4 we keep the definitions of some evanes- 
cents arbitrary: 



El[Ql/2] 

El [Qs/e] 
EiiQj] 



l^lul^fL -f^^Yl^L - (4 - 8s)^^L -f^L 
l^lul-dL ® -f'^Yl^L - (4 + aie)-f^L -ff^L 
It^lvl^R ® I'^Yl^L - (16 + 026)7^^? ® l^L. 
lnluli)L ® YYi^L - (4 + ai£)7^L ® 7'^L 



(6.13) 



Insertions of Qi and Q2 into the diagrams of fig. 5.5 also involve 



Now one finds for N — 3: 
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404 ^ 130ai 
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32ai 
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~3~ 
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88ai 32ai 
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80ai 
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+ 
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3 3 3 

1576 SOai 325i 802 
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3 3 3 3 

1664 28ai 28ai 52a2 



3 / 



3 3 
52 llOSi 

3~ 

-60 + 4ai 

1592 8ai 16ai 
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3 3 3 

-88 + 8ai + 4ai 



1160 16ai 16ai 802 
i _^ i _^ 1 

3 3 3 3 
V -128 - 4ai - 4ai - 4a2 } 



he -f^L^-ff^L. (6.14) 



(6.15) 



88 
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The dependence on 61 cancels from (6.15), although it is present in the individual contri- 
butions to the first two components of 7^ 7. Yet we discover a dependence on the a^'s and 
oi as stated in (4.43) and (4.44). 
For the numerical analysis we will use 

Qi — —8, 02 = —16, 5i = —8. (6.16) 

The first two choices have to made to comply with [44], from which we have taken the 
I AS| = 1 anomalous dimension matrix and the initial condition for the | AS| = 1 coefficient 
vector. 



6.5. T]3 

We describe the initial condition and the RG evolution for the formalism of sect. 5.4.1 

— * 

using D according to (5.36) and (5.38). All formulas refer to the NLO. The LO expressions 
can be simply obtained from them by only keeping the leading terms and using the LO 
running a. 



6.5.1. Initial Condition 

The initial condition for the | AS | — 1 coefficient vector can be taken from [44] : 
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B2 
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x{l%-\\x-x^) 2 
12(1 -x)3 3 



(6.17) 



In (6.17) we have introduced a term involving the second row of the LO |AS| = 1 anomalous 
dimension matrix 7^°^ to allow for the choice ^ Mw 

— * 

Further D in (5.36) involves the initial coefficients of the operators (5+ and 



C+{iJiw) - 5i + B2 + 7rin 



(0), /^vy 



Ml 



w 



C4fxw) = -i^i + i^2 + 7i°^ln^.(6.18) 



Ml 



w 
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For the initial coefficient of C7 we have to perform the matching calculation with the 
diagrams of fig. 5.1 and fig. 5.4. If the evanescent operator of is defined as in (6.13), 
one easily finds: 



C7{lJ,w) = 



Mw ^ 

where F{xt) is the top-dependent part of S{xc,Xt): 

S{xc,xt) = -Xclnxc + XcF{xt). 



(6.19) 



6.5.2. RG Evolution between /xw and fXc 

The 8x8 anomalous dimension matrix 7 has been introduced in (5.38): 

/ 7 \ 



7 = 



T ~ 

r+,7 7 





- 7+ 
7 - 7- / 



(6.20) 



While the |AS| = 1 mixing matrix 7 can be taken from [44] and 7± and 7^7 have been 
given in (5.19), (6.1) and (6.15), we still have to list the anomalous dimension 7 of Qj. 
Here one has to take into account the running of the mass and coupling in its definition: 



7 



7+ + 27^ - 2/3o. 

The anomalous dimension 7^ of the mass has been given in (3.11). 

Now one can build the matrices U^^^ and J^'^^ as described in (3.40) and (3.43). The NLO 
evolution matrix (3.42) from nw down to Hc reads: 

a{iJ,b) 



(6.21) 



^ ^ Q(/^c) -j(4) 

An 



_ o\pw)_-j{5) 

An 



1 ^ (ji^) _ j(4) ^ _T(5) _ ^T(4) 

An V 



An 



(6.22) 



Here the matching correction at the b-threshold has been taken into account, r*^-^^ coincides 
in the upper left 6x6 submatrix with the result of [44]. The 7th and 8th row and column 
are zero, because there is no b-quark involved in the diagrams of fig. 5.5 and fig. 5.13. 
Now after obtaining D^/^dic) in (6.22) we can find the coefficient C'j{iJi,^: 



(6.23) 
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6.5.3. Matching at /ic 

The matching condition at fic has been given in (5.44). Now both the diagrams of fig. 5.4 
and those of fig. 5.12 have to be calculated. One finds for the quantities rij 32 defined in 
(5.43): 



[-41n(me//iJ + (3 + ai/2)]Ti,- 
[-8 ln(mc//Xc) - 4] Tij 
[81n(mc/Aic) + 4] 



for j <2 
for 3 < j < 4 
for 5 < J < 6 



(6.24) 



with the colour factors 



r±i = r±3^r±5 = -{1±N), 
T+j — 1 and T_j = for j even. 



Finally with (5.46) we write down the NLO result for rj^: 



(6.25) 



S {Xc{llc):Xt{llw)) 



,(3) 



1 47r 



4q;(/Xc) 



1 - 



1 ^ 



(1 + 0(a^ 



(6.26) 



Order a terms should be consistently removed here. 

In most cases it is more useful to define 773 such that it multiplies S{x*^^x*^) in the hamil- 
tonian, where = xj^mf) and x\ = Xt{mt). We then mark 773 with a star: 



S{xl,xt] 



[a(/^c)]' 



1 Att 



4 a{iic) 



1 - 



1 



+ 9 H II?^ii,S2Ci(/ic)Cj(/ic 
^ i=+-j=l 



(l + 0(a^ 



(6.27) 



I.e. one has 



rjsS {xc{i^c),xt{f^w)) = r73,S(x*,x*) . 



(6.26) has been checked by expanding all evolution matrices to the first order in a to verify 
that the leading logarithm in S{xc, Xt) is correctly reproduced. Indeed it has been found 
that ri3b{n) = l + 0(a), when the RG improvement is switched off in this way. Further rj^ 
has been expanded around /ic — ttIc to check whether the logarithm in (6.24) cancels the 
dependence of the LO term on Hc in the calculated order. This is also the case. 



6.6 r]i 
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6.6. r]i 



As discussed in sect. 5.4.1 the RG factors in rji only involve the known RG evolution of C+ 
and C_. The only complicated task is the matching at the scale /^c shown in (5.41), which 
requires the finite parts of the diagrams of fig. 5.5 and fig. 5.6. We will be sketchy here, for 
more details on the calculation we refer to [17] and [61]. 
The LO matrix element of Oij reads 



with the colour factors 



-A^ + 1 



(6.28) 



A^-1 



In the NLO one obtains: 



^ [(g52)^°)a?^'V) + {T)^'^njhT + {U)^'^njhu] . (6.29) 



The coefficients multiplying the unphysical matrix elements are of course the same as those 
of the Standard Model amplitude given in (5.13). We split ay"''* in (6.29) into its physical 



part cy and those parts which depend on the infrared regulators or involve the gauge 
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The desired physical parts in (6.30) are found as 
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The matrix element of Qs2 can be found in [15,43]. Writing 

{Qs2 (//)) = {Qs2 (//)>(°) + ^ [{Qs2Y'^a{tJi) + (f')(°)/i^ + (C/)(°)/.c7] (6.32) 



the coefficient a(/x) reads 
' N -1 



c + 



2Ar 



3 log 



Cf + 



N -1 
2N 



2-2 



log(m^//i2) - log(my//2) " 



m3 



iV 



2N 



1 , mimi 



-3Cf -5 



A^- 1 
2A^ ' 



(6.33) 



Now we can solve the matching condition (5.41) for C'gldj.c): 

i,j=+,- 

with 

Again the unphysical terms have canceled from the Wilson coefficient. 
Hence (5.46) gives for rji. 



(6.34) 



(6.35) 



^(3) 

(«(/^c)) + 22 C'i{nc)Cj{nc) 
hj=+,- 



An 



(6.36) 



r]l is simply given by 



m 



*2 



(6.37) 



We close the section on r]i by looking back to the matching at the scale /iw as displayed 
in (5.30) to analyze the findings of [39] and sect. 5.4.1: 

Consider first the diagrams Di and D3. They can only be divided into their |AS| = 1 
substructures by a thrcc-particlc cut. This cut separates two |AS| = 1 subdiagrams with 
four external (anti-) quarks and one external gluon. The cut is flown through by two 
integration momenta, but both refer to loop integrals which are finite after using the GIM 
mechanism, which reduces the UV behaviour of the integral involving the internal up- 
type quark propagators. Hence these diagrams match their counterparts in fig. 5.2 with 
the tree-level Wilson coefficients Cl^^Cf'\ The same reasoning applies to the self-energy 
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insertion in diagram D5. Now the remaining diagrams can be divided by a two-particle cut 
into their |AS| = 1 substructures. Again the divergences of the loop integration associated 
with the momentum flowing through the cut are canceled by the GIM mechanism. But 
now the cut separates a tree-level |AS| = 1 vertex from a dressed one, which involves an 
UV-divergent loop. These UV divergences are responsible for the 0{a) corrections to the 
|AS| = 1 coefficients C±{iJ,w)- In fact the colour singlet diagram Dq does not contribute to 
the |AS| = 1 Wilson coefficient. We have verified that only D2, D4 and contain terms 
requiring C^^\fiw) to match their Standard Model counterparts. Finally only Dj gives a 
term involving In fi]^ / M\.y , which appears in the NLO matching to compensate the LO 
scale ambiguity. This can be traced back to the fact that involves diagram (3) of fig. 4.1 
as a subdiagram, which is the only source of the LO anomalous dimension of the | AS| = 1 
operators. 
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7. Numerical Results and Phenomenology 



7.1. r]3 

At first we will investigate the dependence of rjs and 773 on the various physical parameters 

and on the scales at which particles are integrated out. 

Recall the relevant part of the effective low-energy hamiltonian from (5.40): 

= ^M^2X,XtvlS{xl,x;)b{i,)Qs2{l^). (7.1) 

Wc will start with the discussion of the scale dependence: 773 and rj^ depend on the initial 
scale fiw and on fic ~ itLc, at which the c-quark is removed from the theory. There is also 
a dependence on the scale /ib ~ mb, at which we pass from a five-quark theory to a four- 
quark theory. But /ib only appears in the running a and the dependence is so small that one 
could even set /ib = A*c- Now rjs in (7.1) multiplies S{xc{iJ^c),Xt{lJ'w)), which also depends 
on the scales, and ideally the product is scale independent. Clearly for the discussion of 
this issue the notion of rj^ is more useful, because it multiplies a scale independent function 
and should therefore also be independent of nw and Hc- In sect. 6.5 we have mentioned 
that the LO scale dependences of 773 are analytically canceled by the NLO terms within 
the calculated order. Hence wc can use the remaining scale dependence as a measure of the 
accuracy of the calculation. Since W-boson and top-quark are simultaneously integrated 
out at /X = iJ,w, the range for is 

Mw ^ l^w ^ ruf 

Wc expect the dependence of 773 on ~ ttIc to be larger than the one on /iw because of 
the larger coupling constant involved. 

Next consider the dependence on the QCD scale parameter Aqcd: In the following Aqcd is 
understood to be defined with respect to four active flavours, the corresponding quantities 
in the three- and five-flavour theory are obtained by imposing continuity on the coupling 
at fib and /j^. The world average for a{Mz) = 0.117 [11] correspond to Aj^ = 310MeV for 
lib = 4.6 GeV and Aj^ = 315MeV for /ib = 5.0 GeV. The leading order Aqcd, however, 
differs from Aj^ by an overall /^-dependent factor, see (3. 15) . If one equates the LO coupling 
and the NLO MS-coupling constant at the scale fj, — Mz — 91 GeV, one finds that — 
315MeV corresponds to Aqq^ = llOMeV. If the relation is imposed at the low scale 
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fjL — 1.3 GeV, one finds Aq^d = ISOMeV. This shows one contribution to the error bar in 
LO calculations. It is, however, more pronounced in rji than in 773. 
Let us therefore pick the following set of parameters: 

Ams = 0.31 GeV, A^Oj, = 0.15 GeV, 

me = //c = 1-3 GeV, = 4.8 GeV, (7.2) 

Mw = 80 GeV, mt{mt) = 170 GeV, nw = 130 GeV. 

The value for 77^ corresponding to this set reads: 

rifo = 0.365, rif^° = 0.468. (7.3) 

Hence the NLO calculation has enhanced r/3 by 27%. From the difference of 0.103 between 
the two values in (7.3) 0.019 originates in the change from the LO to the NLO running 
a. The smallness of this contribution is of course caused by the adjustment of Aq^j^ to 
fit the NLO running coupling as described in the previous paragraph. The explicit 0(a)- 
corrections from the NLO mixing and matching contribute 0.084. Let us list the dominant 
sources of the enhancement: At the initial scale nw the coefficients C2, C+ and C_ have 
a size of the order 1, while all other coefficients are almost negligible. The RG evolution 
from to He enhances the coefficient C_ by roughly 75 % because of the negative sign 
of the anomalous dimension of while the coefficient C+ is damped by 25%. Now the 
penguin coefficients C3 to Ce are still negligible at /i — /ic, only Ci(/ic), C2{iJ>c) and Cj^ic) 
are important. In the matching at He the contribution of C+(//c)C'i(//c) numerically cancels 
the one of C^{fic)Ci{fic), because the RG damping of the former is compensated by a 
larger colour factor r+i = 2 having the opposite sign of r_i = —1. Finally CT^fic) ~ 0.7 has 
become large due to the RG admixtures from C2. Hence only C_(/ic)C'2(/ic), Cj^{iic)C2{Hc) 
and C'j{hc) are important. In LO only Cjdic) enters 7)3. Keeping only C'j{hc) in the NLO 
expression, however, overestimates the NLO enhancement by a factor of roughly 1.5, be- 
cause C2(/ic) contributes with a negative sign to 773 (for the standard definition of the 
evanescents: ai = —8), see (6.26). 

Next we will display the dependence of 773 and 773 on the various parameters. In all plots we 
take the fixed quantities from the standard data set (7.2). First consider the dependence 
of on the initial scale as depicted in fig. 7.1. We find a sizeable scale dependence of 
14% in the LO result. It is almost totally removed in the NLO, where we find less than 
3% change in r^g, when fiy/ is varied between 70 GeV and 190 GeV. The situation is not so 

Figure 7.1: Dependence of r/g on the scale hw, at which the initial condition is defined. 

nice in the case of the dependence on ^c- To show the large effects related to the running 
c-quark mass in (7.1) we first display the running of 773 with ^c- In (7-1) it multiplies 
S{xc{[jic)iXt{Hw))i which grows with increasing Xc and therefore falls off with increasing 
He- Hence 773 has to grow with increasing ^c- From fig. 7.2 one realizes that the LO result for 
773 worsens the effect of the running mass in S{xc{Hc)-iXt{Hw))-i while the NLO 773 correctly 
grows with increasing ^c- 
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Figure 7.2: Dependence of rj^ on the scale fic, at which the charm quark is integrated out. 
rjs has to compensate the large running of the charm quark mass. See also fig. 7.3. 

Figure 7.3: Dependence of 773 on the scale Hc- For the scales below 1 GeV one recognizes 
the breakdown of perturbation theory. For a realistic estimate of the scale ambiguity take 
1.1 GeV < //c ^ 1.6 GeV. 



The corresponding plot for r)^ is shown in fig. 7.3. We have intentionally extended the 
range for to the unphysically low value of ~ 0.7 GeV to visualize the breakdown 
of perturbation theory. Varying within the range 1.1 GeV < jic ^ 1.6 GeV yields the 
following estimate of the scale uncertainty: 

0.33 <r]l^^ < 0.40, 0.43 < 773^^° < 0.50. (7.4) 

This corresponds to a reduction of the scale uncertainty from 19% to 15%. One reason 
for the poor improvement is the fact that the NLO running of the mass is stronger than 
the LO one. One reahzes from (7.4) that the scale ambiguity alone is not always a good 
measure of the accuracy of the calculation, because the central value for r]'^^^'^ is not in 
the range quoted for t^^^*^. Yet in the NLO order one may also judge the contribution of 
the uncalculated 0(Q;^)-terms by squaring the calculated 0{a) corrections. By this one is 
lead to the same interval as in (7.4), so that we may consider (7.4) as a realistic estimate 
for 77*^^0. 

Let us now look at the dependence of 7/3 on the physical parameters. From the smallness 
of the initial coefficient Cj we expect the QCD correction to be almost independent of 
ml = mt{mt). The relevant quantity to be discussed is the product ri^S{xl,xl) which 
enters (7.1) and physical observables. It is shown in fig. 7.4. Clearly the LO curve is fiat. 



Figure 7.4: Dependence of r)lS{xl,x1) on ml = mt{mt). There is no dependence on mt in 
the leading log approximation. 

We also display the dependence of 773 alone on mf, see fig. 7.5. One cannot get much physical 
insight from this plot, it only serves to verify that we can treat 773 as a constant function 
of mt in phenomenological analyses. 

On the other hand there is a sizeable dependence of 773 5* (x*, xl) on 7n* = mdmc) as shown 
in fig. 7.6. In 77^ alone, however, the LO dependence on m* happens to be washed out in 
the NLO as can be seen in fig. 7.7. We close this section with a plot of the dependence 
of 7/3 on Aqcd- fig. 7.8 reveals a very moderate dependence on the QCD scale parameter. 
Recall that the actual measurements of a{Mz) corresponds to LO values for Aqcd close to 
the left edge of the plot, while the NLO scale parameter A^^ is close to the highest values 
displayed. 

Now we can easily summarize the result: The largest uncertainty of our estimate of 773 is 
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Figure 7.5: Dependence of 7^3 on m^. 
Figure 7.6: Dependence of ri^S{xl, x^) on m* 
due to the the scale ambiguity stemming from /Xc and we obtain: 

*NLO n +0-03 



*^i.u _ 0.47^""" (7.5) 

" -0.04 ^ ' 



Figure 7.7: Dependence of rj^ on m*. 
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Figure 7.8: Dependence of rj^ on Aqcd. The NLO curve corresponds to Aqcd = ^ms- 
Realistic values are 300MeV < A^jg < 330MeV and lOOMeV < A^^^ < 200MeV. 



7.2. rji 

Again we only sketch the NLO result for rji and refer to [17,61] for more details. The part 
of the effective hamiltonian involving rji has been found to be 

- Xlvtm:%{pi)Qs2{pi). (7.6) 



167r2 " " 

For the set of parameters given in (7.2) one finds: 



r^*LO ^ o_go +0.20 ,NLO ^ +0.21 _ .7^. 

-0.16 -0.23 ^ ^ 



For m* = 1.4 GeV one hkewise finds 



c 



*LO ^ Q_78 +0.17 .NLO ^ 122 . (7.8) 

-0.14 ' '1 -0.18 ^ ' 

Here the the dominant source of the error is again the dependence on fic, which has also 
been varied between 1.1 GeV and 1.6 GeV in (7.7) and between 1.2 GcV and 1.7GeV in 
(7.8). Now the dependence on Aqcd is more pronounced than in the case of r]^. Therefore 
we have included the variation with Aq^d for llOMeV < Aq^j^ < 180MeV in the error bar 
for the LO result. The dependence on Aqcd is displayed in fig. 7.9. The NLO calculation 
has therefore yielded a reduction of the theoretical error from 45% to 33% for m* = 1.3 GeV 
and from 40% to 26% for m* = 1.4 GeV. 



Figure 7.9: The dependence of rjl on Aqcd, ^ms- "^^^ choice m* — lA GeV has been made. 
All other fixed parameters are as in (7.2). 

(7.7) reveals a much larger scale dependence for 77* than we have found for r^gin (7.5). Let us 
therefore investigate the origin of this uncertainty: There are three different contributions to 
rjl in (6.36) stemming from the products C+{^c)C+{f^c), C+{fic)C-{fJ'c) and C^{fic)C-{f^c)- 
The last one is clearly the largest, because of the negative anomalous dimension oi Q-. 
These different sources are displayed in fig. 7.10: One can easily see, that the NLO curves 
related to C+C+ and C+C_ are almost flat as they should be, but they contribute with 
opposite sign to r)l and therefore almost cancel in the sum. On the other hand the part 



100 



7. Numerical Results and Phenomenology 



Figure 7.10: The dependence of the individual contributions to ril on fi^ for ^qcdj-^ms ~ 
SOOMeV. m* = 1.4 GeV has been chosen. The line labeled by (all) refers to rjl- 

stemming from C-C- contains a large residual scale dependence and is identified as the 
source of the large theoretical error in (7.7). 

We have already noticed in (7.7) and (7.8) that 77* exhibits a sizeable dependence on m*. 
We therefore plot function of m* in fig. 7.11. Physical quantities involve ml'^rjl, 

which is shown in fig. 7.12. 

Figure 7.11: The dependence of rjl on m*. For very low values m* ~ 1.1 GeV as favoured 
by some sum rule analyses our perturbative calculation of rji is not very reliable. 



Figure 7.12: The dependence of rilm*'^, which enters the hamiltonian (7.6), on m*. The 
(dominant) contribution of rjl to the short distance contribution to the Kl — Ks-mass 
difference is obtained by multiplying the displayed value for 'q{m*'^ by 1.371 • 10~^^ GeV~^ • 
Bk, where Bk is defined in (5.49) 



7.3. Indirect CP-Violation in - RO -Mixing 

Here we briefly describe the indirect CP-violation in the neutral Kaon system (cf. e.g. [62]). 
The immediate consequence of K^ — K^ -mixing is that the weak interaction eigenstates 
and \K^) do not coincide with the mass eigenstates. Yet if the CP symmetry were 
conserved, the latter would be identical to the CP eigenstates \Ki) and \K2): 



K,) = i= (l^O) + |X0)) (CP=+1), 
= ^{\K')-m) (CP=-1). 



\K'. 

Now the mass eigenstates \Ki) and \Ks) differ from \Ki) and |i^2) by small admixtures of 
the other CP eigenstate: 

\K,) + e\K^) \K,) + e\K^) 



2 



The decay of Kl into the CP-even two-pion state may proceed via a CP-conserving decay 
of the small admixture of \Ki) or via a CP-violating decay of the dominant contribution 
of \K'^. These possibilities are termed indirect and direct CP-violation and are compactly 
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explained by the following picture: 
1 



\Kl) = 



Decay: 



CP=-1 



direct 



CP=+1 

indirect 



7r°7r° (CP=+i) 
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A measure of the indirect CP-violation which is independent of phase conventions is given 
by the ratio of the amplitudes 

In terms of the off-diagonal element of the mass matrix 

Mi2 = {K''\H^^=^\W) (7.10) 

Sk reads 

Sk = -^^(ImMi2 + 2eReMi2), 

where ,^ is a small quantity involving the amplitude of A/ = 1/2 transitions and Am is 
the Kl — Ks-mass difference. Experimentally Ek is very well measured [11]: 

Sk = (2.266 ± 0.023) • lO-V'^/^ (7.11) 
7.3.1. Parametrization of the CKM Matrix 

Since one may redefine the quark mass eigenstates by an arbitrary phase factor, there 
are many physically equivalent forms of the CKM-matrix. We adopt the standard con- 
vention, in which Vudi Kis, and Vtb are real and positive. It is further advantageous to 
have a simple parametrization of the CKM-matrix. Here we use the improved Wolfenstein 
parametrization, which is used in [8]. It reads 



VcKM 

with 



-A - iA^XS 1 - f AA2 

\ AX^{l-p-ir]) -AX^-iAX^T] 1 ) 



( X^\ _ ( A2^ 

P = P 1-y , ^ = ^ 



(7.12) 



(7.12) is an expansion in the parameter A = jT^sl + 0{X^). (7.12) is exact to order A^ and 
contains the phenomenologically important terms up to the order A^. One of the relations 
provided by the unitarity of Vckm is 

Vudv::,+v^v:,+vuv:, = o. 

When expressed in terms of the parameters introduced above, this relation corresponds to 

a unitarity triangle in the (p, r7)-plane with edges (0,0), (1,0) and (p, 77). 
One finds to an accuracy of 0.2% [8] for Xj = VjdV*^: 

ImXt = -ImAc = rjA^X^ (7.13) 

CP-violation occurs for 77 ^ 0. Because of (7.13) we may use ImAt to parametrize all 
CP-violating quantities in the Standard Model. 
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7.3.2. Phenomenology^ 

Inserting (5.40), (5.48) and (7.11) into (7.9) yields the following relation between p and 
rj[8]: 



V 

with 



p)A%Sixt) + Po A^Bk = 0.223. (7.14) 



Pq = [vsSixcXt) - i]ix^] — . 
When A is fixed, (7.14) defines a hyperbola in the (p, r7)-plane. 

Now from (7.12) we notice that A is simply related to Vcb via A = iKtl/A^. We remark 
here that much progress has been made in the recent past to determine |V"cb|. Exclusive 
tree-level b ^ c decays allow the extraction from the endpoint of the spectrum, where 
heavy quark symmetry reduces the problem to the determination of a single form factor. 
Here typically values in the range 0.036 < \Vcb\ < 0.050 are stated [63]. The theoretical 
predictions for inclusive analyses still suffer from large scheme and scale dependences, 
because they are proportional to the fifth power of the bottom quark mass [64] . The latter 
can be determined from T spectroscopy [63] or by the use of heavy quark symmetry from 
inclusive D-decays [64, 65] . The different lifetimes of the charmed hadrons, however, make 
the use of heavy quark symmetry troublesome here. Inclusive analyses give lower central 
values for \Vcb\ than the exclusive method. 

To fix the third edge (p, rj) of the unitarity triangle, we need a further input. It is straight- 
forward to verify that the knowledge of [Kib/Kf,! fixes a circle in the (p, rj) -plane: 

If the phase 6 in the CKM-matrix is the only source of indirect CP-violation in the K^- 
system, the equations (7.14) and (7.15) must have a common solution {p,r]). In general the 
hyperbola will intersect the circle in two points to give two solutions. Yet there is a critical 
set {{mt,Vcb,\Vub/Vcb\, 
Bk)}, for which the hyperbola only touches the circle. This set encodes a lower bound 
on any of the four quantities as a function of the other parameters [8,66]. 
Fig. 7.13 shows the hyperbola for the LO and the NLO result for 773. One can see that the 
hyperbola has moved downwards. Therefore the NLO value for 773 allows for a wider range 
of rritiVch-i \ Vub/Vch\ and Bk than the older imprecise LO result. In this sense the Standard 
Model mechanism for CP-violation has been vindicated. For example a future determina- 
tion (mt, Vcb, 
\Vub/Vch\-i Bk) = (170 GeV, 0.040, 0.08, 0.7) would have been interpreted as a signal of new 
physics, if the LO values 771 = 0.78, 773 = 0.37 had been used. With the NLO results 
771 = 1.32, 773 = 0.47 , however, the hyperbola (7.14) intersects the circle (7.15) and there 
is no conflict with the CKM mechanism of CP-violation. 



^for more details see [67] 
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Figure 7.13: Sk'- The hyperbolas in the (p, rj)-plane for rjt^^'^' = 0.37 and rjt,^^^ = 0.47. 
The other parameters are rjl = 1.32, r)2 = 0.57, rrit = 165 GeV, \ Vcb\ = 0.041, rric = 1.3 GeV 
and Bk — 0.7. The circle corresponds to [Kift/Kbl = 0.08. The NLO value for 773 permits 
two solutions for (p, fj) in this example. 

The effect of the NLO calculation for 773 on the bounds on rut and Bk is shown in fig. 7.14 
and fig. 7.15. 

Figure 7.14: Lower bound on rrit from ek for iKib/Kb] = 0.08, Bk = 0.5, 0.7 and 0.9, 
ril = 1.32, ri^ = 0.57, r^l^^^o = 0.47 and 77^^° = 0.37. 



Figure 7.15: Lower bound on Bk from Ek for rrit = 165 GeV, iV^b/V^b] = 0.06, 0.08 and 
0.10, ril = 1.32, 7]^ = 0.57, rj^^^^ = 0.47 and 77^^° = 0.37. 

Table of Im At 

Since Im At enters all CP-violating quantities, we now give a table of the two solutions for 
this quantity for various values of m.t,Vcb, \ Vub/Vcb\ and Bk- The fixed quantities are as in 
(7.2). The corresponding 77j's are 

771 = 1.32, 772 = 0.57 and 773 = 0.47. 

Both solutions are listed, the entries have to be multiplied by 10~^. n.s. means that there 
is no solution. 

Bk = 0.5 and rrit = 160 GeV: 



IK*! = 


0.034 


0.036 


0.038 


0.04 


0.042 


0.044 




Vut 


- 0.06 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 




Vub 
Vcb 


L = 0.08 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


= 0.1 

\Vcb\ 


n.s. 


n.s. 


n.s. 


n.s. 


1.61 
1.75 


1.37 
1.87 



Bk = 0.5 and mt = 170 GeV: 



\Vcb\ = 


0.034 


0.036 


0.038 


0.04 


0.042 


0.044 


I,!)"! = 0.06 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


1!;^' = 0.08 

\Vcb\ 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


1.43 
1.54 


\Vub\ _ 1 


n.s. 


n.s. 


n.s. 


n.s. 


1.41 
1.76 


1.25 
1.81 
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Bk = 0.5 and m* = 190 GeV: 



|v;.6| = 


0.031 


0.036 


0.038 


0.01 


0.012 


0.011 


\ = 0.06 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


ISI - 0-08 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


1.15 
1.49 


\Vub\ _ n 1 
\Vcb\ ~ 


n.s. 


n.s. 


n.s. 


1.35 
1.6 


1.17 
1.67 


1.06 
1.68 



Bk = 0.7 and = 160 GeV: 



\Vcb\ = 


0.034 


0.036 


0.038 


0.04 


0.042 


0.044 


I!;''! = 0.06 

\Vcb\ 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


I,!;"! = 0.08 


n.s. 


n.s. 


n.s. 


n.s. 


1.13 
1.38 


1.01 
1.39 


\Vcb\ - 


n.s. 


n.s. 


n.s. 


1.14 
1.54 


1.03 
1.55 


0.94 
1.54 



Bk = 0.7 and mj = 170 GeV: 



\Vcb\ = 


0.034 


0.036 


0.038 


0.04 


0.042 


0.044 




n.s. 


n.s. 


n.s. 


n.s. 


n.s. 


1.03 
1.16 


ISI - 0-08 


n.s. 


n.s. 


n.s. 


1.19 
1.27 


1.02 
1.33 


0.92 
1.32 


\Vub\ _ n 1 


n.s. 


n.s. 


1.19 
1.44 


1.04 
1.49 


0.94 
1.49 


0.86 
1.46 



Bk = 0.7 and m* = 190 GeV: 





0.034 


0.036 


0.038 


0.04 


0.042 


0.044 


^^^■'\ = 0.06 

^cb 


n.s. 


n.s. 


n.s. 


n.s. 


0.98 
1.06 


0.86 
1.07 


ll;''! = 0.08 

\Vcb\ 


n.s. 


n.s. 


n.s. 


0.97 
1.23 


0.87 
1.23 


0.79 
1.2 


l^ub] n 1 

\Vcb\ ~ 


n.s. 


1.18 
1.28 


0.99 
1.38 


0.89 
1.38 


0.81 
1.36 


0.74 
1.33 
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Bk = 0.9 and m* = 160 GeV: 





0.034 


0.036 


0.038 


0.04 


0.042 


0.044 


#4 = 0.06 


n.s. 


n.s. 


n.s. 


n.s. 


0.93 
1.05 


0.83 
1.04 


l!;"! = 0.08 


n.s. 


n.s. 


1.08 
1.15 


0.93 
1.2 


0.84 
1.19 


0.76 
1.16 


\Vub\ _ n 1 
IK6l ~ 


n.s. 


1.09 
1.29 


0.95 
1.34 


0.86 
1.33 


0.78 
1.31 


0.72 
1.28 



Bk = 0.9 and m* = 170 GeV: 



\Vcb\ = 


0.034 


0.036 


0.038 


0.04 


0.042 


0.044 


= 0.06 


n.s. 


ri . s . 


n.s. 


n.s. 


0.84 
1.01 


0.76 

0.99 


^^Y"] = 0.08 

\Vcb\ 


n.s. 


n.s. 


0.96 
1.14 


0.85 
1.15 


0.77 
1.13 


0.7 
1.09 


\Vub\ _ n 1 
\Vcb\ - 


n.s. 


0.99 
1.27 


0.87 
1.29 


0.79 
1.27 


0.72 
1.25 


0.66 
1.21 



Bk = 0.9 and = 190 GeV: 





0.034 


0.036 


0.038 


0.04 


0.042 


0.044 


l!;"! = 0.06 

\Vcb\ 


n.s. 


n.s. 


n.s. 


0.8 
0.93 


0.71 
0.92 


0.65 
0.89 




n.s. 


0.94 
1.04 


0.81 
1.07 


0.72 
1.05 


0.66 
1.02 


0.6 
0.98 




0.97 
1.15 


0.83 
1.2 


0.75 
1.19 


0.68 
1.16 


0.62 
1.13 


0.57 
1.09 



7.4. Kl — Ks— mass difference 

We will now discuss the short distance contribution to the Kl — Kg -mass difference: Ne- 
glecting the small imaginary parts of the CKM elements, one finds 



(Am) 



SD 



6^ 



mKfKBK (Re \cYnil% + 2(Re A^)(Re \c)M^S{xl, x*M 



+ {Re\fM^S{xl)r]l. 



(7.16) 



Here the relative size of the three terms in the brackets is roughly 100 : 10 : 1. Hence the 
Kl — Ks-mass difference is dominated by 771. This contribution has already been tabulated 
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in [17,61], so tfiat we comment on tfie corrections stemming from the second term involving 
773 here: Since 



to an accuracy of 2%, the term with 1]^ depends on \Vcb\ and p. Its average contribution 
to (Am)5/)/((Am)exp) is about 0.06. Hence one can obtain the full {Am)sD from the 
tabulated values for the first term in (7.16) in [17,61] by multiplying with 1.1. With the 
central value for rjf^'^ listed in (7.8) our short distance calculation therefore reproduces 
68% of the observed mass difference for rric — 1.3 and 73% of {Am)ex:p) for rric = 1.4, if the 
1/A^-result = 0.7 is used. Taking into account the large theoretical error in (7.8) the 
short distance contribution from (7.16) is between 52% and 95%. 

In most textbooks the Kl — Ks-mass difference is termed to be dominated by poorly cal- 
culable long-distance effects. Indeed, with the old LO result for rji and with the old smaller 
values for Aqqj-, the short distance part of Am is less than 50%. Long distance effects come 
from the operators with only light quarks in (5.28), which are also present below the charm 
threshold and may generate meson poles in the low-energy matrix elements. But a long 
distance dominance is clearly a puzzle, because by power counting these contributions are 
suppressed with Aqqy)/'>tiI with respect to the short distance part, because the coefficient 
of the leading dimension-six operator in (5.40) is proportional to (for a discussion see 
e.g. [68]). This fact has even stimulated speculations about new physics contributing to 
the Kl — Ks-mass difference (e.g. in [69]). With the new NLO result for rji the short dis- 
tance part is well above potential long distance contributions. Keeping in mind the large 
uncertainties in the problem it appears to be hopeless to find evidence of new physics in 
the Kl — Ks -mass difference. 



Re A; 




(7.17) 
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